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(N 

t^- | The validity of Weinberg's two sum rules for massless QCD, as well as the 

six additional sum rules introduced into chiral perturbation theory by Gasser 



X 



o 

Q\ ■ and Leutwyler, are investigated for the extended Nambu-Jona-Lasinio chi- 

ral model that includes vector and axial vector degrees of freedom. A de- 
Q^i tailed analysis of the vector, axial vector and coupled pion plus longitudinal 

Q_i! axial vector modes is given. We show that, under Pauli-Villars regulariza- 

tion of the meson polarization amplitudes that determine the spectral density 
functions, all of the sum rules involving inverse moments higher than zero 
are automatically obeyed by the model spectral densities. By contrast, the 
zero moment sum rules acquire a non-vanishing right hand side that is pro- 
portional to the quark condensate density of the non-perturbative ground- 
state. We use selected sum rules in conjunction with other calculations to 
obtain explicit expressions for the scale-independent coupling constants U 
of chiral perturbation theory in the combination Tj + ln(m^/;U 2 ), to evalu- 
ate the Das-Guralnik-Mathur-Low- Young current algebra expression for the 
electromagnetic mass difference of the pion, and to compute the pion elec- 
trical polar izability oie- The sum rule calculations set an upper limit of 
ole < ce/(87T 2 m n f%) w 6 x 10~ 4 fm 3 on this quantity, that is independent 
of the quark mass up to C(m 2 /m 2 ), and only depends on the fundamental 
constants: pion mass and pion decay constant /„-, in addition to fine 
structure constant a 1/137. 
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I. INTRODUCTION 



Our present knowledge of the low energy hadronic spectra and the properties thereof 
is based on the premise that chiral symmetry is a good approximate symmetry of the un- 
derlying quantum chromodynamics (QCD) Lagrangian. This has found verification time 
and again in understanding the nature of the low-lying octet of pseudoscalar mesons in 
particular. Calculational methods have been developed primarily in two complementary di- 
rections. On the one hand, the current algebra tools of the 1960's |1|] plus the two sum rules 
of Weinberg |||| attempt to address the properties of the QCD Lagrangian per se in the 
perturbative regime. On the other, the chiral perturbation theory (CHPT) approach 
which is based on an effective Lagrangian that only contains mesonic degrees of freedom, 
addresses the low-energy non-perturbative regime of QCD. The latter approach introduces 
terms into this Lagrangian in an expansion in momenta. The unknown coefficients of these 
terms in each order of the expansion, which in the two flavor case amount to ten in total, 
have to be fixed by experimentally measurable quantities. 

An effective model approach of describing the low energy sector ab initio in terms of quark 
degrees of freedom has also been studied for some years now in the so-called Nambu-Jona- 
Lasinio (NJL) model [pHTB||. In its original two flavor formulation, the minimal version of this 
model contains scalar and pseudoscalar fields in a chirally symmetric combination, and the 
model only has three parameters: the (strong) coupling G\, an averaged current quark mass 
m = (mJJ + m°)/2, and a cutoff A that sets the scale of the theory. Mesons are constructed as 
collective excitations within this model, and the Goldstone theorem is fulfilled in the chiral 
limit. The spontaneous breakdown of chiral symmetry leads to a dynamically generated 
quark mass. In addition, it can easily be demonstrated within the framework of this model 
that the Goldberger-Treiman and Gell-Mann-Oakes-Renner relations are fulfilled, and that 



the quark axial current is conserved. In addition, Dashen's theorem [13] that prevents 
a chiral neutral pion from acquiring a mass shift in the presence of an electromagnetic 
field, also continues to be obeyed []15j| . These properties, which are a consequence of chiral 
symmetry, continue to hold in an extended version of the NJL model or ENJL model, as 
we will abbreviate it fl8||L3l|. In this extended model, the vector and axial vector degrees of 
freedom are explicitly included, and an additional (strong) coupling strength G 2 is required 
in order to couple these modes into the Lagrangian. It is imperative in demonstrating that 
the expected properties of chiral symmetry hold, that a suitable regularization scheme be 



employed under which the relevant Ward identities continue to hold p5|JT6[] . 

The primary purpose of this paper is to investigate the sum rules of both the Weinberg 
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type || and also those introduced by Gasser and Leutwyler |J. These sum rules involve, 
inter alia the vector as well as axial vector spectral densities, and thus it appears necessary 
in regarding a quark model, to have such degrees of freedom explicitly supported. Hence 
our interest in using the ENJL model for this purpose. However, the application of this 
model to this problem is fraught with difficulties. On the one hand, the identification of 
the vector and axial vector meson modes with the physical p and a\ mesons is especially 
problematic, and we will return to comment on their analysis somewhat later on in this 
introduction. We turn first rather to the issue of the spectral functions themselves. We 
note that either version of the NJL model constitutes an effective field theory that does not 
confine the explicit quark degrees of freedom that it contains. Consequently one expects 
the strength distributions of the meson spectra, which contain interacting quark- ant iquark 
pairs, to behave unphysically in the high energy domain. This domain is naturally important 
if the sum rules are to work. However, in spite of these deficiencies, one can show explicitly 



17 1 that under the Pauli-Villars (PV) regularization of the amplitudes, Weinberg's first sum 



rule 

ds 



Jo s 

is obeyed exactly, while the second one satisfies a modified equation, 

ds(p^s) - pf(s)) = -m$1>) ^ . (1.2) 



Both statements hold to leading order in the number of colors. Here f n is the pion decay 
constant, and m and (ipip) are the quark mass and scalar quark condensate in the ground 
state. These conclusions also hold for the minimal version of the NJL model. This is perhaps 
all the more surprizing, since the minimal version does not support either vector mode as 
a fundamental excitation, so that it contains no resonance structures of any kind in these 
channels. 

Besides the two Weinberg sum rules, we are able to demonstrate that the six additional 
sum rules introduced by Gasser and Leutwyler are also obeyed by the spectral densities of 
the ENJL model. However, these results come at a price: as pointed out in fl8|| , for example, 
the Pauli-Villars regularization procedure effectively introduces additional spinor fields that 
may correspond to indefinite metric sectors of the Hilbert space. In the present case, the use 
of the PV procedure leads to negative and therefore unphysical spectral densities at large 
momentum transfers. However, this occurs in such a fashion as to give an overall result that 
is in accordance with the sum rule. We thus interpret the, perhaps unexpected, continued 
validity of the sum rules under these circumstances as a consequence of the underlying chiral 
symmetry and its faithful representation in the Pauli-Villars regularization method. 
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The results obtained from the ENJL model fall naturally under scrutiny with regard to 
their accordance or lack thereof with chiral perturbation theory. The retention of the sum 
rules in both approaches suggests that there could be considerable overlap in these results. 
In fact, the ENJL model provides a dynamical basis for calculating the empirical Gasser- 
Leutwyler coupling constants of chiral perturbation theory. We evaluate these constants 
in both the NJL and ENJL models, and express them as far as possible in terms of the 
physically measurable parameters f n and g^, the quark axial form factor. One finds that one 
additional physically non-measurable quantity, the (constituent) quark mass m, also occurs 
in these expressions, giving them this single parameter dependence. However, since the 
calculations of all quantities within the ENJL model have only been performed to leading 
order O(N ) in the number of colors, one has to identify the calculated coefficients with 
the Gasser-Leutwyler constants with some care. The reason for this is that the physically 
meaningful scale-independent constants U of Gasser and Leutwyler are not well-defined in 
the chiral limit, where they diverge logarithmically with vanishing pion mass. This comes 
about due to the appearance of chiral logarithms in CHPT that are to be anticipated on 
general grounds |19|| . In either version of the NJL model, such chiral logarithms appear as 
a part of the (1/N C ) non-chiral corrections that necessarily include mesons in intermediate 
states Thus we must identify the coefficients obtained from the model calculations 



with the renormalized /['s of Gasser and Leutwyler, that have a well-defined chiral limit, 
and not directly with the li as is commonly done (see PTJ and further references cited 



therein^]). We do this here. In this way, predictions can be obtained for the combinations 

k + Mml/v 2 ) ~ 11 (1.3) 

Here the mass /i is a scale parameter introduced via the dimensional regularization of the 
chiral perturbation theory amplitudes |J. One can explicitly validate this approach [|20 



the coefficients that have been evaluated to 0(N C ) in either version of the NJL model are 
finite in the chiral limit, and acquire chiral logarithms (at the scale fi = 2m « m a of the a 
meson mass) upon incorporating the meson loops. 

The results obtained in this manner are in reasonable agreement with the known empirical 
values of the Vs, after including the logarithm as shown. As already noted, this agreement, 
however, depends explicitly on the calculated value of the dynamically generated quark mass 



We note that the good agreement claimed in previous calculations can be traced back to the fact 
that the chiral logarithm is a slowly varying function. 
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as obtained from the gap equation. An important exception to this remark is the result for 
the electric polarizability, a^, of charged pions. This coefficient can be obtained directly 
from Holstein's sum rule in conjunction with the ENJL expressions for the spectral 
densities and the pion radius. One obtains the simple form 

487r 2 /^ 8n 2 f2 m w 8n 2 f 2 

where a = e 2 /47r pa 1/137, and the other symbols have their usual meanings. This is an 
extremely interesting result for several reasons. Firstly it coincides exactly in form with the 
CHPT expression for the polarizability in terms of the U to one- loop order ||23|| . Secondly, the 
terms coming from the dynamically generated quark mass m cancel exactly in the difference 
[Enjl _ ienjl^ i eac ij n g t a result that only contains physical constants. Thirdly, only 
the axial meson degree of freedom contributes to the polarizability through g^; there is no 
contribution from the p meson. Fourthly, since g\ < 1, this establishes a parameter free 
prediction for the upper limit on the ENJL value of ag, at least to 0(ctN c ). Using the 
estimated value = 0.75, one has the numerical result, ag = 3.40 x 10 _4 fm 3 . This answer 
can be compared directly with the prediction [f22j of chiral perturbation theory to one-loop 
order, af pt = 2.68 x 10" 4 fm 3 , that uses the empirical values of the Vs. Experimentally, the 
value of cve is still very uncertain. Current analysis of the available data gives values of Qg 
ranging between 2.2 ± 1.6 to 20 ± 12 x 10~ 4 fm 3 p^-H- 

One can also evaluate the current algebra result of Das et al. p9| for the electromagnetic 
splitting of the charged to neutral pion masses using the ENJL spectral densities. The Das 
et al. formula links the splitting AmJ = m 2 ± — m^ to a modified sum rule for the spectral 
density difference p\ — pf, integrated over all momentum transfers of the virtual photons 
that are responsible for the electromagnetic self-energy of the pion. The resulting ENJL 



expression for this mass squared difference can be written in the compact form [17 



AmJ = 3*e 2 / ^^(^^(^^(g 2 ) (1-5) 



(27r) 4 g 2 ' 

where the F's are pion and quark form factors that are defined in Section II of the main text. 
The special assumption that the form factors in Eq. ( |1.5| ) describe point-like mesons without 
internal structure, leads to the result originally derived by Das et al. Further evaluation 
of this expression is postponed to Section IV. There we also give the set of gauge invariant 
pion polarization diagrams to 0(aN c ) that corresponds to this result. The electromagnetic 
gauge invariance of the mass difference expression is thereby made explicit. 

As has been alluded to in an earlier part of the introduction, the determination and 
interpretation of the vector and axial vector excitation modes that are supported by the 
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ENJL model is not straightforward, but remain essential ingredients for the construction of 
the spectral densities. We have therefore investigated these modes anew. The particular 
innate difficulty in these calculations arises from the fact that both the p meson and the 
axial vector meson a\ have empirical masses that lie above the 4m 2 threshold for free quark 
pair creation, a non-physical process from the point of view of confinement. The existence 
of this threshold means that the complex plane for the momentum transfer variable squared 
s has a unitarity cut along 4m 2 < s < oo. Accordingly the ENJL model masses have to 
be sought as the poles of the relevant meson propagators that have been continued through 
this cut onto the contiguous or "second" Riemann sheet. 

Taken at face value, these facts seem to make a model of this nature of dubious value for 
the vector modes. In fact, the situation is far more complicated than this. A careful analysis 
of the vector polarization shows that this function has two isolated poles on the second sheet 
for our parameter choice: one on the real axis below 4m 2 , and a second complex pole in the 
lower half plane of that sheet. The real pole corresponds to a p mode containing virtually 
bound qq pairsQ. The p spectral density function receives contributions from both of these 
poles. A single pole approximation to describe the resonance behavior of the vector strength 
in terms of a single mass plus a width parameter is thus incomplete. While this density still 
displays a single-peak structure that can be placed around the observed physical p mass 
for an appropriate choice of parameters, the extent to which the virtual state strength 
intrudes^ into the energy domain of the peak is parameter-dependent. Thus one cannot 
unambiguously associate the position of the peak with the mass squared of the p meson as 
being given by either the virtual bound state pole, or the real part of the complex pole. We 
give full details of this problem in Section III. As such, the interpretation of the p meson, like 
the 7r meson, differs vastly from the conventional view of being simply a bound state of two 
constituents. Note in passing that, simply taken on its own, the p meson is of considerable 
interest currently, given its special role in the interpretations of the observed excess in the 



2 The occurrence of this state in the ENJL model was first pointed out in |kJ. Increasing the 
coupling strength G2 of the vector plus axial vector terms in the Lagrangian causes this mode to 
move through the branch cut onto the physical sheet, to become a true bound state for the model 
p meson. 

3 The enhancement of the neutron-proton scattering cross-section in the singlet channel over that 
in the triplet channel due to the presence of a virtually bound singlet state in the former, is a case 
in point [31]. 
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lepton spectra in heavy-ion collisions in the intermediate mass range |32| . 

By contrast, the axial vector polarization displays a single complex pole in the lower half 
plane of the second sheet, and the axial vector spectral density has a single peak structure 
that can be associated with this pole. The positions of peaks in both the vector and the 
axial vector spectral functions are basically controlled by the vector coupling strength G^- 
While the value of this coupling can be fixed via qa at G 2 — (1 — (?a)/(8/^) ~ 3.61GeV~ 2 
[from Eq. ( |2.85| )], plus the known value of /„- ~ 93MeV and an assumed value for = 0.75, 
the actual peak positions are fairly robust against moderate changes in G 2 - Taking the 
resulting peaks to define the p(770) and ai(1260) meson masses, one obtains m p = 713MeV, 
m ai = 1027MeV for the vector masses of the ENJL model, to be compared with the observed 
masses of 768.1 ± 0.6MeV and 1230 ± 40MeV It is interesting to observe that the 

predicted ratio m ai /m p = 1.44 rs \/2 for the above value of G2 is remarkably close to the 
original Weinberg estimate of exactly \/2. 

The values for the model masses are satisfactory from an experimental point of view. One 
thus has to try to reconcile these results with the fact that neither spectral function contains 
its dominant physical decay channel, p — > hit or a\ — > pre respectively, but rather describes 
the (hypothetical) p — > qq or a\ — > qq decays allowed by the non-confining ENJL Lagrangian. 
This is a long-recognised deficiency ||. On the other hand, the model spectral densities 



continue to saturate the first sum rule, Eq. (14.). From these facts alone, one can argue with 
some justification that the positions of the peaks in these distributions are more significant 
than the particular manner in which they decay, and therefore do lead to meaningful masses. 
One should recall in this regard that, in Weinberg's original discussion of the problem, the 
vector meson mass ratio was determined by using delta function distributions for the vector 
densities that ignored all decay channels, but with strengths that allowed the sum rules to 
be satisfied. More generally, the fact that both Weinberg's and Gasser and Leutwyler's sum 
rules are satisfied by the model spectral density functions in spite of their undesirable decay 
properties, suggests that results based on integrals over these densities do nevertheless have 
physical significance. 

An analysis of the ENJL model commences with the construction of the irreducible polar- 
izations that are now coupled in the isovector pseudoscalar and longitudinal axial channels. 
This is done in Section II. We are not the first authors to carry out such an analysis (see e.g. 



T3|j). Hence this section is in the nature of a summary of those aspects of the problem that 



are necessary for an understanding of the sum rules, and it also serves to set our notation. 
From a knowledge of the irreducible polarization functions, the new meson propagators in 
the extended model can be constructed. This is done in such a fashion that the generaliza- 
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tion from the NJL model becomes obvious, and so that the NJL results are simply recovered 
on setting G 2 to zero. Central to the derivations is that, as far as possible, quantities are 
expressed in terms of physical variables, the pion decay constants f n and the pion quark 
coupling strength, g nqq , in the ENJL model. One observes that these quantities undergo 
a simple multiplicative scaling f n = y/gXfp and g nqq = g p / ' s/Qa via the quark axial form 
factor g^ with respect to their NJL values f p and g p in such a way as to preserve the chiral 
theorems, like the Goldberger-Treiman relation. In fact, apart from such scaling, the main 
effect of including the vector and axial vector degrees of freedom is to drastically redistribute 
the associated NJL spectral densities into two peaks that serve to determine the masses of 
these modes in the model. Conservation of the axial vector current can also be demonstrated 
explicitly. Here one finds that the ENJL model leads to non-trivial vector and axial form 
factors F v (q 2 ) and Ga(o 2 ) for the quark currents such that F v (0) = 1, but Ga(0) = gA < 1- 
There is of couse an overriding caveat on all of these statements: the issue of regular- 
ization. In addition to being non-confining, the NJL model and its ENJL extension are 
both non-renormalizable. Thus a cutoff A in their mass spectra is necessary in order to 
obtain finite answers. One anticipates that the magnitude of this cutoff will be similar for 
either version of the model, but this parameter has in any event to be of order of the chiral 
symmetry-breaking scale A x ~ lGeV of QCD. Otherwise the NJL Lagrangian could not 
reasonably fulfill the role of a low energy effective model theory. However the problem is 
more complicated than this. One has to also ensure that the regularization scheme em- 
ployed respects the chiral Ward identities satisfied by those amplitudes that express the 
consequences of the chiral symmetry of this model Lagrangian. The Pauli-Villars regular- 



ization scheme 0] fulfulls these requirements [|15| , |T6f : in particular, it preserves those chiral 
Ward identities between regulated amplitudes that are necessary for a demonstration of the 
sum rulesf]. We take up the PV regularization of the polarization amplitudes in some detail 
in Appendix A. There we show that the PV scheme can be implemented in conjunction with 
either unsubtracted, or once-subtracted dispersion relations that lead, not surprizingly, to a 
very different asymptotic behavior of the polarization amplitude in question. This feature 
has a direct bearing on the ability of the associated spectral density, which is given by the 
imaginary part of the polarization, to satisfy the sum rules. 



4 This is not to say that the Pauli-Villars regularization scheme will guarantee all Ward identities. 
A discussion of the related questions that arise in the case of the PV regularizarion of the a-model, 
for example, can be found in |||]. 
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The spectral density functions carrying specific quantum numbers are next introduced 
in Section III with the aid of the appropriate current-current correlators and the Kallen- 
Lehman spectral representation. We establish the analytic properties of the polarization 
functions. Then we analyse the analytic behavior of their corresponding spectral density 
functions under PV-regularization, and show by means of contour integration techniques 
that all eight sum rules, the two of Weinberg, and the six of Gasser and Leutwyler, are 
obeyed by ENJL spectral densities regulated in this manner. The only exceptions that 
occur are for the zero-moment sum rules, like in Eq. (|1.2p , where a finite value instead of 
zero appears on the right hand side. We use these results in conjunction with the non-chiral 
expansions for the physical pion mass and decay constant of Section II to identify several of 
the coupling parameters of chiral perturbation theory. In parallel with these developments, 
we also explore the analytic properties of polarization amplitudes that have been obtained 
from once-subtracted dispersion relations prior to regularization. Such amplitudes are shown 
to give rise to Landau ghost poles at space-like q 2 in the n, p, and a\ channels. Some of the 
consequences of nevertheless admitting such unphysical propagators are briefly discussed. 

In Section IV, we discuss the use of ENJL spectral densities in conjunction with sum rules 
to compute two electromagnetic properties of the pion: the electric polarizability coefficient 
of the charged pion from Holstein's sum rule [^], and the — tt° mass splitting from the 



current algebra sum rule of Das, Guralnik, Mathur, Low and Young [£{J. We also provide 
an alternative derivation of the mass splitting in terms of Feynman diagrams that leads to 
exactly the same result. An explicit demonstration of how Dashen's theorem is satisfied by 
the ENJL model electromagnetic self-energies for the neutral pion is also given. A summary 
of results and conclusions is presented in Section V, while Appendix A contains additional 
calculational details of a more technical nature that do not appear in the main text. 

II. REVIEW OF FORMALISM 

Since the properties of Nambu-Jona-Lasinio type Lagrangians have already been the 
subject of extensive analysis in the recent literature |7|-|T3|, we use this section to briefly 



summarize only those aspects that will be required as background for the specific develop- 
ments to follow in Sections III and IV. 
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A. Two-point functions 



In this section, we discuss the two-point correlation functions from which the spectral 
densities we wish to study are obtained. Since the calculational procedure is the same for all 
spectral densities, let us denote by J a (x) a generic operator with SU{2) isospin and Lorentz 
indices summarized by the index a. Then the associated time-ordered two-point correlators, 
or polarization functions, are given by 

d 4 x e iq - x (0\T{J a (x) J fe (0)}|0) (2.1) 

in momentum space, where T is the usual time-ordering operator. These currents are as- 
sumed to be bilinear forms of the two-flavor quark fields i/j(x) and their Dirac adjoints, 
i.e. 

J a (x) = i){x)V a ijj{x), (2.2) 

where the vertex T a characterizes the nature of the excitation. Using the notation of Ref . [H , 
the bilinear forms P, S, V, A for isospin currents J a are 



+ 



TT 



P a (x) = ip(x)i'j 5 T a -ip(x) 0~ 

S a {x) = ^{x)r a ^{x) 0++ a 

V£{x) = ^(xh^^x) l -- p 

Al(x) = ^(x)^ 5 ^(x) 1++ a x 

with r a , a = 1 2 3 being the 577(2) Pauli isospin matrices. The spin, parity and charge 
conjugation quantum numbers J PC of the excitations they lead to relative to the invariant 
vacuum, as well as representative low-lying meson states that carry these quantum numbers, 
are also indicated. In addition to these isospin currents, it is useful to consider the isoscalar 
densities 

P°(x) = tp(x)i^(x) 0~+ r] ^ ^ 

S°(x) = ^(x)i(}(x) ++ fo or a 

The factor of 1/2 that is included in Eq. (|2.3|) in the definition of the vector and axial vector 
currents is purely convention. If we further assume that the dynamics of these quark currents 
are governed by a two-flavor £7(1) X SUl(2) x SUr(2) chirally symmetric Lagrangian of the 
Nambu-Jona-Lasinio type, then both the vector and the axial vector (Noether) currents will 
be conserved. If this is the case, then the Lorentz and and isospin structure of the expression 
Tl) is known ahead of time. If we write 
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/ 



d 4 x e^(0\T{r(x)Jtm\0) = -\^IMI (2-5) 



where J£(x) is given by either or A a ^ then both correlation functions become purely 
transverse tensors in their Lorentz indices and diagonal in their isospin indices, i.e. 



U^ ab (q 2 )=U J - J (q 2 )T,J ab , (2.6) 

where = (g ^ — q^q u ) with q^ = q^j \[q 1 . In a similar fashion, the scalar and pseudoscalar 
two-point functions have the structure 

J d A x e^ x (0\T{J a (x) J b (0)}\0) = -itl J \q 2 )5 ah) (2.7) 

where J a (x) is given by either S a or P a for the isovector currents, or by S* or P° for the 
isoscalar densities. In the latter case, there is no isospin factor. Note that in Eqs. (|2.5|) 
to ( |2.7| ), we have introduced the full polarization functions Il JJ (q 2 ). In order to achieve a 
uniformity of definition for these, we have introduced a factor 1/4 in Eq. ( f2.5| ) to compensate 
for the factor 1/2 carried in Eq. (|2.3Q in the current definition. 

It is important to bear in mind that the transverse form given by Eq. ( |2.6j ) holds in the 
first place for the exact vector or axial vector polarization in a chirally symmetric theory, and 
therefore must also hold for any approximate calculation of either function. Therefore any 
approximate treatment has to be implemented in such a manner that current conservation 
is maintained at each level of approximation: otherwise the physical consequences of the 
symmetries contained in the Lagrangian could become distorted. A case in point that is 
discussed again below is the mixing of the longitudinal axial vector and pseudoscalar pion 
modes already at the lowest order of approximation in the minimal NJL model that is 



necessary to ensure axial current conservation in the chiral limit |21]. 



B. Effective chiral Lagrangian 

In order to calculate the two-point correlation functions, we need to specify the La- 
grangian £ governing their dynamics. For this purpose, we identify £ with the two-flavor 
Nambu-Jona-Lasinio Lagrangian, that has been extended to include vector and axial degrees 



of freedom in a chirally invariant fashion J8, 12 , 13] . We write 



£enjl = $[$-rh]if> + Gi (H>) 2 + (V'iTsT^) 2 - G 2 {^^f + 0/^75 r</0 2 (2.8) 
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where the isovector and axial vector terms have to appear with a common coupling constant 
G2 in order to retain chiral invariance 5 . In this expression, the ip's are iso-doublet, color- 
triplet Dirac quark fields and rh is the averaged current quark mass m = (mP u + m°)/2. We 
also note that £enjl breaks the £7^(1) symmetry in maximal fashion due to the 't Hooft 



term [3l| that it contains [10 and which models the QCD instanton interaction. The G's 
should be assumed to scale like N^ 1 if the correct QCD scaling properties for the quark and 
meson masses are to be maintained. 

In the following, it will be convenient to refer to results that are obtained with G2 either 
present or absent as pertaining to the ENJL or NJL model respectively. Then, apart from 
the pionic and sigma-like excitations with J PC = h and ++ that are collective modes of 
the NJL Lagrangian built out of coherent superpositions of strongly interacting qq pairs of 
binding energy 2m and zero respectively 0, the ENJL Lagrangian also supports additional 
collective modes with J PC = 1 and 1 ++ that coincide with the quantum numbers of the 
p and a\ meson. Although we will continue to refer to these modes as "meson" states of the 
corresponding name for brevity, to what extent this identification with the actual physical 
mesons of the same quantum numbers is appropriate, is a matter for later discussion. 



C. Schwinger-Dyson amplitudes to leading order 0{N C ) via diagram summation: the 

RPA approximation 

Due to the specific structure of the Lagrangian in Eq. (|2.8|) , the Schwinger-Dyson [SD] 
integral equations for the full polarization functions can be solved explicitly J7| in terms of 
the lowest order single quark loop irreducible polarization diagrams. These individual loops 
are denoted as H JJ ' together with the appropriate additional spinor and flavor structure. 
They may be decomposed according to their Lorentz structure as 

II^V^cU (2.9a) 
[U^ A (q 2 )T^ + Yi A L A {q 2 )L^}5 ab (2.9b) 



5 One could equally well include additional chirally invariant terms in this Lagrangian, like the 
separately chirally invariant isoscalar vector and the axial vector terms (carrying the quantum 
numbers of the u> and /1 mesons). However, since the important collective low- lying excitations 
contained in £njl are the pions which are also isovector, the isovector vector mesons are the 
relevant degrees of freedom for our purposes. 



12 



IC(<? 2 ) = H AP (q 2 W ab , & = (2.9c) 

<V) = n^(g 2 )5 afe (2.9d) 

U s J(q 2 ) = U ss (q 2 )8 ab . (2.9e) 

The two remaining polarizations Il PoPo (g 2 ) and II °(g 2 ) generated by the density operators 
P° and S* are scalars in both Lorentz and isospin space. Here = is a longitudinal 
tensor such that T^ IU + L^ IU = g^ iu . As has been indicated explicitly, the irreducible po- 
larizations are all functions of q 2 . The full polarization functions that were introduced in 
Eqs. (|2.5|) to (|2.7[) and are denoted by Tl JJ ' have the same invariant structure as the irre- 
ducible one-loop amplitudes that generate them. Note that the vertices that these one-loop 
diagrams connect are indicated by the following superscripts, 

V^^^Ta, A^7 M 7 5 r a , P -> ry 5 r a , S -> r a , P «T5, So -> 1. (2.10) 

All such one-loop diagrams involve traces over all quark variables, which make them 
proportional to iV c . Thus their products with G\ or G 2 are independent of N c . Hence 
the RPA approximation for the polarization amplitudes, which involves a selective summa- 
tion of diagrams to all orders in the (strong) interaction parameters G\ and G2, will only 
contain interaction effects to leading order 0(1) in the expansion in powers of the inverse 
color variable 1/N C if, as here, we replace the exact irreducible polarization by its one- loop 
approximation. 

We can now write down separate SD integral equations for the pure vector (p), pure 
transverse axial vector, and coupled longitudinal axial vector-pseudoscalar (at, 71") channels. 
We can do this for either the scattering amplitudes in these channels or the polarization 
diagrams. If we do the latter, then one sees from the structure of £njl that the full vector 
polarization is simply given by the geometric progression of one-loop irreducible amplitudes 
connected via the interaction in this channel, 

nZ v = n%y - 2GM v uV = — u r v (2 n) 

T T T T i + 2G 2 n^' v ; 

and is automatically transverse and strictly 0(N C ) since the denominator is 0(1) in the 
number of colors. We refer to this in the following as the random phase approximation, 
or RPA because of its similarity with the diagram summation technique used in discussing 
collective modes in nuclear physics. The axial vector and pseudoscalar polarization functions 
are by contrast coupled by the amplitude n^ fc that allows for transitions between the 
longitudinal part of the axial vector and the pseudoscalar pion modes. This feature was 
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first observed in the early field theoretic formulations of Weinberg's model |36[ where it led 
to a finite pion field renormalization. The situation here is very similar although there are 
no elementary meson fields present in the Lagrangian. Taking this off-diagonal polarization 
amplitude into account, one has 

K A = K A + n#(-2G 3 )^n# + nf (2G?o^, (2.12) 

after dropping the now redundant isospin labels. Note that U AP I1 PA = Yl AP Yl PA q^q v is 
longitudinal. Hence one can split the solution of the above equation into parts that are 
purely transverse, and parts that are purely longitudinal: 

T\ AA 

U AA = U AA - 2G 2 U AA U AA = i - 2 g (2.13) 

and 

fi AA = u AA - 2G 2 n AA n AA + 2G 1 n AP fi PA . (2.14) 

The pseudoscalar polarization function can likewise be found as the solution of the pair 
of coupled SD equations, 

n pp = n pp + n pp (2G x )n pp + n PA (-2G 2 )n AP , (2.15) 



while 



fl P % = Tl p % + n PA f(-2G 2 )fl AA + U pp (2G 1 )fl p \ 



or 



fl PA = u PA + Tl PA (-2G 2 )fl AA + n pp (2G 1 )fl PA , (2.16) 

since U PA ' U I1 AA = 11^11^^ projects out onto the longitudinal part of the axial polariza- 
tion. Thus we end up with a single equation for the transverse axial polarization, but three 
coupled equations ( |2.14| - |2.l5| ) for the pseudoscalar polarization and the longitudinal part of 



the axial polarization. These latter equations can be written as a single matrix equation in 
the P, A space as 

n = n + nKri, (2.17) 

where the boldface symbols have the matrix structure, 

, n pp U PA \ ( 2G X , 

n= „ , K= . 2.18 
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The II has the same matrix structure as II. While the solution to this matrix equation may 
be found without difficulty |J, it is more instructive to proceed as follows. We split up the 
matrix of interaction strengths as 



K = K! + K 2 = + | , (2.19) 




and introduce the auxiliary matrix of polarization functions II, distinguished by carrying a 
hat, that are generated by the axial vector part of the interaction acting in isolation. This 
auxiliary matrix satisfies the equation 

n = n + riK 2 n, (2.20) 

with solution 

n = (1 - nK 2 ) x n. (2.21) 

Subtracting the matrix equations satisfied by 11 and A respectively and using the solution 
for II given above, one finds that 

n = n + nK 1 n, (2.22) 

with solution 

n = (l-nK!) ifr (2.23) 

This result shows that the role of the longitudinal axial mode can be completely absorbed in 
the determination of the auxiliary polarization matrix II. The effect of this is to renormalize 
the one-loop polarization functions as follows, as may be seen by evaluating Eq. Q2.21|) : 



yrPAyrAP 

Y\ AP 

n AP = -t^ttaa (2-25) 



1 + 2G 2 n 

tjAA 
t AA ll l 



L 



ni = irtsF' (2 ' 26) 

where the argument q 2 has been suppressed for brevity. Physically this renormalization 
comes about due to the partial summation to infinity of one-loop diagrams connected by the 
longitudinal piece of the axial vector interaction. Together with the explicit solution for the 
transverse axial polarization given by Eq. (|2.13|) above, Eqs. Q2.21 ) and (|2.23|) provide us with 
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a complete description of the polarization matrix. It will turn out that the longitudinal part 
of vanishes identically in the presence of the pseudoscalar interaction term contained 
in £njl (in the chiral limit). This is not true of Et^f that continues to have a longitudinal 
component. 

Of the remaining polarization functions, only the isoscalar one, Yl SoS ° is modified by the 
interactions contained in £njl- Using the RPA again, one finds 

All that is required now are the explicit forms of the matrix elements of the polarization 
at the one loop level to complete the solution for the various polarization functions. These 
are given below. 



D. Calculation of the one-loop polarization functions 

The one-loop diagrams corresponding to the isospin currents of Eq. ( p.3[ ) are given by 



(2.28a 

= * J ^Tr[ llx r a S(p + q) lv r h S{p)\ (2.28b 

K A ,ab = * J -0j-Jr[ W r a S(p + qh^ 5 r b S(p)] (2.28c 

K P ab = * J J^kTr[w5T a S(p + q)i^r b S(p)} (2.28d 

nfif = i I -0^Tr[t 75 r a S(p + q)il 5 r b S(p)} (2.28e 

nff = i J 4^ T rKS(p + q)r b S(p)}, (2.28f 

where the trace Tr operation runs over all intrinsic variables, and where S(p) = (p , —m+ie)~ 
is the quark propagator. The corresponding isoscalar diagrams are 

n P ° P ° = i J -0-Tr[i l5 S(p + q)i l5 S(p)} (2.29a 

U s ° s °= i f AV-Tr[S{p + q)S{p)\. (2.29b 



We first consider the chiral limit, rh = 0, and start with the scalar and pseudoscalar polar- 
ization diagrams, as well as the (h^l) coupled diagram. It has been shown [p|-|TT|j that all 
these diagrams may be written in terms of the mean field (or Hartree) expressions for the 
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pion weak decay constant f p , the irqq coupling constant g p and pion form factor Fp in the 
chiral limit. One hasQ 

n pp = (2d)- 1 + S^q 2 = Tl p ° p ° (2.30a) 
m 2 

n AP = 2z/ p 2 (g 2 ),/4 = (n^)* (2.30b) 

' V J71 



n SoSo = + ^|^(g 2 - 4m 2 ) = U ss , (2.30c) 

after introducing the abbreviation 

f 2 p (q 2 ) = f 2 p F P (q 2 ). (2.31) 



Here 



/ 2 = -4ziV c m 2 /(0) (2.32a) 
^ 2 = [-^/(O)]" 1 (2.32b) 
F P (g 2 )=/(g 2 )//(0), (2.32c) 



where I(q ) is the common loop integral 



l^h = I tSiTTZ I 35 LiO ^T - (2-33) 



(27r) 4 [(p + g) 2 — m 2 ] [p 2 — m 2 ] 

A suitable regularization of this divergent integral is discussed in Appendix A. Notice that 
the mean field pion decay constant and pion-quark coupling satisfy the Goldberger-Treiman 
relation f p g p = m at the quark level, where the quark mass m arises from spontaneous chiral 
symmetry breaking in the groundstate of £njl- This mass is self-consistently determined by 
the single quark self-energy £ as m = S(m), the so-called "gap" equation JFj. In the Hartree 
mean-field approximation, one finds for either the NJL model or its ENJL extension that 

d 4 p 



m*Z H = 2G 1 J *LTr[iS(p)] 



(2tt)< 

lQG^m [ f P - ' (2.34) 
J (27t) 4 p l — m l + ie 



6 The expressions for the scalar and pseudoscalar amplitudes also contain surface terms |}7|] that 
arise because of the momentum variable shifts that have to be introduced to simplify them us- 
ing the quadratically divergent gap equation. The surface terms are removed after Pauli-Villars 



regularization as discussed in more detail in Ref. [15|. 
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since only the scalar interaction in £njl is operative in determining This equation leads 

to the relations IT pp (0) = n 5oSo (4m 2 ) = (2d)' 1 that have been used to simplify II PP and 
U s s in Eqs ( p^a|) and (g^D§). 



We now turn to the vector and axial vector polarization diagrams of Eqs. Q2.28b| ) and 



( p.28c| ). Commuting jsS(p) = S(— p)^ under the trace in the latter expression leads to 



K»M) - K»M) = 4f 2 p F P (q 2 ) gil J ab . (2.35) 

One also notes that n^Jf a6 is purely transverse, as may be established by contracting with 
and using the identity </= S~ 1 (p + q) — S~ 1 (p). Then 

<fKU = * J 704^r [S(p) lu - S(p + q) lv ] 5 ab = 0, (2.36) 

provided that the necessary shift in the momentum integration variable in the divergent 
integral on the right is allowed. If so, one can conclude that 

n£V) - nr(<? 2 ) = nf A (g 2 ) = 4/ p V), (2.37) 

by inserting the axial and vector polarizations in terms of their transverse and longitudinal 
parts from Eqs. (|2.9b|) and (|2.9a| ), projecting out with the help of the relations T^T^ 



T^g^ = 3 and L^L^ = 1, and using Eq. (|2~3T 



The relations (|2.37|) , which are in the nature of Ward identities, are crucial for proving 
the Weinberg sum rules as well as establishing the relation between the Feynman diagram 
structure and the Das et al. formulation for the electromagnetic pion mass splitting in 
the ENJL model. Since all these quantities involve formally divergent integrals at the one- 
loop level, their regularization has to be implemented in such a fashion that Eqs. ( |2.35| ), 



( |2.36| ) and Eq. ( |2.37|) on which these relations are based, remain valid between the regulated 



quantities as well. We discuss this later. The explicit form of Uj> v {q 2 ) after regularization 
is given in the appendix. It is shown there that n^ y (0) = 0. Thus one sees from Eq. ( |2.37| ) 
that 



Ur (0) = nr (0) = 4/;. (2.38) 

The ingredients needed for constructing the polarization matrix II are now available: 
the pseudoscalar and mixed polarization diagrams appear in Eqs. ( |2.30a| , |2.30b| ), and nf A is 
already known from Eq. ( |2.37| ). The resulting one- loop polarization in (P, A) space is given 
by the hermitian matrix 

n= t (2G 1 )- 1 + / P 2 (g 2 )i -2if}tf)J£\ 
{ ^(g 2 )^ 4/ p V) )' 
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After renormalization by the longitudinal axial mode according to Eq. (|2.21|) , this becomes 
ii | " rV '"' ~ " ,! : ' ' V "'" I ■ (2.40) 



(2G 1 )- 1 + /. 2 (g 2 )^ -2^V)V5 



Notice that the so/e effect of the longitudinal mode is to replace the factor f p (q 2 ) of Eq. (|2.31|) 
by 



/;(<f ) = #w 



(2.41) 



It is shown later [Eq. ( |2.79| )1 that the chiral limit q 2 — > of this expression gives the pion 
weak decay constant f n for the ENJL model. Since 



riKi) 



where the determinant D is 



v- 1 



1 



(2.42) 



P = £>et(l - nKi) = 1 - 2G l U pp (q 2 ^ 

the matrix elements of II = (1 — nKi) _1 II now follow as 

fl pp (a 2 ) 



(2.43) 



n 



pp 



1 - 2G 1 U pp {q 2 



iif 2 M 2 " 



n 



AP 



1 - 2G 1 U pp {q 2 ) 

u PA = {fi AP y 



where 



fir = o, 



n pp (g 2 ) = (2G0 



-i 



(2.44) 

(2.45) 

(2.46) 
(2.47) 

(2.48) 



is the renormalized version of the pseudoscalar polarization of ( p.30a| ) . Thus II has a zero 
in the lower right hand corner: the full longitudinal polarization vanishes as it must, to 
ensure axial current conservation when the longitudinal axial vector (oll) mode mixing with 
7r is properly taken into account. Note that this conclusion also remains true for the minimal 
NJL model, where G 2 = . 
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E. Propagators 



Together with Eqs. ( |2.11| ), ( [2.13| ) and (|2.27|) , the above suite of equations completes the 
solutions for the transverse vector and axial vector, scalar, and the coupled pseudoscalar- 
longitudinal axial vector polarization functions for the ENJL model in the RPA approxima- 
tion. We now turn to the task of constructing the associated propagators for these various 
modes. 



1. The p and transverse a\ mode propagators 

The scattering amplitude of a quark-antiquark pair arising from the exchange of either 
vector meson mode has the form 

(IT)i 2 \-DtUl 2 )} ( T f U* ( 2 - 49 ) 

where 

r M / = 7 V or 7 M 7 5 r a , (2.50) 

for the vector or axial vector channels respectively. The D^. ab = D J J,8 a b defines the corre- 
sponding meson propagator. In the RPA approximation, it is built up from the interaction 
lines —2iG 2 and 2iG\ that connect 0, 1, 2,.... polarization bubbles to the initial and final 
vertices. In the vector channel, the sequence of contributions for —iDYJf reads 

- *D™ = (-2iG 2 )g^ + (-2iG 2 ){ll^ v T lw + (-2G 2 )(I^ v )% p T^ + ■■■ }(-2iG 2 ) 

since the irreducible vector polarization is purely transverse. Hence the transverse and 
longitudinal parts of the rho propagator are 



The corresponding expression for the axial vector propagator is complicated by the presence 
of the coupling amplitudes II^ P . In place of Eq. ( p.51| ), one has the series 

- iDtf = {-2iG 2 ) 9lxv + (-2iG a )i{n^ + n^(-2G 2 )^'n^ + • • • 

+U^ p (2G 1 )U^ A + ...}(-2iG 2 ). (2.53) 
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The series in the curly brackets is just the series for the full axial polarization given 



in Eq. Q2.12 ). However, as this is again purely transverse, see Eq. ( |2.47|) , the results for 
the transverse and longitudinal parts of —iD AA are, according to Eq. ( p. 13 ) , identical with 
those of -i£>J7 with Wf v replaced by IL$ A . Thus 



= 1 { +2G%^ ~ lDtA = { ~ 2tG2) - (2 ' 54) 



The longitudinal part of this propagator also has non-vanishing off-diagonal elements, as we 
discuss below. 



2. The coupled ir — mode propagators 

As we have seen, the pion and longitudinal axial modes 7r and are coupled through 
the off-diagonal elements in EL The corresponding SD equation for the resulting propagator 
is thus a matrix equation involving the full polarization matrix EI, 

% i % 



iK + *K-{n + Kn + KnKn+ }«K 



= iK + zK~fiiK, (2.55) 
i 

which in view of Eq. ( |2.23| ) has the solution 

-iD = zK + zK(l - IIK^riK. (2.56) 
Insert the known value of the inverse matrix to obtain 

- !D -*™(-«c)> (2 ' 57) 

where A, B, C are the real matrix elements 

A = 1 (2.58) 
B = 4G2fitf)j£ (2.59) 

C = 2G 2 f*(q 2 )^, (2.60) 
and where the renormalized propagator for a chiral pion, 
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has been factored out. Here 



+ SGWpV)] 5 ■ (2-62) 

As shown below, the chiral limit of this function gives the renormalized nqq coupling con- 
stant g nqq (0) = g^qq. The auxiliary constants f p and g p have thus disappeared from the 
problem, to be replaced by their renormalized values f w and g nqq . Although this will be 
justified later by examining the conserved axial current, one already sees that these renor- 
malized values also satisfy the GT relation f^g-nqq = m. 

The matrix D enters into the qq scattering amplitude and has to be diagonalized in 
order to identify the independent tt and longitudinal axial eigenmodes and their coupling 
constants to the quarks. To this end, we introduce the real symmetric matrix 

. A B 

M — (2.63) 



B C 



that has real eigenvalues 



K,a L (<l 2 ) = + C) ± -y/(A - Cf + AB\ (2.64) 

where the signs correspond to the identification of the re and eigenmodes according to 
the expected limits A^- — > 1, \ aL — > as G 2 — ► 0. The orthogonal transformation matrix 
that renders M diagonal is 

/ cos 9 - sin 9 \ 2B , 

U= , tan2fl= — — . 2.65 

\sin£ costf ) (A-C) K } 

The scattering amplitude can then be written as follows 

/ :(/-,vr.-/ #-,vr)M| \ l) r 

i #7 5 r 



pp 



where 



= (^75r , -i # 7 5T)U{U T MU}U T \ I 5 }d 

V 1 ^ T J 

= T 7T + T aL , (2.66) 



T n = [(cos 9- #sin #)i 75 r] <g> [(cos 9+ #sin 9)i l5 r] X n D (q ) (2.67) 



and 
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T aL = [(sm6+ $cos6)i l5 T} <g> [(em 9- #cos^)i 75 r] X aL D (q ). (2.68) 

The separate contributions of the propagating pion and longitudinal axial eigenmodes are 
now self-evident. However, since the eigenvalue \ aL ~ q 2 as q 2 — > according to Eq. ( j2.64| ). 



the product X aL D pp (q 2 ) does not have a pole at q 2 = 0. Thus the longitudinal axial scatter- 
ing amplitude does not represent a propagating physical particle carrying mass. Nevertheless 
this amplitude is an essential adjunct to the propagating pion eigenmode, which does repre- 
sent a physical particle, as far as conservation of the isospin current is concerned, as well as 
the interaction of the pion with electromagnetic fields. These comments will be elucidated 
in the later sections. By contrast, from the structure of D pp (q 2 ), it is evident that the pion 
propagator part of T does have a pole at q 2 = 0: the Goldstone mode has been preserved 
by the RPA. However, as already noted before (see e.g. [§]), the coupling constant of the 
quarks to the (Goldstone) pion as well as its weak decay constant are renormalized, and 
an additional pseudovector coupling constant, g pV} is induced by the it — a L coupling in the 
ENJL Lagrangian. One confirms this by examining the residue of T w at the pole q 2 —>■ 0. In 
this limit, one finds from Eqs. Q2.64|) and (|2.65|) that 

K -> 1 (2.69) 
cos#^l (2.70) 

AgA (2.71) 



y^q 2 m 



Then can be recast as 



T, - [(g^- rff^Tfcr] ® i(g mq + ^)i7fer]i (2.72) 



in the vicinity of the pion pole, that allows one to identify the pseudoscalar and induced 
pseudovector pion-quark coupling constants g nqq and g pv as 



and 



g nqq = lim <^ 9 (gVW)c°s0 = g p [l + 8G 2 f 2 }^ = (2.73) 

? J 9 A 



sin 9 



g pv = lim 2mg 7Tqq (q 2 )J\(q 2 ) —= = g vqq (8G 2 fl) = g* qq (l - g A )- (2.74) 
q 2 ^0 ^jq^ 



The quark axial form factor g^ that is identified in Eq. (|2.85 ) to follow, has been introduced 



to simplify the final expressions for these two coupling constants. 
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F. Chiral symmetry invariance to 0(N C ) in the RPA 

Apart from satisfying Goldstone's theorem in the RPA, the ENJL model also obeys the 
remaining chiral theorems as realized by the Goldberger-Treiman relation and axial current 
conservation. The Goldberger-Treiman relation can be established by calculating the one 
pion to vacuum matrix element of the axial current that defines the pion weak decay constant 
f n . In the quark model, this matrix element is determined by connecting the vertex for the 
incoming pion, in this case i{g-n qq — 4 g pv /2m)i , y5T c/ , via a quark loop to the axial current 
vertex 7^75 r a /2. Since the induced pseudovector term in the nqq interaction only couples to 
the longitudinal part of the axial polarization, the resulting single quark loop diagram, can 
be re-expressed in terms of the U AP and Ii A A . The details are summarized below. One has 



But 



n ^(? 2 ) = 2i ^- and KU^W = 4f 2 p (q 2 )q»5 ab , (2.76) 



in view of Eqs. (|2.30b|) and ( 2.37Q . Hence the pion decay constant satisfies the relation 



P 

U = —(gnqq- 9pv) ■ (2-77) 

m 

However 

(g-nqq ~ 9pv) = ^TrggfiU (2.78) 



follows from Eqs. ( 2.73 ) and ( 2.74Q . Inserting this difference into the expression for /„., one 
finds 



U = Vuf P = f P [l + 8G 2 / p r 3. (2.79) 
Thus fa and g^ are again connected by a Goldberger-Treiman relation, 

Ugnqq = ^/g~Af P -^= = m. (2.80) 



9A 

We also establish the axial current conservation for the Lagrangian given by Eq. (p.8[) 
and identify the quark isovector axial form factor. The expected modification of the quark 
current due to the exchange of transverse axial ai plus the tt and longitudinal axial 
eigenmodes is depicted in Fig. [TJ Since the current operator represented by the sum of these 
diagrams is to eventually be sandwiched between initial and final quark momentum states p 
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and p' where q = p' — p is the momentum entering the vertex, we can freely use the identity 
^7 5 = 2m75. Translating the Feynman diagrams in Fig. |1| with the help of Eqs. ( f2.76| ), 



|2.67| ), ( |2.t)8D and ( 2.54j ) in terms of the polarization loops and propagators they represent, 



one finds, after some calculation, that 

J»5 = 7M7sy-— —%^7-^{ X Acose-y= sm6) + A o (sm0 + -= cos0) }- 



q 2 PM 2 ) ' V vV ' J 2 

7 ^7^75-. (2.81) 



l + 2G 2 n^ A4 " 7 75 2 



The contents of the curly bracket is given by 

4m 2 „ 4m 



{ • • • } = A + ^T C - = 1 - 8G ^V) = T~^hr ( 2V (2-82) 



see Eq. ( |A8(j| ) of Appendix A. This just cancels out the factor fp(q 2 )/f 2 (q 2 ) in view of 
Eq. ( p. 41 ) and we are left with 



/ _2m75 2G 2 U AA „ 

J fj,5 1^7^75 2 ^ 1 _i_ O/^ ttAA j m^7 7s J 



g 2 1 + 2G 2 U AA "V 2 

/ q,q u 2G 2 U AA x „ r» 

« 2 l + 2G 2 n^H 775 2 



(i + 2G 2 n 



T a 



AA I P" 
2±l T 

= G , A (9 2 )^7 iy 75y. (2-83) 

Thus the axial current is purely transverse and consequently conserved, q^J^ = 0. However 
the quarks have acquired an axial form factor 

associated with the tt + and a\ exchanges that contribute to the current. At zero mo- 
mentum transfer this gives the renormalized axial coupling as 

9a = G A (0) = [1 + 2G 2 n^(0)]- 1 = [1 + SG,/ 2 ]- 1 = [1 - 8G 2 / 2 ], (2.85) 



in view of Eq. (|2.38|) followed by Eq. Q2.79 ). It will prove useful to introduce a renormalized 
form for the axial form factor 

F A (q 2 ) = G A (q 2 )/G A (0) = g A l G A {q 2 ) (2.86) 

that has the property that F A (0) = 1. One notes the essential role of the axial ai degree 
of freedom in obtaining all of these results. In particular, the presence of the diagonalized 
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version of the longitudinal axial meson propagator in Eq. Q2.81Q is crucial for axial current 
conservation. One also observes that the quark's axial coupling to the current has been 
reduced from unity down to g& < 1 by the associated meson clouds. The same remark is 
true of the weak decay constant that is also reduced below its mean field value to f% = gAfp, 
as shown in Eq. (|2.79| ). Thus one has to identify the renormalized f n and not f p with the 
physical value of ~ 93 MeV. 

Vector current conservation follows in the same way, except that here p° exchange replaces 
the 7r and axial modes. Otherwise the calculation is the same as above. However, the "bare" 
current and meson exchange current are separately conserved in this case. Translating Fig. [I] 
once more, one finds 



i + 2G2nV T ^ lU 2 ~ (i + 2G 2 m: y ) T ^ 7l/ 2 ~ if Ws 2 ) : W 2 ' 



^2.87) 



after invoking the identity ^=0 that is appropriate when calculating (p'\J^\p). Thus is 
also conserved, and the vector form factor of the quark is 

^ = 1 + <2 - 88) 
In this case, there is no renormalization at zero momentum transfer, Fy(0) = 1, since 

nlf v (o) = o. 



G. Explicit symmetry-breaking and chiral perturbation theory 

The results obtained so far all refer to the chiral limit of vanishing current quark mass. 
We close this section by briefly indicating their extension to the case m ^ that is required 
for the discussion of some of the sum rules in Section III. As a byproduct, this extension will 
also allow us to make contact and compare with the results of chiral perturbation theory 
(CHPT) that has been extensively used as a theory for studying hadron properties in the 
non-perturbative regime of QCD P, P3"|,|3~8" . 



When a non-zero current quark mass is included, the gap equation for the new quark 
mass that replaces the m of Eq. ( |2.34|) is 

m* -m = E(m*) « E H (m*), (2.89) 

that has a revised solution m* ^ m. The changed gap equation alters the structure of the 
pseudoscalar polarization from that given by Eq. ffTlgD to 
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evaluated at the new mass m*. The new index nc explicitly denotes the use of non-chiral 
functions. Carrying out the same set of matrix multiplications as those leading to Eqs. ( |2.44| )- 
( p.47| ), it is found that the non-chiral expressions for II PP and U AP have exactly the same 
form as their chiral counterparts, but with fl pp replaced by IT PP . Furthermore, the longitu- 
dinal polarization now no longer vanishes since the associated axial current is not conserved. 
The final expressions for the non-chiral case are thus 

Tl pp (a 2 ) 

n pp = nc {q ' (2 91) 

u - i - 2G 1 n pp (^ ) ^ JL > 

frAP _ 2if 2 (q 2 )\f^ 
nc 1 - 2G 1 n pp (g2) ^ • > 

= ; :" JXrl ( 2 - 93 ) 



that do not permit further simplification. All amplitudes now have to be evaluated for a 
quark mass m = m* given by the root of Eq. fl2.89p . 

The propagator matrix representing the mixing of the pion and longitudinal axial vector 
modes has the same form as before, 




<D«=-i£>£V) _ I (2.95) 



with A and B unchanged except for being evaluated at m = m*, but C contains an additional 
term that is proportional to m, 

A = 1 (2.96) 

C = 2G 2 [^g 2 - (2G 1 )- 1 ^[1 - SGWSfo 2 )]]. (2.98) 
The propagator for a non-chiral pion, 

has been extracted. The pion pole has thus been shifted away from the Goldstone mode at 
q 2 = to some finite value at q 2 = m 2 that gives the mass of the ENJL pion. To lowest 
order in the current quark mass, this pole lies at 
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2 mm 



mi^^— 2 +0{m% (2.100) 
The corresponding pole for the pion mass of the minimal NJL model lies at 



mm 



m p = ^aJ~2+°(™)- ( 2 - 101 ) 

The quark mass in these expressions now refers to the mass obtained from the gap equation 
in the chiral limit, m — 0. Using the mean field result, m = — 2Gi(tptp), that connects m 
to the condensate density in either the NJL or the ENJL version of the model in the chiral 
limit, one can recast these equations as 

m 2 J 2 = m 2 p f 2 = || = -m(^) + 0{m 2 ) (2.102) 

to confirm that the Gell-Mann Oakes Renner current algebra relation |39| is satisfied by 
either version of the NJL model §|,|l(| . 

The diagonalization of the propagator matrix proceeds as before, and one regains 
Eqs. ( [2.67D an d ( |2.68| ) with the rotation angle calculated using the revised values of A , B , C 
and with D pp replaced by D pp '. The non-chiral analogs of the pseudoscalar and induced 
pseudovector coupling constants of Eqs. (|2.73|) and (|2.74j) are related to the residue at the 
pion pole of D pp (q 2 ). If we use capital letters to distinguish the non-chiral coupling con- 
stants from the chiral ones, then 



G« qq = lim [dIi P n P /dq 2 }-^\(q 2 ) cosfl (2.103) 

and 

G pv = lim 2m*[dU pp /d q r^ Jxtf) (2.104) 

The modified quark mass m* enters into all these expressions in a very complicated way. 
However, since the current quark mass (m ~ 3 — 5 MeV) is expected to be of order of (2—3)% 
of the quark mass that is generated dynamically via the chiral gap equation, a perturbative 



approach suggests itself [^T|] that is analogous to the procedure used in CHPT. This approach 
seeks to express the non-chiral amplitudes in terms of their chiral counterparts to first non- 
vanishing order in the current quark mass m. Since the relevant physical values of q 2 ~ m 2 
are always of 0(m), which in turn are of the same order as the shift 5m = (m* — m) away 
from the chiral mass, all amplitudes have to be expanded simultaneously in q 2 and 5m about 
zero. Basic to these expansions are the following results: 
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(i) The mass shift. The mass shift 5m to 0(m) is obtained by solving the non-chiral gap 
equation ( |2.89| ) to this order. One finds |2TJ 5m m m/8Gif 2 , a result that can be rewritten 
in either of the two equivalent forms 

5m or 5m = ^ L , (2.105) 

4m 4m 

that lie between 20 and 13 MeV depending on the parameters used, see Table |. The physical 
origin of the mass shift m — > m* that is induced can be understood directly in terms of an 
additional self-energy change due to the a mode being exchanged between the quark and 
an external scalar field So = ipip, that is associated with the current quark mass in the 
Lagrangian. The diagrams entering the gap equation that determine m* to 0(m) are shown 
in Fig. The additional quark self-energy that is due to the non- vanishing current quark 
mass is given diagramatically by Fig. [| (b) and (c): translating these Feynman diagrams, 
one obtains 



i5m = -im + (-im)[-iIl SoSo (0)][-iD SoS ° (0)] 
—im 



1 -2G 1 n 5 o5o( )' 
since the a propagator is given by 



(2.106) 



- iD SoSo (a 2 ) = = — (2 107) 

in terms of the scalar polarization of Eq. (|2.30c| ). The a exchange has generated an effective 
scalar form factor [1 — 2Gin s °' So (0)] _1 that modifies the bare current quark mass vertex in 
an analogous way to what was found for the vector and axial vector current vertices, and 
its structure as a geometric series in G\ confirms that the calculation of 5m is still non- 
perturbative in the strong coupling constant. Inserting the zero momentum transfer value 
of the combination 1 - 2G , iII 5oSo (0) = 8Gi/ p 2 from Eq. (|2T30cl) , one regains Eq. ( [2~T05l) . 



(ii) Non-chiral corrections to the pion-quark coupling constants. The revision of the pion 
polarization to order 0(m 2 ) is obtained by expanding n£f to lowest in q 2 and 5m. To do 
so, first calculate the variation of the basic integral I{q 2 ) in Eq. ( |2.33| ) with q 2 and m 2 . We 
temporarily write I(q 2 ,m* 2 ) in order to emphasize the two independent variables involved 
in this function, and obtain 

^* 2 ) = '<°- m2 > + (4^-(W + ''-' (2 ' 108) 

with the help of the relations 



29 



9/(0, m 2 ) i 1 dl(0,m 2 ) _ i 1 

dq 2 (47r) 2 6m 2 ' dm 2 (47r) 2 m 2 ' 

that are easily established. Here m* 2 = m 2 + 5m 2 + 0(m 2 ), where 5m 2 = 2mSm = \m 2 p1 
from Eq. ( p. 1051 ). Hence / 2 (g 2 ) = —AiN c m* 2 I(q 2 ,m* 2 ) in Eq. (|2.31|) has the expansion 

m m *) = fp[1 + _J_ _ J=L + j£ + OM] , (2 , 09 ) 
from which one concludes that f 2 {q 2 ) of Eq. (|2.41|) behaves like 



fA^m* 2 ) = / 2 [1 + ^q 2 + ^ - J|m 2 + 0(m 2 )]. (2.110) 

We have set N c = 3 in both of these expansions, and used the relations in 2 = gAm^ and 
/ 2 = f 2 /gA to eliminate the NJL values of the pion mass and decay constants in favor of 
those pertaining to the chiral ENJL case. Then one finds that the pseudoscalar polarization 
has the expansion 

rrPPf 2 *2\ 171* — 171 2 f 2 g\ 2 3g A 2 g A {\~QA) 2 . A 2 .. 

2&im* wr on on In 2m 

(2.111) 

Accordingly, the derivative of this expression at the ENJL pion mass m w that enters the 
determination of the pseudoscalar and pseudovector coupling constants is 

2 9 
Ttl Tfl 

[dK c p (q 2 , m* 2 )/dq\^ ml = g£[l - g 2 A (3g A - 2)^ - g A (l - g A )^ + 0(m 2 )]. 

(2.112) 

Now, the final determination of the coupling constants, and thus also the non-chiral correc- 
tions to the pion weak decay constant, also requires the expansion of the root A 7r (m 2 ) as 
well as sm9 and cos# to lowest order in m 2 . To obtain these, we need the expanded forms 
of the coefficients A, B and C of Eqs. Q2.98|) . One then finds that 



Tfl I 

A w (m 2 )^l + (4G 2 / 2 ) 2 ^, ^K(ml) cos 0«1 

(2m*/m 7T )^X 7r (m 2 7r )sme « 8G 2 / 2 (m 2 , m* 2 ) (2.113) 

to lowest order in w, after dropping terms of 0(rh 2 ). Here, the function f 2 (m 2 , m* 2 ) is given 
by the expansion in Eq. ( p. 110] ) evaluated at q 2 = m 2 . Hence the coupling constants take 
the form 

2 2 

in in 
G-nqq = <fcr M [l + gfagA - 2 )^tp + 9a(1 ~ 9a)^ + 0(rh 2 )} (2.114) 
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and 

2 9 

m Tn 
G pv = g„ m {l - g A )[l - Zgl^tp + 9a(1 + 9a)^ + 0(m 2 )]. (2.115) 

As a check, notice that this last relation correctly reduces to ( |2.74| ) in the chiral limit. 

(hi) Non-chiral corrections to the pion mass and weak decay constant. The values of the 
physical pion mass and weak decay constant M n and F n , expressed as expansions in powers 
of m, can now be found. The weak decay constant is obtained by evaluating the right hand 
side of Eq. ( |2.75| ) at q 2 = m 2 instead of zero, while replacing g nqq — g pv by G nqq — G pv and 



m by m*. One obtains 



m 2 f 2 

= + -^(4*n*±& ~ *&) + 0(m 2 )], (2.116) 



to the indicated order. The corrections to the ENJL pion mass of Eq. ( |2.100| ) follow upon ex- 
panding the defining equation for the value of the pole in the non-chiral propagator D^(q 2 ) 
to higher powers than the first in fn: 

m*m 



Mi 



2G 1 /2(M2,m* 2 ) 

m 7r / q 2 fir 

32n 2 f%\ m? 



m l\ l ~ i^f^^—M^A - 1) - 4^(3<m - 1) ) + 0(rh 2 )] (2.117) 



after using Eq. (|2.100 ) to eliminate the combination mm/ (2G\f 2 ) 



2 

ml. 



(iv) The Gasser-Leutwyler coupling constants of chiral perturbation theory. The coefficients 
of the m 2 terms in expansions like Eqs. (|2.116|) and ( p.H7| ) that have been underscored with 
braces can be used to identify the ENJL model predictions for two of the coupling constants, 
l 3 and Z4 of CHPT that were introduced by Gasser and Leutwyler ||. However, as they 
stand, these equations actually determine the renormalized coupling constants Z[, and not 
the scale-invariant versions U of Gasser and Leutwyler, which are the physical parameters 
to be extracted from experiment. The price one pays for introducing the scale invariance 
in li is that the latter coupling constants then all diverge logarithmically like — In m 2 in 
the chiral limit. Such "chiral logarithms", that are expected on general grounds |19| |, only 



appear in the NJL model or its ENJL extension when non-chiral 0(1/N C ) corrections to the 
mean field results are included. As we have not incorporated such corrections into our mean 
field treatment of the problem, we must identify the underscored quantities not with Zj, but 
with the combination 
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^ + ln(m 2 //x 2 ) ~ /[, (2.118) 

which is well-defined in the chiral limit. Here fi is the as yet arbitrary scale parameter 
of CHPTQ This way of including the chiral logarithms has been verified explicitly in the 



calculation of the pion radius parameter Iq in ||20|| , and for the pion electromagnetic mass 



difference in | 15|| . The same conclusion holds for the remaining coupling constants. We thus 
have 

If NJL + Hml/fu 2 ) = {8vr 2 4(2^ - 1) - 4g A (3g A - l)}g A (2.119) 

Tf NJL + Hml/fJ) = {4tt 2 ^ - 3g A }g 2 A . (2.120) 

These expressions reduce to the values given previously (in a somewhat different form) 
for the NJL model (g A — > 1, f 2 — > f 2 and m 2 — > m 2 ): 

l» JL + ln(m 2 //i 2 ) = 8n 2 S- - 8 (2.121) 

m z 

ff JL + Hml/fi 2 ) = 4tt 2 4 - 3. (2.122) 



In Eqs.( p.ll9|) and ( |2.120|) , the combinations on the left are specified as functions of two 



physically measurable quantities, the pion decay constant as well as the axial form factor 
g A . The only unknown factor is the quark mass m. Numerical values for these quantities 
will be discussed in the next section, after the parameters of the model have been fixed. 

H. Parameters. 

The determination of the parameters appearing in the effective Lagrangian follows the 
standard procedure of attempting to reproduce some of the experimental properties of the 
mesons and quarks that can be calculated from the model. In the case of the ENJL model 
with the chiral symmetry-breaking current quark mass included, one has the two interaction 
strengths, G\ and G2, an averaged current quark mass m, and a regulating cutoff A, to 
determine, i.e. four parameters in all. On the other hand, it is clear that the /j's derived 
above to 0(rh 2 ) only refer to the solution for m from the gap equation (|2.34j ) in the absence 



7 The arbitrariness in the choice of the scale fi is reflected in the slow logarithmic dependence on 
it and the overall smallness of the logarithmic correction in any event. This explains the good 
agreement of previous calculations with the empirical values of the l^s, see \21\. 
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of a current quark mass. This is the general situation. We may thus ignore the effect of rh on 
the quark mass, and instead regard rh as interchangeable with the pion mass m v in view of 
the GMOR relation ( |2.1(J2| ) rather than as an additional parameter. This then leaves three 
parameters to fix: G\ , G2 and A. In order to determine these, we use the physical value 
F-x « 93MeV of the pion decay constant^ to fix coupled with the inferred value of the 
condensate density per light quark flavor of (uu) ~ (dd) ~ (— 200MeV) 3 to (— 300MeV) 3 , 
and a typical value for the quark axial form factor to set these parameters. Choosing 
the value of g& is, according to Eq. ( [2.85D , equivalent to choosing a value for the interaction 
strength Gi multiplying the vector fields, once is fixed. We take gA = 0.75, since this 
leads to a correct estimate of the nucleon axial form factor g% = (5/3) (7.4. ~ 1.26 using the 
naive quark model formula ||5T]. Then G 2 is fixed at G2 = (1 — <7a)/(8/^) — 3.61GeV~ 2 . 
The ratio fl/m 2 follows from Eq. (|2.32a|) plus f 2 = gAfp- We use the Pauli-Villars (PV) 
regularization scheme for this calculation. The PV regulated version of Eq. (|2.32a|) is given 
by Eq. (|A34| ) in the appendix. One has 



K 3 9A ^ „ Mi 3g 



m 



2 



a=0 



A 2 

L(x) = 21n(l + x) -ln(l + 2x), x = — , (2.123) 

m 2 

for A^ c = 3, after inserting the standard choice |1| C a = (1, 1, —2) and M 2 = m 2 + a a A 2 with 
a = (0,2,1) for the PV regularization[|. Thus / 2 /m 2 , and therefore also the combination 
Zj+ln(m 2 //i 2 ) are only functions of the one remaining dimensionless combination x = A 2 /m 2 , 
i.e. the ratio of the square of the regularization scale to the spontaneously generated quark 
mass at that scale. Thus one may write 



8 A more accurate value is F w = 92.4MeV, see Ref. Q. 

9 As an aside we comment that, when regulated in 0(4) momentum space, the corresponding 
expression for L(x) is 

Lo(i)( x ) = m (l + x) - x/(l + x) 

Here x = Ap^/m 2 and A ( 4 ) is a new cutoff. The scaling AL 4 j — > (2 In 2) A 2 essentially maps 
this function onto L(x), i.e. Lqu-\[(2 ln2)x] ~ L(x) to a very high degree of approximation, see 
Fig. 13 of Ref. |R|. The results to follow are thus essentially unchanged in the 0(4) regularization 
scheme, apart from requiring a larger cutoff. However, for the reasons already outlined in the 
Introduction, we prefer to use the PV regularization. 



33 



lf NJL ( x ) + \n(m 2 Jfx 2 ) = 6g 2 A (2g A - l)L(x) - Ag 2 A (3g A - 1) (2.124) 
lf NJL (x) + Hml/fx 2 ) = (3L(x) - 3)g% (2.125) 

which reduce to the especially simple expressions 

% JL (x) + ln(mj//i 2 ) = 6L(x) - 8 (2.126) 
If JL (x) + Hm 2 p /ij 2 ) = 3L(x) - 3 (2.127) 

for the NJL values of these coefficients. 

For fixed f n and g A , Eq. ( |2.123| ) determines how m must vary with the ratio x. The 
corresponding value of the strong coupling G\ that is required to reproduce this value of m 
follows from the gap equation that m must satisfy. Using the PV regulated version of the 
gap equation given in the appendix as Eq. (|A4|) for iV c = 3, one finds 



G r i = ^[2xln{i±^}-L(x)]. (2.128) 



Knowing G\ and m, the quark condensate density per light quark flavor is then given by 

— 777 

- (uu) = ~{dd) = - 7r + 0{fn). (2.129) 
4Gi 

The variation with x = A 2 /m 2 of the quark mass m and condensate density is shown in 
Fig. |^ for the ENJL model. The input parameters are as stated above: /„- = 93MeV and 
g A = 0.75. If we insist that f p must take on the value of the physical pion decay constant of 
93MeV for the minimal NJL model case also, then both these quantities simply scale down 
with ^[gX: m NJL = y/gXm ENJL and (qq)]{j L = \fg~A {Q0.)^enjli so tnat A njl = ^/Qa^enjl, 
and {Gi}njl = 9a 1 {Gi}enjl- It is interesting that the p meson degree of freedom takes no 
part in this scaling. One notes that both m and the quark condensate are slow functions of 
x: a fourfold change from x = 4 to x = 16 can accomodate an uncertainty in the empirical 
value of the condensate density that lies between — (300MeV) 3 < (uu) < — (250MeV) 3 . A 
comparison of the parameters and the numerical results for the two cases is given in Table | 
for x = 16. 

In Tables || and |T|, we have compiled the analytic expressions and predicted numerical 
values for the Gasser-Leutwyler scale-independent coupling parameters as extracted from 
Eqs. ( |2.124j ), ( |2.125| ) for l 3 and U for both the NJL and ENJL models, and compared these 



values with experiment. In addition, we have included the values of Z5, Iq and I7 (for which 
there is no barred version) and which are calculated from Eqs. Q3.39 ) and Ql.7|) as well as 
Eq. ( |3.52|) in the sections to come. 
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The predicted values for these parameters have been calculated using the same input 
data as for Table |. The logarithmic contribution is shown separately in this table. This 
contribution has been calculated at the scale \i = 2m a ~ 2m = 485MeV or 528MeV respec- 
tively in these two cases, that is suggested by the structure of the pion radius calculation 
p0| . Gasser and Leutwyler fix p = M v = 549MeV at the rj meson mass. From this table, 
one notes (i) that either version of the model gives very similar predictions for the coupling 
constants of chiral perturbation theory, and (ii) that the agreement with the empirical data 
is quite acceptable, the more so since no attempt was made to fit any of the Vs directly. 
The only exception to this is I3, which is more than twice as large as the suggested empirical 
value which, however, carries a large error. 

III. SPECTRAL DENSITY FUNCTIONS AND SUM RULES 

In this section, we construct the spectral density functions that give the meson mode 
(mass) 2 strength distributions in the various channels, and discuss their sum rule properties. 
Six such spectral functions can be defined and there are eight sum rules linking them that 
we discuss. We first introduce the spectral functions, and we then evaluate them in the 
ensuing subsections within the ENJL model. 

We start with the isovector vector and axial vector polarization given by Eq. ( p.5|) and 
consider their spectral densities p J . We write 

^ J rf 4 xe^(0|\/;(x)V/ 6 (0)|0) = plVX-T,,)^ (3.1a) 
i / d 4 xe^(0\A;(x)Al(0)\0) = p?(q 2 )(-T^)5 ab + p£(q 2 )L,J ab . (3.1b) 

The subscripts and 1 indicate the spin projection carried by these spectral densities. Since 
the normalized transverse and longitudinal projectors and where + = 
have been used, these identifications of the vector and axial vector densities agree with those 
introduced in Weinberg's original paper on sum rules, and equal q 2 times the vector densities 
of Gasser and Leutwyler [Q. 

In a seminal paper ||, Weinberg proved two QCD sum rules that now bear his name, 
and which are obeyed by the inverse moments of these spectral densities in the chiral limit 
in which both the vector and axial vector currents are conserved: 

/ -{p V As)- P ^s)) = fl + 0{m) (3.2a) 
Jo s 

roc 

/ ds{ P X{s)- P t{s)) = 0{m). (3.2b) 
Jo 



35 



The first sum rule is exactly fulfilled in the chiral limit, where, as our notation suggests, f 2 
represents the value of the physical weak decay constant F 2 = f 2 + 0(m) in the limit of 
vanishing pion mass. The second sum rule contains additional QCD corrections |T],^[. These 
relations assume (see Ref. Q for a recent overview) that the associated polarization functions 
satisfy unsubtracted dispersion relations. Then, with the aid of the Kallen-Lehmann spectral 



representation ||Iq| , one can write 

with J = V, A, after setting s = q 2 . The spectral density distributions for the vector and 
axial vector polarizations are then 

p{(s) = -i-/m[n^( S )] , and p£(s) = ^Im[ti A L \s)]. (3.4) 

In addition to the Weinberg sum rules for the vector and axial vector densities, further 
sum rules have been derived by Gasser and Leutwyler for higher inverse moments of these 
densities, as well as for the moments of the scalar and pseudoscalar densities, that are specific 
to expansions in the current quark mass m in the context of chiral perturbation theory ||. 
To introduce these sum rules, we use the same definitions as these authors in the scalar - 
pseudoscalar sector and write 



2vr 
1 

1 

1 

27 



d 4 xe^<0|S (:r)S (0)|0> = p s {q 2 ) (3.5a) 

d 4 x e iqx (0\P a (x)P b (0) |0> = p P (q 2 )5 ab (3.5b) 

d 4 xe iqx (0\S a (x)S b (0)\0) = p S (q 2 )S ab (3.5c) 

d 4 xe iqx (0\P°(x)P°(0)\0) = p p (q 2 ) (3.5d) 



with 



p J (s) = -Im[fl JJ {s)] (3.6) 

7T 

for the scalar and pseudoscalar correlators: here J = S,P identifies the character of the 
relevant density operator. The tilde over the scalar isovector and pseudoscalar isoscalar 
spectral densities in Eqs. (|3.5c| ) and ( |3.5d ) serves to distinguish them from the isoscalar and 
isovector densities of the same spatial character in the first pair of equations. The density 
operators appearing in these correlation functions have been defined in Eqs. ( |2.3| ) and 
Then Gasser and Leutwyler show by standard techniques ]IJ that 
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00 d s I 

~{p V As) - p?(s)) = —(k - 1) + 0(m) (3.7a) 

S~ 407T Z 

roo A a f% 

l ^) = ^| + 0(m), (3.7b) 

to the indicated order in m. Two further sum rules that require a non- vanishing current 
quark mass are 

roo A a 

-m\p s {s) - p p (s)) = -flml + 0(m 2 ) (3.8a) 
o s 

oo Ac 

—rh 2 {p s {s) - p p {s)) = -2m% + 0(m 3 ). (3.8b) 
o s 

As noted previously, f% and m 2 refer to the lowest non-vanishing order expressions for the 
pion decay constant and pion mass, considered as a function of m; l§ and If are simply two 
further parameters introduced by CHPT. These authors also show that 

00 

ds(p s (s) - p p (s)) = 0(m) (3.9a) 

roo 







roo 

/ ds(p s (s) - p p (s)) = 0(m), (3.9b) 
Jo 

that mimic Weinberg's second sum rule for the scalar minus pseudoscalar densities. 



A. Vector densities and their properties 



Recently it was shown []T7 that Weinberg's first sum rule is obeyed exactly, and the 
second sum rule obeyed in modified form, by the vector and axial vector densities generated 
by both the NJL model and its ENJL extension, provided that these densities were calculated 
from polarization functions that were PV-regulated. Here we expand on these, perhaps 
unexpected, results and also extend the discussion to cover the remaining sum rules stated 
above. 

We first examine the vector and axial vector spectral densities for the NJL model. These 
correspond to densities at the one-loop level of approximation. Their difference comes di- 
rectly from the imaginary part of the Ward identity relation in Eq. ( 2.37Q after Pauli-Villars 
regularization. The result is 



(o)v, x (0)A/ n N c m 2 Am? 2 N c m 2 A 4 

Pi (s)-pi (s) = -^r {y ~ IT ( s ^ PV ^72"' s ^°°' ( 3 - 10 ) 



which is derived in detail in the appendix as Eq. ( A55 ), and where the PV notation is defined. 



Without the PV-regulating instruction in place, this difference reaches the constant value 
N c m 2 /(4:7c 2 ) at large s, and neither sum rule would converge. With it, this difference is 
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sufficiently convergent to give finite answers for both sum rules. By direct integration, one 
in fact finds that 

r-(pr(s)- P r(s))=f P , (n JL) , (3 .n) 

Jo s ' 

while 

r ds(pf )v (s) - p ( i )A (s)) = -m(^> + . (NJL). (3.12) 
Jo 

Details of this calculation can be found in Eqs. ( |A57| ) and ( |A59| ). Thus the PV-regulated 
vector spectral densities of the minimal NJL model continue to obey the sum rules even 
though these one-loop densities do not contain any vector or axial vector resonances struc- 
tures of any kind. One sees that the first Weinberg sum rule is obeyed exactly, with f% 
replaced by its mean field equivalent f£ appropriate to the NJL model. In particular, one 
realises that the PV regularization of the NJL spectral densities has introduced precisely 
the correct amount of additional negative or "unphysical" spectral density to satisfy the 
first Weinberg sum rule exactly^. The second sum rule is obeyed in a modified form that is 
connected with the existence of the quark condensate in the NJL ground state. Only for the 
case that the condensate density vanishes, does the second Weinberg relation hold exactly. 
The value of the right hand side of Eq. (|3.12|) is typically (250MeV) 4 in the broken phase. 

The situation for the ENJL generated spectral densities is quite different. Now quark- 
quark interactions are taken into account that effectively screen out the undesirable high 
energy behavior of the one-loop densities, as is shown below. This in turn means that two 
possibilities are available: (i) to PV-regulate both the real and imaginary parts of the irre- 
ducible polarization functions, which implies that these satisfy an unsubtracted dispersion 
relation involving PV-regulated spectral densities. Polarizations and associated quantities 
that are calculated in this scheme are denoted as "barred" functions in the following, (ii) 
Alternately, one can regulate only the real part of the irreducible polarization functions, 
leaving the imaginary part and thus the spectral densities unregulated. Such a polarization 



10 We also remark that the same sum rules for the NJL densities are obtained if the the polarization 
functions Tlj, v (s) and Hj< A (s) that generate them are regulated via an 0(4) cutoff in momentum 
space. However the problem then is that the 'Ward identity", Eq. ( 2.37| ), that lies at the root of 
both sum rules does not necessarily hold without introducing additional assumptions to deal with 
the non- uniqueness of regulating quadratic divergences under 0(4). This difficulty is discussed 



further in Refs. |37]] and 15]. We do not employ the 0(4) regularization method in this paper. 
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function satisfies a once subtracted dispersion relation containing unregulated spectral den- 
sities. These quantities are denoted as "unbarred" in what follows. Details of the schemes 
are given in the appendix. There it is shown explicitly that, from the point of view of 
dispersion relations, the distinction between the barred and unbarred amplitudes can be 
traced back to using either the PV-regulated or unregulated one-loop spectral densities in 
the dispersion relation calculations given in Eqs. ( |A48| ) and (|A52|) . 

We now examine the ENJL-generated spectral densities as regulated in the barred scheme 
and show in more detail that the same form of results as was found for the NJL model in 
Eqs. (|3. 11| ) and ( |3.12|) also hold in the ENJL case. The vector and axial vector densities are 
obtained from the corresponding polarization functions that include all interactions in the 
RPA approximation, as recorded in Eqs. ( [2.1 1| ) and ( |2.13| ). Then 

™ = -><ttwm) < 313 > 



The barred polarization functions are given in Eqs. (|A46b| ) and ( |A54a ) of the appendix. 
Carrying out the indicated operations, one finds that 

pY(s)=p? )V (s)\F v (s)\ 2 (3.15) 
pf(s) = p? )A (s)\G A (s)\ 2 , (3.16) 



where Fy and Ga are defined in Eqs. ( |2.88| ) and ( |2.84j ) with IT replaced by its regulated 



version ft. Thus the densities of the ENJL model are simply those of the minimal NJL 
model, modified by the modulus squared of the corresponding quark current form factor. 

Plots of the ENJL vector and axial vector densities are shown in Fig. |] using the typical 
set of parameters given in the second line of Table | that are appropriate in the chiral 
limit (we refer in the following to this set of parameters as the standard set). The different 
threshold momentum dependence of the spectral densities in Fig. |4] comes about via the 
behavior of the factors p± V or p^ A . These factors behave like p 2L+1 = p or p 3 , where 
p = (l/2)Vs — 4m 2 is the decay momentum for the decay of the spin one bosons with 
J PC = 1 or 1 ++ respectively into a qq pair. The L dependence comes about since the 
parity and charge conjugation quantum numbers of the qq system are given by P = — (— 1) L 
and C = (— 1) L+S . The opposite parity and charge conjugation in the two channels then 
requires that the spins of the qq pair be aligned to S = 1 in both cases, but also to carry 
internal angular momentum L = 1 in the second case. One sees that the presence of the 
quark form factors drastically redistributes both spectral densities, causing a pile-up of 
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strength at low momentum transfers s = q 2 by suppressing the high energy tails of p^ V 
and p^ A . However, at large values of s, both spectral densities become negative due to 
the asymptotic behavior of the regulated one loop densities p^ V and p{ . This is shown 
explicitly in Fig. ||] for the vector density of both the ENJL case and the (one loop) NJL case. 
The strong redistribution of the vector strength due to the inclusion of the vector mode is 
especially clear in this figure. In this case, the redistribution of strength takes place in such 
a manner that the integral over p\(s)/s always vanishes. 

In Fig. we have shown the difference p^ — p[ A \ plotted out to values of s that extend 
beyond both PV-induced thresholds, to illustrate the extent of modification introduced by 
the vector and axial vector degrees of freedom. Again one sees that there is a pile-up of 
strength at low momentum transfers. 

If one were to follow the second, unbarred regularization scheme, then the PV-regulated 
one-loop NJL spectral densities of Eqs. ( |A50|) and (|A56|) are replaced by their unregulated 
versions 



miv, \ N c s r / 4m 2 , 1 4m 2 9ni , 

Pl (S) 24^ { V 1 " ^ (1 + 2— )6{s ~ Am )} (3 - 17) 



and 



in Eqs. ( |3.15| ) and ( ^.16[ ). This effectively removes the bars on Fy{s) and Ga(s). Since 

Pi(s) ~ P A ( S ) ~ n \ h2 » |s|^oo, (3.19) 
s(ln \s\) z 

from Table 0, one then obtains spectral densities for the ENJL model that, on the one 
hand, always remain positive definite, instead of behaving like — 1/s, but which on the other 
hand, are essentially indistinguishable from those given in Fig. [| over the resonance regions 
on the scale of this figure. Later we discuss some consequences of making this choice of 
regularization instead. 

If one associates the peaks in either the regulated or unregulated spectral densities with 
the masses squared m 2 and m 2 ai of the p and a\ modes of the ENJL model, then for our 
standard set of parameters, one finds m p = 0.713GeV and m ai = 1.027GeV from Fig. |j. 



This is in fact in reasonable qualitative agreement with the measured values |33[ of these 
meson masses, which lie at 0.768GeV and 1.230GeV respectively. However, since both 
spectral densities peak above the qq decay threshold at 4m 2 , which is unphysical from the 
point of view of confinement, the identification of these peaks with the meson mode masses 
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can be questioned. On the other hand, we show below that Weinberg's first sum rule is 
satisfied exactly by the associated vector minus axial vector densities that exhibit these 
peak values, so that the integrated strength is exactly reproduced in accordance with the 
tenets of chiral symmetry, on which this sum rule is based. In fact, in Weinberg's original 
treatment [[|, the meson masses were determined by saturating this sum rule with delta 
function strength distributions for the vector and axial vector spectral densities that simply 
ignored the presence of any decay channels as a first approximation. Since the ENJL spectral 
densities also saturate the first sum rule, it is thus reasonable to use the peak values of these 
spectral densities (that have been broadened by replacing the physical p — > nn and a\ — > pn 
decay channels by the unphysical qq decay channel in each case) to identify the associated 
meson masses in the ENJL model, with, as we have seen, not dissimilar predictions for these 
masses. 

The masses of these modes should on the general basis of propagator theory also be 
related to the poles of the corresponding propagators. Below we investigate the analytic 
structure of the form factors Fy{s) and Fa(s), or equivalently the transverse vector and 
axial vector propagators D)f v (s) and (s) given in Eqs. Q2.52| ) and ( |2.54| ). 



1. Pole structure of the p and a\ propagators 



One knows on general grounds |TB| that any complex poles in the scattering amplitudes 
or propagators Dj, v (s) and D^ A (s) that correspond to physical resonances must be located 
on their second (or "unphysical" ) Riemann sheets of the cut s-plane, for reasons of causality. 
From Eqs. ( |2.52j ) and ( |2.54j ), one sees that the poles of these functions as a function of s are 
governed by the roots of the combinations 

l + 2G 2 n£ y (s) and l + 2G 2 L T ^ 4 (s). (3.20) 

We will find that for a range of coupling strengths G± and G2, including those in our standard 
parameter set of Table |, the axial vector mode has a single pole in the lower half plane of the 
second Riemann sheet, while the vector mode has two poles on that sheet: a real root below 
4m 2 , corresponding to a virtual bound state, and a second complex root in the lower half 



plane. The positions of the two p poles and the a\ pole are given numerically in Table [TV 
for the standard ENJL parameter set. Details of this analysis now follow. 

We discuss the vector mode first. According to Eqs. ( |A46b[ ) and ( |A46c| ), the polarization 
function for this channel is given by 

Wis) = -\S- 2 s + + (3 - /)( V / 7coth- 1 ff - 1)} PV , (3.21) 

3 m z bn z 3 v v 
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where / = 1 — 4m 2 /s. This form of the regulated vector polarization function satisfies an un- 
subtracted dispersion relation exactly of the form given in Eq. ( |3.3| ), if the PV-regularization 
instruction is left in place on the second factor. If this instruction is removed (we indicate 
this by removing the bar) then Wf v (s)/s satisfies a once-subtracted dispersion relation. A 
detailed derivation of these statements is found in the appendix. 

The analytic properties of n^ y (s) are determined by the behavior of the function 

J(z) = v /7coth- 1 v /7, / = 1 -p 

that is defined in Eq. ([All]) , with z = s/4m 2 . Considered as a function of a complex variable 
z, this function is single-valued in the cut 2-plane with the cut extending from z = 1 to oo 
along the real axis. The schedule of forms for J(z) and its analytic continuation J(z) onto 
the contiguous Riemann sheet connected through this cut is given in Eqs. (|A20|) through 
( |A21|) . These expression show that both J(x) and J(x) are real in the interval < x < 1 
on the real axis. Plots of the real part of J(z) and J(z) as a function of the real variable 
z = s/4m 2 , are shown in Fig. [7|. 

The corresponding behavior of the real part of Wf v (s) is shown in Fig. |8| as a function of 
the real variable z = s/4m 2 on the first and second sheets of the cut s-plane. For s > 4m 2 
this function has to be continued onto the the upper lip of the cut in accord with the 



Feynman prescription |Tj| . The values of the polarization along the upper lip of the cut on 
the first sheet join smoothly with their analytic continuation onto the lower lip of the cut 
on the second sheet. At the branch point, these values bifurcate into two branches of the 
function that assume different values along the sector < z < 1 of the two sheets. The 
values on the first and second sheets are reached by allowing the variable z to pass either 
infinitesimally above or infinitesimally below the branch point through the cut to reach this 
sector of the real axis. As the vector polarization in Eq. (|3.21 ) also has contributions from 



J's at different arguments s/4M 2 , a = 1,2 due to the PV-instruction, further cuts starting 
at z = 1 + A 2 /m 2 and z = 1 + 2A 2 /m 2 have to be introduced. However, since A 2 /m 2 
is typically ~ 16, both of these cuts commence too far away to the right of the origin to 
essentially alter the behavior of lVf v (s) as determined by J(s/4m 2 ) alone in the interval 
< s/4m 2 < 1. Hence the cusp-like structure where the values of ReTl^ v (s)/4m 2 and its 
analytic continuation meet at s/4m 2 = 1 in Fig. || is a straightforward reflection of the 
behavior of Re J(z) in Fig. [f[ 
The roots in z = s/4m 2 where 

1 + 2G 2 Tl vv (s p ) = 0, (3.22) 
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determine the poles of the p propagator. We now verify that the vector mode has two poles 
on the second sheet, one on the real axis below 4m 2 , and a second complex pole in the lower 
half plane. As pointed out in J3U||, the real pole corresponds having the quark- ant iquark 



pairs making up the vector meson bound into a virtual, (or "antibound" ) state [^6] . 

We briefly discuss this state first: according to Fig. the function Wf, v (z)/4m 2 is real 



and negative on both of its branches for real z in the interval < z < 1 . Hence Eq. ( |3.22| ) 
always has one real root lying in this interval on either the first or the second sheet depending 
on where the ratio Wf v (z)/4m 2 crosses — I/8G2ITL 2 . The solid circles mark the intercepts 
of Wf v (z) / 4m 2 with the value of this quantity as the interaction strength G2 increases 
through the indicated set of values. The value of m has been kept fixed at 264 MeV for 
this illustration. These intercepts determine the position of the real pole of the vector mode 
propagator lying in the interval < s/4m 2 < 1 as a function of the interaction strength. 
One sees from the figure that as G2 increases from zero there is only a second sheet real 
pole, corresponding to a virtual bound state mentioned above for the qq pairs making up 
the p. The binding energy of the virtual state decreases to zero as G2 increases through 
the critical value = — l/2n^ y (4m 2 ) = 7.41 GeV -2 in the present case. After this, the 
pole reappears on the first or physical sheet where it represents a true bound state of 
the qq system whose binding energy again increases towards the limiting value of 2m with 
increasing G2, which it reaches^ as G2 — > 00. We have plotted the motion of the pole versus 
the interaction strength G2 in Fig. |^ with G2 as ordinate, in order to best illustrate the 
turn around of the binding energy as G2 passes through its critical value. For our standard 
parameter set given in Table |, this root corresponds to a virtual bound state at the pole 
position s p = (0.166) 4m 2 = 0.046 GeV 2 , and therefore of mass ^/s^ = 0.215 GeV. 

By contrast, the axial vector propagator has a single complex pole, which we discuss 
together with the complex pole of the vector mode. The axial vector polarization can be 
obtained directly from the Ward identity relation (|2.37|) and the regulated form of the pion 



form factor Fp(s), that is given in Eq. ( A35b ). One has 



U^(s)=U^(s)+4f 2 F P (s) (3.23) 

where 



11 As a consequence, the ENJL groundstate always stays stable against the development of an RPA- 
like instability where becomes pure imaginary, which would be the signal for the formation of 
a "vector" condensate. 
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N r m 2 , ,. , _j 



= [1 - i ^{2 v //coth- i ^ - 2} py ]. (3.24) 



Contrary to Eq. (|3.22|) , the condition 

1 + 2G 2 n£ A (s ai ) = 0, (3.25) 

that determines the poles s ai of the ai propagator, does not have any real roots on either 
sheet, since ReF(s) is positive on both sheets, and overrides the negative vector polarization 
along the real axis. Then fly (s) is real and positive for < s < 4m 2 and no solutions of 
Eq. ( p.25| ) for real z are possible. However, both Eqs. ( |3.22| ) and ( |3.25| ) have complex 



roots located on the second Riemann sheet which may be found by replacing both of the 
polarization functions in these equations by their analytic continuations through the cut. 
In Fig. the variation of the imaginary versus the real part of the second sheet poles of 
the p and a\ propagators obtained in this way is shown with increasing interaction strength. 
One observes (i) that the p meson pole has a real part that always lies below 4m 2 , for the 
indicated range of G2, and (ii) that the imaginary parts of both complex poles are of the 
order larger or equal to their real parts. 

Since the poles themselves are of order unity in units of 4m 2 , and this is much smaller 
than the cutoff scale 4A 2 , we also want to make the point that the numerical results for the 
p and a\ meson poles, as well as the values of the polarization functions in their vicinity, 
are actually insensitive to the replacement of these functions in Eqs. ( |3.22| ) and ( |3.25 ) by 



their "unbarred" counterparts n^ y (s) and U AA (s). In particular, they are insensitive to 
the actual value of the cutoff employed. This is illustrated in Table 0, where the two 
p poles and the single a\ pole are given in units of 4m 2 for our standard parameter set. 
The differences are insignificant. One can see the reason for this more clearly by making 
comparative plots of the vector and axial vector polarization functions generated via their 
unsubtracted and subtracted dispersion relations. From Fig. |ll], one observes that these 
two versions of both the vector and axial vector polarizations are essentially identical in the 
interval — 1 < q 2 /Am 2 < 1, where the influence of the additional branch points introduced 
by using PV-regulated densities have not yet made themselves feltQ. 



12 It should be pointed out that in a low energy expansion ||.13,42] where the non-analytic terms 
are dropped from TOf v (s) and n^(s), the ENJL model leads to approximate meson poles that are 
real. One then finds delta function spectral densities of strengths gy and g 2 Ai that automatically 
obey the Weinberg relations gy = g 2 Ai and g v m v 2 — g\ m A 2 = f 2 (as distinct from the sum rules). 
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One can state this result in yet another way: removing the PV-instruction from the 
second term of the vector polarization function Wf v (s) in Eq. ( |3.21| ) is tantamount to 
removing the explicit cutoff dependence in the vector polarization function, while leaving it 
buried implicitly in the quark mass m in the first term. The values of the roots themselves 



To be precise, 

4 1 6gAm' 2 ' 4 1 6gAm 



from Eqs. (3.21) and ( |3.23 ). This approximation replaces the exact curves for U)f v (s) and Hj, A (s) 



shown in Fig. 11 by tangent lines at s = 0. This completely suppresses all meson dynamics associ- 
ated with the branch cut. In particular, all information regarding the properties of the second sheet 
poles of the p and a\ meson modes is lost. Clearly such an approximation is qualitative at best, since 
it simply ignores the unphysical decay channels that reflect the non-confining nature of the ENJL 
Lagrangian. Its use leads to meson-quark coupling constants g~ qq = g^qq = f 2 /Q9A~m 2 as given by 
the common derivative (—l/4)dH/ds for V and A, and meson masses rriy = 6gAm 2 /(l — gX) and 
m\ = 6m 2 / (1 — gA)i where m\/m v = gT 1 . The factor 1/4 arises in calcuating the coupling con- 
stants if we take the standard form 43] of the interaction Lagrangian as Ci n t = 9pqq('ip'y^T b ip)p^ or 
9aiqq{^p7^l5^T b ' l ( ; )0'i-b f° r coupling the quark currents to the meson fields. Because of the approxi- 
mations, these meson modes have delta function distributions as their spectral densities, of equal 
strengths gy = g\ = (m v / ' g P qq) 2 = f 2 m v /(6gAm 2 )- Inserting this information into Weinberg's 
sum rules, one finds that his two relations quoted above are satisfied (this is not yet a check of the 
sum rules themselves for the ENJL model, as the approximations involved are inappropriate). The 
universality relation g pqq = f p between the p meson coupling constant and its decay constant f p , 
reflecting vector meson dominance of the electromagnetic current, is also recovered from the 
vacuum to one meson state matrix element |jj (n p is the vector meson polarization) 

that governs the p° — > e + e~ decay. After inserting the values of g pqq and fl^fm^), one finds 
9pqq = fp = f^V^BAm 2 that also can be written in the form f p f n = my(l — gA) x ^ 2 - For 
gA = 1/2, this is exactly the KSRF relation 44]. Then one also recovers Weinberg's original 
estimate myi/my = v2 for the mass ratio, as well as the p meson mass my = y6m of the bosonized 
N JL model Ef| . For the standard parameter set of Table | this low energy expansion limit gives 
m v = 1.12 GeV (0.768), m A = 1.29 GeV (1.230), and f p = 6.02(5.10) Q. The experimental values 
are in brackets. 
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are then solely determined by the physical input parameters f n and qa once m (which 
of course contains the cutoff implicitly) has been calculated from the gap equation. The 
differences between the numbers appearing in the first and second rows of Table [TV] thus 
reflect the extent of the explicit cutoff dependence of the roots. However, we learn another, 
equally important feature, from these plots: both the vector and axial vector polarization 
functions that arise from subtracted dispersion relations become negative for large space- 
like q 2 . For the standard parameter set, this has occurred at —q 2 /4m 2 ~ 10, or y/ — s = 
1.7 GeV, but it is clearly a general property. This can be traced to the different asymptotic 
behavior displayed by the barred and unbarred polarization functions, U.(s) ~ — (1/s) ln(— s) 



versus H(s) ~ sln(— s), see Table [VI]. The former expression stays positive for space-like s, 
vanishing at infinity, while the latter function becomes large and negative. This feature in 
both polarization functions generated via once-subtracted dispersion relations means that 
the unbarred versions 1 + 2G 2 U^ v (s) and 1 + 2G 2 ~tt A , A (s) of Eqs. ( ggH) and (ggp each 
develop an additional real root on the physical sheet at space-like q 2 , corresponding to a 



Landau ghost [18| and as shown in Fig. O, i.e. 



1 + 2G 2 U^ V (-M 2 g ) = , and 1 + 2G 2 Il AA {-M 2 ig ) = 0. (3.26) 

We remark in passing that the same is true of the unbarred combination 1 — 2GiIl pp (s) 
that determines the pion pole. Thus the associated p and a\ propagators, as well as the 
pion propagator, contain Landau ghosts if the polarization functions are left unbarred. For 
the standard parameter set, the vector ghost poles occur at the nearly degenerate positions 
s pg = -18.30(4m 2 ) = -5.11 GeV 2 and s ai9 = -18.04(4m 2 ) = -5.05 GeV 2 . The pion ghost 
occurs at s wg = — 13.48(4m 2 ) = — 3.77GeV 2 . The magnitudes of the ghost poles are always 
of C(4A 2 ). 

2. Pole approximations for the vector and axial form factors 

In order to understand the connection between the poles of the form factors and the peaks 
in the spectral densities, we now approximate the former by single pole approximations and 
compare the results with the exact form factors. Their single pole approximations are given 
by 



with residues R p = [2G 2 T\- { ' )VV [sp]}' 1 and R ai = {2G 2 Tl^ AA (s ai )]- 1 . The prime indicates 
the derivative with respect to s of the analytically continued polarization function. Here 
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s p can be either the virtual bound state pole, or the complex pole of the p propagator. 
In principle, both poles influence the shape of the spectral function, and Fy(s) should be 
displayed as the sum of the contributions coming from these two p poles: the mass of the p is 
not uniquely determined by either one. However, the virtual pole, which for our parameter 
set lies at at s p = 0.166(4m 2 ), is not close enough to the branch point at 4m 2 to exert a 
strong influence on Fy(s), so that one expects the complex p pole to dominate in the present 
situation. One should emphasize, however, that this circumstance is parameter-dependent: 
a virtual bound state just below 4m 2 could have a dominant influence on the vector strength 
function. The residues themselves are complex numbers with moduli \R P \ = 0.75 GeV 2 and 
\R ai \ = 0.64 GeV 2 for the standard parameter set. 

A comparison between the modulus squared of the exact form factors and their approx- 



imations using only the complex p and a\ poles is presented in Figs. |13] and [14]. The form 
factors |Fy(s)| 2 and |GU(s)| 2 are faithfully reproduced by their single complex pole approxi- 
mations for our parameter set. The resultant discrepancies between the peak positions in the 
spectral densities of Fig. |] and the pole positions in Fy(s) and Ga(s) are shown in Table [V]. 
These come, not surprizingly, from the distortion introduced by the kinematic threshold 
behavior of the one- loop densities pf^ V and p\ . These factors act multiplicatively to force 
the spectral densities to vanish at the qq threshold, and thus move what would have been 
a pole peak to somewhat higher energies. This is particularly the case for the p resonance 
pole since, as seen in Fig. |T^, its real part is "invisible" below the branch point. In fact, if 
one ignores the threshold behavior of the one-loop densities as a first approximation, they 
become proportional to s and would produce peaks in the spectral densities at the moduli 
of the pole positions \s p \ and |s ttl |. Using the energies given in Table |V|, one obtains the 



estimates y\s p \ = 0.63GeV and y |s 01 | = 0.93GeV, to be compared with the spectral density 
peak positions at 0.71GeV and 1.03GeV. 



B. Weinberg sum rules 

We now investigate the validity of Weinberg's two sum rules for the ENJL generated spec 
tral densities. Taking the difference p\(s) — Pi(s) directly from Eq. ( |3 . 1 5[ ) minus Eq. (R7TBJ) 
one finds that 



Pi(s)-pf(s) 



1 

Im 

47T 



l + 2G 2 nr(s))(l + 2G 2 II™(s)) 



(3.27) 



or 



Pl(s) -Pf(s) = -Im[flF P {s)F v {s)F A {s)}. 

7T 



(3.28) 
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in view of the Ward identity Eq. ( |2.37|) plus the expressions in Eqs. (|2.88|) and ( |2.86|) for the 



vector and axial vector quark form factors Fy and Fa- Notice that the renormalized pion 
decay constant f% = gAfp of Eq. ( f2.7!J| ) appears quite naturally in this result. The inverse 
moment integrals of this difference that are required for investigating Weinberg's sum rules 
can be evaluated in closed form using the methods of contour integration. Their evaluation 
depends on the fact that the barred amplitudes in 

f(s) = f*F P (s)F v (s)F A (s) (3.29) 

that have been constructed using barred polarization functions, are all analytic functions of 
s in the cut plane, with the cut extending along the real axis from 4m 2 to oo. This is so 
because the poles of the vector and axial form factors only occur on the second sheet, and 
Fp(s) has no poles. Hence the integrals of f{s)/s and f(s) taken around any closed contour 
C on the first or physical sheet that excludes the cut (and the pole at s = in the first case) 
are zero. We break up the closed contour into C = Cq + as shown in Fig. [TJ|, where 
Co is indented to exclude the cut plus any poles as necessary, and is closed by the contour 
Coo at infinity. Then the vanishing of the contour integrals of f(s)/s and f(s) leads to the 
relations 



/ 2 -Im[f(s + ie)} = /(0) - — / -f(s) (3.30) 

7T J Am, 2 S 2m JCoo S 



and 



lr 00 If 

- / dslm[f(s + ie)] = -— dsf(s) (3.31) 

7T Jim 2 ' 2m JCoo 

for the integrals involving Im[f(s+ie)} taken along the upper lip of the cut. FromEq. ( gjgp , 



it is clear that the left hand sides of these two relations are exactly the integrals that we 
need for the two Weinberg sum rules. In the first case, /(0) = / 2 Fp(0)-FV(0)Fa(0) = fl 
from the definition of f(s) and the normalization of the three form factors. Moreover, 

/00--^, |s|-oo, (3.32) 

according to Eq. ( |A36| ) , coupled with the information from Table that Fv(s) , GU(s) — > 1 
asymptotically. Hence the contour integral along the great circle will not contribute in the 
first case, but does so for the second, since the asymptotic form of f(s) shows that there is 
a pole at infinity on of residue mtyif)). Consequently 

/•OO A a 

\ -(/%(*)- />}{*)) = ft, (ENJL) (3.33) 
Jo s 
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da(pX(s) - p? (s)) = -m(V^) ^ , (ENJL). (3.34) 

Thus the first sum rule is satisfied exactly by the spectral densities that have been obtained 
by PV-regulating the one-loop densities with the concomitant introduction of negative spec- 
tral density at high s; the second sum rule again assumes a modified form. 

As has already been remarked, the ENJL spectral densities become well-behaved asymp- 
totically if one replaces the one-loop densities in Eqs. ( |3.15| ) and ( |3.16| ) by their unregulated 
versions. This is equivalent to removing the bars on the form factors in Eq. ( |3.29| ), which 
changes the asymptotic behavior of f(s) from what it was in Eq. (|3.32|) to 

f(s) -> — - - / , \s\ -> oo, (3.35) 
s z ux[— s) 



according to Table [V]]. This is sufficiently convergent to prevent the contour integrals along 
Coo from giving any contribution to either sum rule. However this comes at a price: now 
both the vector and axial vector form factors have developed ghost poles at space-like s 
that bring in contributions of their own, and both of the above sum rules will be modified^ 



13 This immediately begs the question of by how much. The answer to this question depends on 
what one chooses to understand by "spectral density" . If we simply mean the unregulated versions 
of Eqs. ([TL5|) and d3T6| ), then 

/ , -(pY(s) ~ Pi(s)) unreg = fl + R P9 + Ra l9 = 0.96/2, 

Jim 2 S y 

either via contour integration, or simply by direct numerical integration, where R pg and R ai g are the 
residues of the function f(s)/s of Eq.( 3.29j ) at the two ghost poles at s = —M pg and —M 2 ig . Thus 



the unregulated densities, which as we have remarked always remain positive definite, saturate 
96% of Weinberg's first sum rule for our parameter set. This equation is, however, open to an 
alternative interpretation. Writing the unbarred version of f(s) in Eq. ( |3.29j ) in the equivalent 
form 

r 1 i _ C pg —C ai g 



8G 2 L l + 2G 2 n^(s) l + 2G 2 nf 4 (s) J s + Mj g s + M% ig ' 

where C pg = {16G| [dU^ v /ds]_ M 2 j" 1 = 0.0648 GeV 4 and C aig = {16G| [dU AA /ds]_ M 2 j" 1 = 

P9 a l 9 

0.0623 GeV 4 are the spectral strengths of the ghost modes, the residues R pg and R ai g of f(s)/s 
are —C pg M~ g and C aig M~ ig . The spectral densities themselves are 5-function distributions at 
space-like s of negative strength, 

pY(s) g host = -Cpg 5(s + M 2 pg ) and p A {s) ghost = -C aig S(s + M% ig ) 
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by these ghost modes. Thus the unregulated ENJL densities that are otherwise physically 
well-behaved have introduced different unphysical elements into the sum rules again by a 
different route. 



C. Gasser-Leutwyler sum rules 

The remaining sum rules can be evaluated by the same method. We start with the 
Gasser-Leutwyler sum rule in Eq. fl3.7a| ). This has a double pole at the origin, so that the 



derivative /'(0) of f(s) at the origin appears on the right hand side of Eq. (|3.30| ) instead of 
f(0). There is also no contribution from C^. Thus 

/ 2 -(P\(s) ~ Pi(s)) = f"[F' P (0) + *y(0) + F' A (0)l (3.36) 

Jim 2 S z 

where the primes indicate the derivatives of the form factors. Straightforward calculations 
yield 



if we assume the standard (— Mg) — it prescription for ghost poles too. Then, by transferring the 
residue contributions from the right hand to the left hand side of the sum rule and including them 
under the integral, the sum rule may be rewritten in the equivalent form 

00 ds , v 



~ (Pi ( s ) Pi ( s )) ghosts+unreg ~ fir 

so that the sum rule violation is uniquely attributable to the presence of the ghosts. Exactly the 
same arguments yield 

roc 

J^ 2 ds( P Y(s) - pt(s)) unreg = C pg - C aig = (224MeV) 4 



or 



ds( Pl (s) - p x {s)) ghosts+unreg = 

for the second sum rule, since C pg — C aig is just the negative of the sum of the residues at the two 
ghost poles, f(Mpg) = C pg and f(—M% ig ) = —C aig , where f(s) is the product in Eq. ( p.29| ), and 
there is still no contribution from the contour at infinity. Thus if one widens the definition of 
spectral density to admit contributions from the ghost modes, both forms of Weinberg's sum rules 
are recovered without modification. These results suggest that the unphysical ghosts represent a 
legitimate part of the meson mode spectrum of the model propagators when subtracted dispersion 
relations are used to construct them. 
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^(0) = ^[l + O(m 2 /A 2 )] 

^(o) 1 ' l ~ g ^ 



6m 2 g A 

Fa(0) = -^(1 - 9a) ~ ~ 9a)[1 + 0(m 2 /A 2 )], (3.37) 



leading to0 



°° 2 ^(pI(s) ~ P f(s)) = f| + £~/—^), (3-38) 
4m 2 s z qti z bm z gA 



after setting SG 2 f 2 — 1 — g A again, and dropping the inconsequential 0(m 2 /A 2 ) correction 
to Fp(0) that appears through the explicit cutoff dependence of the pion radius parameter 
in Eq. ( |A40| ) . On comparing with Eq. ( |3.7a| ), one identifies the ENJL value for l 5 from the 
expression 



lf NJL + Hml/v 2 ) = l + Qg 2 A + 8tt 2 4( L ^), (3-39) 

m 9a 



after adding in the chiral logarithm as in Eq.( p.ll8|) again. Then one finds 1§ NJL = 12.8 for 



the standard parameter set. The empirical value is l§ = 13.9 ±1.3 ||. 

Similar arguments suffice for checking the next sum rule, Eq. (|3.7b|) , except that here we 



have to include non-chiral effects to order 0{m) in the calculation of p$ = (l/47r) Im ff^ 4 



from II^.„ C of Eq. (|2.94j) . Note that has acquired a pole at the physical pion mass due 



to explicit chiral symmetry-breaking by the current quark mass. The behavior of (1/4)I1^ C 

,2 i /n^2\ 



tAA 

V L;nc 

and the contribution to the spectral density in the vicinity of the pion pole at s = m 2 +0(rh 2 ) 
are 

\^fnc « ~ Ca 2 \ • i P$(s), P oie = C aL 5(s - m 2 ) , C aL = /> 2 , (3.40) 
4 s — +ie 

while the behavior on is 

iftfc «^ = -=*£, H^oo. (3.41) 

4 ' m F s s 

In obtaining these results, we have used the Goldberger-Treiman and GMOR relations, 
Eqs. fl2.80|) and (|2.102|) , that are satisfied by the ENJL coupling constants and pion mass, 



14 Since the ghost poles are so massive, using unregulated densities in this case would add a 
completely negligible contribution of C P gM~^ — C a \ g M~^ g ~ 10 -5 to the right hand side. 
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plus the information provided in Table |V|. The pion ghost pole depicted in Fig. [12] does not 
occur, since we are using barred polarization amplitudes. 

We also need the residue of (l/4s 2 )Il^ c at the double pole at s 2 = 0. A short calculation 
gives this residue as 

\ dg4 L;ncts=0 8 7T 2 m 2 / 2 (m 2 ) 8 7T 2 m 2 1 8tT 2 / 2 ' m 2 ' K ' ' 

after using the defining equation ( |2 . 1 1 7] ) for the pion mass up to order 0(m 2 ). Deforming 
the contour in Fig. [TJ| so that the pion pole on the real axis at s = m 2 < 4m 2 is also 
excluded, one finds 



oo w , ; / * + 2 



/7s r°° r/s fl f 2 —C! f 2 

%pt{*) = / ~M(sU ole + pt(s) cut ) = ^ + {4 + -^} = 4 (3.43) 
os 2 Jo s mi mi mi mi 



there being no contribution along in view of Eq. ([3.41 ). We notice that the sum rule 



is saturated by the pion pole strength alone: the contribution (in curly brackets) to the 
integrated spectral density along the cut vanishes. This just reflects the fact that the PV- 
regulated longitudinal axial density of the ENJL model again includes exactly the correct 
amount of negative density at large s for this sum rule to be obeyed. 

The sum rules that involve the isoscalar-scalar minus isovector-pseudoscalar and 
isovector-scalar minus isoscalar-pseudoscalar densities depend on the analytic behavior of 
the polarization differences 

fjSoSof 8 ) - f[ pp ( s ) = n«f ( g ) -HnfO) -4f 2 F P (s) \s\ -> oo 

u - {s) LL " {S) (i-2G 1 n^o( s) )(i_ 2Gl n^( s )) 4/ ^ Hsj ' |s| °°' 

(3.44) 



nc 



obtained from constructing Eq. (|2.27 ) minus Eq. (|2.91 ) and 



nff( S ) - n^( S ) = -4/ 2 F P ( S ), (3.45) 



on subtracting the Eq. ( |2.3(JaD from Eq. ( ^.30c| ). This is done after formally upgrading the 



proper polarizations to include non-chiral contributions too. As these effects will actually 
fall away to the order to which we are working, this upgrading is inessential, except for 
n^ c p . One notes that both differences behave asymptotically in the same way as the axial 
vector minus vector polarization of Eq. ( |2.35| ). Explicit forms of the proper polarizations 
appearing in these expressions have been given previously: n p ° p ° in Eq. ( g30g ), Tl s ° s ° = Tl ss 
in Eq. (|2~30cD , and n p c p in Eq. ( ^90|) . 
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In contrast to Eq. ( 3.45Q , the difference in Eq. (|3.44 ) contains poles at s 



ni 



2 



4m 2 + 0(m) and s = m 2 + 0(rh 2 ), given by the zeros of the two denominators, that 
determine the masses of the a and it meson modes. This structure comes about from 
summing all proper polarization diagrams linked by the interactions that are present in the 
isoscalar and isovector channels, a mechanism that is absent in Eq. (|3.45|) . The residues at 
these poles are given by the negatives — C a and — CV of their spectral strengths in the usual 
way: 

p s (s) apole = CJ(s-ml), CV = f|| (3.46) 
p P (s) npole = CJ(s-ml), Or = f|f. (3-47) 

The sum rules that are still outstanding now follow without difficulty: for sum rule ( p.8a|) , 
only the residue of the integrand at zero contributes to the contour in Fig. [15| (that is suitably 
indented to avoid both the a and the 7r poles), since the residues at these poles just cancel 
against the contributions from p s (s) apo i e and p p (s) 7rpo / e , and there is no contribution along 
Coo. The same remarks hold for evaluating the sum rule (|3.8b|) , except that the meson 
poles are absent. The contribution of the full scalar polarization n 5o5 ° to the residue in the 
first case is always of order m 2 IT 5o5o (0) ~ 0(rh 2 ), and can thus be neglected to the order 
0(m) to which we are working. By contrast, the pseudoscalar piece already contributes at 
the 0{m) level, m 2 II pp (0) ~ 0(rh). This difference arises due to the special nature of the 
quasi-chiral pion mode (mass) 2 which itself vanishes linearly with m, while the sigma mode 
(mass) 2 approaches a constant. Thus 

I" -m\p s {s) - p p (s)) = mm s ^(0) - ^(1 - ™) = -flml + 0(m 2 ), (3.48) 
Jo s ZGi m 

where the last step follows from GMOR. There is also a much more direct way of evaluating 
this sum rule using the relation 



f4{8) = ^P P (s), (3.49) 
that is readily established from the definitions of p$ and p p (s). One then finds that 

roo poo 

—rh 2 p p (s) = ds(pQ (s) 7rpoie + Po(s) cut ) 

= C aL + {- C aL + m 2 / 2 } = mlfl (3.50) 

after noticing that there is a contribution from the contour in this case, in view of 
Eq. ( p.41| ) that just cancels against the residue of (l/4)Ilf.^ c at the pion pole it contains. 
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Thus the contribution of the integral along the cut (in curly brackets) is zero, the entire 
integrated strength being determined by the pion pole in this way of calculating the integral. 
Taking the negative of this result and dropping the scalar density contribution for the same 
reason as before, we regain the sum rule ( 3.481 ). 

Exactly the same arguments hold for obtaining the sum rule in fl3.8b|) . Since there are 
no meson poles in this case, the value of the integral is given by 

/ -^ 2 (p*( s ) - p p (s)) = -4m 2 / 2 F P (0) = -Am 2 f p . (3.51) 
Jo s 1 ' 

Equating this with —Im^l-j in conformity with Eq. (|3.8b|) in order to define I7, we have 



2m 2 f 2 

h = ~ (0.5 to 1.4) x 10" 3 , (3.52) 



7T 



depending on whether the ENJL or NJL parameters are used from Table |. We have used 
the GMOR relation to fix m at 3.0MeV or 4.5MeV respectively, taking = 135MeV. The 
empirical value is ~ 5 x 10 -3 0. Finally, 

/ ds(p s (s)- p p (s)) = (a-a) + {-a + a + 4m(^)} = 4m(^) (3-53) 
Jo 

and 

ds{p s (s) - p p {s)) = 4m(i/nl>)> ( 3 - 54 ) 



using by now familiar arguments. Note that there is a contribution to both integrals coming 
from the contour that in fact determines their common value. 

In summary, we have extended the previous study [17]] by showing that all eight sum 
rules of interest to chiral perturbation theory are obeyed by the PV regulated ENJL spectral 
densities, except for modifications to the three zero moment sum rules that assume a value 
proportional to mfyip) instead of zero. Since the values of these sum rules only involve 
combinations of f n and m^, they continue to be obeyed (with /„. — > f p and m n — > m p ) by 
the PV-regulated NJL spectral densities, which do not contain resonances of any kind from 
the p — at sector. Thus, while the vector and axial vector modes produce resonant spectral 
density functions, such resonances are not a prerequisite for the sum rules to hold. We also 
comment on the modifications of the zero moment sum rules that occurred in Eqs. ( p,12|) , 
( |3.34j ), (|3.53| ) and fl3.54|) . Possible modifications of the original form of Weinberg's sum 
rules, Eqs. (|3.2a|) and ( |3.2b| ), have been investigated by Nieh and Jackiw PTfl , who developed 
criteria for obtaining these sum rules from a specific Hamiltonian. This in turn leads to the 
speculation (that remains to be explored along the lines of the Nieh- Jackiw approach for 
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example) that the modifications found here might perhaps be examples of a more general 
property associated with the dynamical breaking of chiral symmetry with mass generation 
in the groundstate, and which would hold for any Hamiltonian that energetically favors the 
formation of a finite condensate. 



IV. SUM RULES AND THE ELECTROMAGNETIC PROPERTIES OF THE PION 
A. Pion electromagnetic form factor, radius, and polarizability 

1. Form factor and pion radius 
The electromagnetic form factor, F^q 2 ), of the pion is given diagrammatically by Fig. [16. 



The diagram includes the p° exchange renormalization of the isovector component of the 
quark electromagnetic current as in Fig. [I]. This modifies the jqq vertex, shown as an open 
circle, as follows, 

- zeJ^ = -ie{\jl + JJ) - -ze[e qlfl + il#V) \v™(<?)tJ^\. (4.1) 
The charge carried by the quarks has been denoted by e q = ~(g + r 3 ), in units of e. The 



corresponding analytical expression for F w (q 2 ) has been derived as Eq. (|A74|) in Appendix 
A, section 5: 

F n (q 2 ) = {g A F P {q 2 ) + (1 - g A ) + (1 - g A ) 2 ±^ V T V {q 2 )}F v {q 2 ). (4.2) 

If the composite structure of the Goldstone n meson is suppressed by setting the contents of 
the curly brackets to unity, one obtains the vector dominance model form: F n (q 2 ) « Fv(q 2 )- 
On the other hand, setting g A = 1 returns one to the NJL expression Fp(q 2 ) for the form 
factor. The pion charge radius is identified from the low q 2 behavior of F v (q 2 ): 



- 1 + ^ - ^ + ^1 + • - 1 + \f(rl) + • ■ ■ (4-3) 



The contribution from the derivative of the vector form factor to this expression has been 
displayed separately as 

2^(0) = (1 - 9a) I ^9 Am 2 = l/m 2 V} (4.4) 

to show that it has exactly the structure of the vector dominance model when expressed in 
terms of the p meson mass my as given by the low energy expansion of the ENJL model. 
We thus obtain the following expression for the charge radius of a chiral pion, 
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w^^w-i*^- 1 '- (4 - 5) 

For qa = 1, this reduces to Tarrach's expression |48| for (r 2 ). We have introduced the 
definition of the Z 6 CHPT coupling parameter in the last step. As a case in point, we remark 
that inclusion of the leading order non-chiral contributions to F^^q 2 ) would add a chiral 
logarithm of the form 

~16^ ln(m ' /4m2) (46) 
to the middle term of the above expression for (r 2 ). This term, which arises from the 
leading non-chiral contribution of the dumbbell diagrams to F n (q 2 ), has exactly the correct 
coefficient to combine with Zg to give the prediction 

k + InK/V) = 1 + 12g 2 A + 8vr 2 4(^^) (4-7) 

at the scale \x = 2m. The right hand side of this equation gives the ENJL value for the 
renormalized coupling constant Zg up to a conversion factor ||. However, there is no need 
to include this factor to compute Zg, since l 6 is the physically interesting quantity and this 
may now be computed directly. Numerically, the calculated radius comes to 0.36 fm 2 for 
our standard parameter set, and leads to the value Z 6 = 16.2. The experimental value is 
(r 2 ) = 0.44 ± 0.03 fm 2 p9|j50|] , which in turn leads to the empirical value Zg = 16.5 ± 1.1. 



The behavior of the ENJL form factor is compared with the experimental data |5T| for 
\F n (q 2 )\ 2 in Fig. |17] at both space-like and time-like q 2 . Given that no parameters have 
been re- adjusted, one notes that the overall agreement is fair. Thus is particularly true at 
small space-like q 2 where the behavior is quantitatively the same as the CHPT prediction 
|512|] (which, however, uses the experimental pion radius as input). This agreement simply 
reflects the reasonable value of (r 2 ) = 0.36fm 2 for the ENJL prediction of the pion radius. 
For larger q 2 , and more especially in the time-like region, the same comments as were made 
in connection with the continuum structure of the vector and and axial vector spectral 
densities should be borne in mind here too: The continuum behavior of the calculated \F V \ 2 
for q 2 > 4m 2 in Fig. O is determined by the qq decay channel, not the physically expected 
p — > mi channel. Thus the ENJL form factor embraces contributions in the time-like region 
from a decay channel that would not be available in a confining theory. 



2. Polarizability of charged pions. 



Using the current algebra sum rule of Das, Mathur and Okubo [53~| for the axial structure 
constant that enters into the tc — > ez/7 radiative decay matrix element, Holstein has derived 
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the following expression for the electric polarizability of charged pions f22 

a , o, a r°° ds 



a E 



, (rl)-^- r%{ P X{s)-pt{s)). (4.8) 



The first term corresponds to the "classical" contribution from a charged particle of radius 
squared (r^) [Q- Taking the pion radius squared from the Eq. ( [4.5|) and the value of the 
integral from Eq. ( |3.38| ), one obtains (we express the result directly in terms of the ENJL 
l^s, and take their difference from Table [D]) 

cle = —jAr^ NJL - l! NJL ) = Ah— = 3-4 x 10-fm 3 , (4.9) 



if gA = 0.75. In computing the numerical result in Eq. ( f4.9| ), we have also used m n = 135MeV 
and /„- = 93MeV in addition to a ~ 1/137 as input. 

The first expression coincides exactly with the form of the result for ag that is obtained by 
treating the forward scattering of soft photons by pions in chiral perturbation theory at the 



one-loop level |P3|| . There are three additional noteworthy features about this expression: 
(i) since the radius contribution and the contribution from the integral over the density 
functions occur with exactly the same numerical coefficient when expressed in terms of the 
chiral coupling constants, all factors common to both Z 6 and Z5 cancel in the difference. 
This feature leads to a complete cancellation of the contribution from the p meson degrees 
of freedom, a feature that has also been obtained by Holstein [^] via another method, 
that of saturating the sum rule with delta function distributions in conjunction with the 
KSRF relation, (ii) the quark mass m has fallen away due to this cancellation, and (iii) 
the quenching effect of the a\ meson degree of freedom only enters via the axial form factor 
gA- Thus Eq. ( |4.9| ) is a prediction of the ENJL model that only contains physical quantities 
that are known in principle. In particular, it is a prediction for an upper limit on ag, 

a 1 

ole < ^T77 ~ 6 x 10 " 4 fm3 - ( 4 - 10 ) 

The result ( [4.9|) , and the upper limit in particular, can be compared directly with the 
prediction of chiral perturbation theory at one loop that uses the empirical values of the 
coupling constants I5 and ^6 to obtain [f22| .f23[| 

for the polarizability of the charged pion. One can also compare our result with the estimate 
a E = (5.6 ±0.5) x 10 _4 fm 3 that has been obtained from Eq. ( |4.8|) by using the experimental 
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pion radius as input, and estimating the integral involving the vector and axial vector 



densities using QCD sum rule methods |55 



The experimental value of a# is still very poorly known. The data analysis for the 
combinations ag± /3m, where (3m is the magnetic polarizability, yields values for from 
2.2 ± 1.6 |2|J2|] through 6.8 ± 1.4 |26| to 20 ± 12 |27||, in units of 10" 4 fm 3 . In addition, 



these numbers depend somewhat on whether the constraint of good chiral symmetry |56| 
ole + Pm = is implemented or not. Thus one may conclude that the chiral perturbation 
theory prediction and the ENJL model sum rule prediction agree qualitatively with each 



other, and that neither result agrees with the experimental data to date. We refer to [28| for 



an overview of the present state of experiment and comparison with theoretical calculations. 

B. Pion electromagnetic mass difference 

In a historically important paper from the point of view of current algebra techniques, 
Das, Guralnik, Mathur, Low and Young [29] exploited Weinberg's two sum rules, Eqs. (|3.2a|) 



and ( p.2b| ) for massless QCD to obtain an expression for the pion electromagnetic mass 
difference in terms of vector and axial vector spectral densities. These authors showed 
that the mass squared difference between the charged and uncharged chiral pions (therefore 
having common rest mass zero) may be written in the form 

Ami = ml, - ml, = X*f? j |^ £ ^0*00 " (4-12) 

to 0(a) in the electromagnetic coupling a = e 2 /4n « 1/137. By saturating the Weinberg 
sum rules with low-lying vector and axial vector mesons, they obtained the estimate Ami ~ 
(3a/27i)m 2 p In {m 2 ai /m 2 ) = (37.85MeV) 2 in terms of the p and a\ meson masses, taking m p = 
770MeV and the Weinberg ratio m ai = y2m p for the a± mass. The current experimental 
value is AmJ = 1261MeV 2 =(35.51 MeV) 2 . This calculation represented one of the first 
successful applications of QCD sum rules. 

By contrast, the ENJL expression for the spectral density difference appearing in 
Eq. ( p.28 ) allows us to evaluate the Das et al. expression without making any further 



assumptions as to the resonance content of the spectral density functions. Since, according 
to Eq. (|3.30|) , the required integral over s is given by the residue f(s) at q 2 instead of zero, 



where f(s) = f 2 F P (s)F v (s)FA(s), one immediately has the result 

d 4 q 1 



= 3*e 2 / 4\l~,F P {<l 2 )Fv{<l 2 )Uq 2 )- (4-13) 
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A formula essentially equivalent to this expression has also been derived using a QCD ef- 
fective action approach ||13|| . Below, we give yet another derivation of this expression using 
Feynman diagrams. In the limit qa = 1, where Fy and Fa are both unity, it reduces 
to the expression given previously [15] for AmJ for the NJL model, while making the ad 
hoc assumption that all the mesons are elementary without internal structure of any kind, 
Fp(s) = 1, Fy 7 A ~ (1 — q 2 /m 2 pai )~ l leads one back to the Das et al. estimate. 

Futher evaluation of Am 2 has to be done numerically. However, since the product of PV- 
regulated form factors in the integrand of Eq. ( [4.1 3D goes like — 1/q 2 according to Table [VI|, 
the integral over the photon momenta diverges logarithmically, and some form of photon 
momentum cutoff is required. This contrasts with the special case of approximate evaluation 
by Das et al. mentioned above, where the asumption of delta distributions for the vector and 
axial vector densities that satisfy both Weinberg sum rules guarantees a convergent result. 
A logarithmic divergence also results when their approach is extended to cover physical (i.e. 
non-chiral) pions ||58|| . 

Performing a Wick rotation to Euclidean photon momenta, one finds 



9 3a 



f' h dQ 2 F P (-Q 2 tF v (-Q 2 )F A (-Q 2 t. (4.14) 

JO 



This transformation is permissible since the poles of the integrand all lie on the second 
Riemann sheet and thus do not interfere with the deformation of the integration contour. 
The behavior of Am 2 versus X 2 h is shown in Fig. [18]. A cutoff of X 2 h = (1.27GeV) 2 leads to 
the observed value AmJ = 1261 MeV 2 of this quantity, if we use the standard parameter set 
of this paper to compute the form factors in Eq. ( |4.13 ). This value of the photon momentum 
cutoff requires some comment. In the 0(4) regularization scheme (which we do not employ) 
the Euclidean quark momenta are limited by the cutoff Ao(4), so that the photon momenta 
are restricted by Q 2 < AAq/^ in Eq. (|4.14|) . Since one has the rough equivalence [16 



Aq/ 4 n « (21n2)A 2 ~ (2.5 GeV) 2 , this restriction is met by our photon cutoff. Furthermore 
this cutoff also agrees satisfactorily with the scale fi 2 = (0.95 GeV) 2 quoted in J0|] at which 
a smooth matching |57| between Eq. ([4.14|) (for a different set of input parameters however) 
and the short distance QCD contribution to Am 2 can be accomplished. 



C. Derivation via Feynman diagrams 

While the derivation of the expression just given for Am 2 via spectral densities is tech- 
nically very concise, the underlying physics determining this splitting is less transparent. 
We remedy this by using the diagram method. One can then verify the ENJL expression 
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Q4.13Q for AjtjJ explicitly by identifying the additional set of Feynman diagrams for the pion 
polarization that include single photons in intermediate states. This method has already 
been described in [[TJJ for the case of the minimal NJL model. The diagrammatic method 
offers some advantages over the sum rule approach. The electromagnetic gauge invariance 
can be made explicit, as well as the consequences of the residual Ua{\) x Z7y(l) chiral sym- 



metry as expressed by Dashen's theorem |TJJ], that keeps the neutral chiral pion mass at 
zero in the presence of its own internal electromagnetic field. It confirms that the Das et 
al. formula in conjunction with the ENJL spectral densities is correct to 0(aN c ). Also, it 
is straightforward to extend this method to cover physical pions carrying non-zero mass in 
the absence of the electromagnetic interaction. An attempt to extend the approach of Das 
et al. directly to cover the case of non-zero pion mass is reported in ||5S|| . 

In order to identify the relevant electromagnetic diagrams for calculating the electro- 
magnetic mass splitting of the pion, we include the electromagnetic interactions via minimal 
coupling in the ENJL model, and study their effect to 0(a). The additional diagrams are 
obtained in the usual way by dressing the skeleton pion polarization diagrams with single 
photon self-energy and vertex corrections. The diagrammatic calculation of AmJ to 0(a) 
starts out with the observation that only those diagrams involving the exchange of a single 
photon that are different in the isovector and isoscalar channels will contribute. This im- 
mediately excludes all quark self-energy diagrams in the pion polarization where the photon 
line dresses either a quark, or an antiquark line. Such diagrams (which are however essential 
to maintain gauge invariance) contribute equally to the charged and neutral pion channels 



and thus cancel in the mass difference calculation. They are shown in Fig. [19] for later 
reference. The diagrams we thus require are of two types: (i) scattering diagrams where 
the quark- antiquark pair making up the polarization loop interact via photon exchange in 
intermediate states, and (ii) meson-pole, or "dumbbell" diagrams where the n + ai,p and 
a\ meson modes are exchanged together with a photon between two 'jtttt vertices. These 
diagrams are presented in Fig. In addition, p° exchange causes a renormalization of 
the isovector component of the quark electromagnetic current as in Eq. ( |4.1| ). The modified 



vertex itself is again indicated as an open circle to which the photon is attached as in Fig. [16. 

Both sets of diagrams in Fig. are vertex corrections to the basic one-loop irreducible 
diagram that determines the pion polarization to 0(aN c ). At first sight, the appearance of 
such two-loop diagrams containing mesons in intermediate states is surprizing, since they 
formally contain higher powers in the meson quark coupling constants than the single loop 
diagram. However they are of the same order 0(aN c ) as the single loop contribution to 
the electromagnetic self-energy of the pion, and therefore of 0(a) relative to the pion self- 
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energy in the absence of electromagnetic interactions. This comes about since the meson 
propagators scale like 1/N C . This is a special feature arising from the linked nature of the 
solutions for the quark and meson self-energies in the NJL model. 

Let Ii^ M (q 2 ) denote the contribution from the sum of diagrams in Fig. 20 in the channel 



of isospin projection T3. Here T3 equals ±1 for ^ and for it . Then Am\ is given by (15 



Ami = -<4 9 {n± M (0) - U° EM (0)} (4.15) 



one can 



for chiral pions. Using an extension of the arguments presented in detail in [15[ 
again show that this expression is fully gauge invariant. We can thus describe the photon 
propagation in any convenient gauge. We use the Feynman gauge, 

-iD^) = -i 9 —. (4.16) 

The translation of all the Feynman diagrams in Fig. ^ is deferred to Appendix A, section 
6. Here we simply summarize the results for the two classes of diagrams: 

(i) Scattering diagrams. The scattering diagrams give rise to a contribution to the difference 
( |Og ) that is of the form [Eq. ([MOD I, 

{H± M (0) - U m M} scatt = -*e 2 Jr / J^Fptf) 

[ 7^^(? 2 )W) ■ (4-17) 



m 2 J (2tt)V 

(ii) Meson pole, or dumbbell diagrams. By contrast, the dumbbell diagrams only contribute 
to the charged channels. There is no contribution for T 3 = 0, due to the electromagnetic 
vertices they contain. One finds for the sum of the tt + at, contributions that [Eq. (|A88 )] 



n^ aL) (0) = f j^Frtf) (4.18) 



w? J (27r) 4 g : 
while 

<>(o, - 3 - 2 l / j^k^ 2 ^ 1 WMw) } (419) 

gives the contribution from the a\ dumbbell diagram [Eq. (|A94|) 1. In obtaining the latter 



result, we have inserted the specific form of the a± propagator into Eq. ( |A94j) . Finally we 



note that there is no contribution from the analogous dumbbell diagram with a p meson 
exchange replacing the a\ vector meson exchange in the chiral limit. The reason for this is 
given in more detail in Appendix A. 
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We now add Eqs. ( f4.17|) through (|4.19|) . The 7r + dumbbell contribution cancels 



against a similar term in the scattering contribution, and we are left with a result for the 
self-energy difference required in Eq. ( |4.15| ) that may be written as 

-iAM - nWO)} = 3</^/^M 5 *)W){l - / + G ;|gg 2) }. 

(4.20) 

At first sight, this expression gives the impression that the contributions from the p° degrees 
of freedom enter quadratically through the vector form factor Fy(q 2 ). This impression is 
false. Using the Ward identity of Eq. ( 2.37]) in reverse, one sees that 



r SG 2 f}F P tf) 1 _ l + 2G 2 nr(g 2 ) _ , , y (2) ( , . 

I 1 " 1 + 2G 2 I#V) / " l + 2G 2 n^(g 2 ) ~ ^ V v(9 } ' ( } 

Thus one power of Fv(q 2 ) is removed. Furthermore, since g 2 qq = gp/g A , the strong inter- 
action physics part of the factor in front of the integral becomes f 2 g 2 /{m 2 gA) = 1/gA 
that supplies precisely the correct factor to convert Ga{q 2 ) into the axial form factor 
F A (q 2 ) = G A (q 2 )/g A . Hence 



Ami = 3te2 I ^4 Fp{q2)Mq2)FA{q2) 



We thus recover the same expression for Am 2 as given in Eq. ( 4.13| ) by a completely different 



route. All that remains is to specify the regularization procedure. 

Having obtained expressions for the sub diagrams making up the pion polarization, it 
becomes a simple matter to identify the underlying physics of the pion mass splitting in 
more detail. To this end, let us introduce the electromagnetic self-energies E Ts for the 
charged (T3 = ±1) and neutral (T3 = 0) pions. This self-energy, which has the dimensions 
of a (mass) 2 , is obtained from the total pion polarization by including the nqq coupling 
constant ig nqq at each external pion vertex: 

S Ta = (^ 99 ) 2 n| 3 M (0). (4.22) 

Consequently, a negative pion polarization corresponds to a repulsive self-energy, and vice 
versa. Now this n|f M (0), and therefore E Ta , consists of the sum of all the polarization 
diagrams appearing in Fig. [H| plus Fig. |2(| Thus we have reinterpreted this II^ 3 jy(0) slightly 
in order to avoid introducing yet another symbol. Let us write out this sum for the two 
charged states of the pion separately. In doing so, we notice that the sum of diagrams from 



Fig. [L!|, that are all in the nature of self-energy insertions on the quark and antiquark lines, 
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are common to both charge channels. This is so because the quarks and antiquarks enter 
symmetrically. This is not true for the diagrams in Fig. |20|. Here diagram (a) of this figure 
represents the interaction via photon exchange of the qq pair in the polarization loop. For 
the like-flavored quarks making up the 7r°, this is an attractive interaction; otherwise it 
is repulsive. Finally the dumbbell diagrams in Fig. ^(b) only contribute in the charged 
channel due to the electromagnetic vertices that they contain. We continue to work in the 
chiral limit. Then the electromagnetic (mass) 2 of the 7r° is 



r(°+ 6 )frv\ _l rr( c + d )/ 

(scatt 



2 -(i 9mq r ngJ'(o) + ngJ?(o) +fe ? ) 2 {n° M (o)} f 



:: - ■ , t , , 

Fig. 19 Fig. |20|(a) 

= Z self + Y? scatt . (4.23) 

Further progress is made by calculating the actual values of all these pion self-energies. This 
is completed in Appendix A7. The scattering contribution represented by the last term is 
known from Eq. (|A75|) . It is negative as anticipated (in the Feynman gauge), and exactly 
cancels the quark self-energy contribution given by the sum in the first bracket. One finds 
by direct calculation the result given in Eq. ( |A108|) that 

Zseif = -Kcatt = ieWj (^{ 4tr ( e ?) + \(*vtf) ~ 1)}^(9 2 )- (4-24) 

Thus m^o = 0: the neutral pion does not acquire any electromagnetic mass at all from 
its own internal electromagnetic interactions in the chiral limit. The increase in its mass 
due to the electromagnetically heavier quarks it contains is compensated for exactly by the 
attractive interaction between these quarks. Notice that this is a gauge invariant statement. 
While the individual terms appearing in Eq. ( |4.23j ), are not all gauge-invariant, their sum 
is. More formally this result, which is just Dashen's theorem in operation, has to follow, 
once one realises that the inclusion of the photon gauge field via minimal coupling breaks 
the U(l)v x SUl(2) x S77#(2) chiral symmetry of the original ENJL Lagrangian in Eq. ( |2.8|) 
down to Ua{^-) x LV(1)- The massless 7r° is simply the Goldstone mode realization of this 
residual symmetry. As we have shown, this result is satisfied in an entirely non-trivial way 
by the ENJL model. It is also a proof that our set of polarization diagrams to 0(aN c ) is 
correct and complete. 

By contrast, the charged chiral pion picks up a non-zero electromagnetic mass: 



mU = (W n^(o) + i4 c + d) (o)l + (^) 2 {n| M (o)} 



:catt 



Fig. 19 Fig. gCj(a) 

+(^) 2 fn^ aL) (o) +nK(o)l, (4.25) 



Fig. 2* fc) 
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or 



m n ± - E se i f + S scatt + E dumbbells (4.26) 

using a self-evident notation for the various contributions to the self-energy of the charged 
pion. In this expression, both the quark self-energy contribution as well as the scattering 
contribution, are positive, the latter being so because the EM interaction between qq pairs 
of different flavors (ud or du) is repulsive. Thus the last term arising from the sum dumbbell 
diagrams must be negative, and, as we will see, strongly so, in order to cancel out most of the 
repulsive self-energy associated with the first two terms. We naturally regain Eq. (|4.15|) again 



by forming the difference m 2 ± — m^. As we have seen, only the charged pion contributes 



to this difference. In Fig. 21, we display the various contributions to the electromagnetic 



(masses) 2 of the 7r° and 7r ± , using the same parameter set as for evaluating Eq. ( 4.13| ). In 



examining this figure, it should again be borne in mind that the individual shifts are gauge- 
dependent (we have used the Feynman gauge). Only the total electromagnetic (masses) 2 , 
and therefore their splitting, are gauge-invariant quantities. 

These results only hold in the chiral limit of vanishing pion mass. Their extension to 
the non-chiral case of physical pions requires re-evaluation of the same set of diagrams with 
the external pion placed on its mass shell at k 2 = m 2 already in the absence of its own 



internal electromagnetic field. This is a tedious task, even for the minimal NJL model |L5 



and no less so in chiral perturbation theory [55]. However, the main conclusions reached 
by including non-chiral effects [TJJ carry over without change into the present calculation, 
and can be restated quite simply. The expression (|4.12j ) accounts for the major part of the 
electromagnetic mass shift. In addition, the assumption of a finite current quark mass in 
Eq. (|2.8| ) breaks the chiral symmetry explicitly. This gives the pions of the model a finite 
mass m 2 ~ rh as per Eq. ( ]2.100| ). The sub leading (in m 2 ) contribution to the electromagnetic 
splitting from this source comes from expanding the scattering and dumbbell diagrams in 
Fig. |20] in terms of m 2 . This gives rise to 0{am 2 T ) correction terms in the neutral channel, 
but an additional chiral logarithm ~ am 2 In m 2 appears in the charged channels. The latter 
term, which is identical to the correction found in CHPT |33], arises specifically from the 
pion dumbbell diagram in Fig. It dominates the non-chiral correction. Since the chiral 
limit of this diagram, as given by Eq. ( |4.18| ), is the same for both versions of the NJL model, 



we can take the precise form of the chiral logarithm correction from the former calculation 



[15 1 and write 



{ Am 2 } non _ chiral = Am 2 - ^ ln(m 2 /m 2 ) + 0(am 2 ) + • ■ • (4.27) 
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for the ENJL model too. The chiral logarithm correction in this equation coincides exactly 
with the result obtained in chiral perturbation theory [R|[] (at scale /i 2 = m 2 ). Numerically, 
the chiral logarithm and 0(am 2 ) terms constitute typically 10 % and 1 % of the leading chiral 
mass splitting term. Unequal current quark masses, or TT°r] or rf mixing in the Ul(3) x Ur(3) 
flavor extension of the ENJL model (including the 't Hooft term |3f| to remove the unwanted 
Ua(1) symmetry), are additional sources of chiral symmetry breaking that contribute to 
Am 2 . Both effects contribute at the few percent level. One knows this from a direct 



estimate [15 of the up-down current quark mass difference contribution ~ (rh u — mi) to 



Am 2 , and from a rather general discussion given in |Kj for the mixing effects. Bearing these 



results in mind, one is led to regard the mass of the charged pion, and especially the 7r ± — 7r° 
mass difference, as an essentially electromagnetic phenomenon. 

V. SUMMARY AND DISCUSSION 

In this paper, we have established that the two sum rules of Weinberg for massless 
QCD, plus the six sum rules introduced by Gasser and Leutwyler in connection with chiral 
perturbation theory, are all obeyed by spectral densities that are generated by the Nambu- 
Jona-Lasinio model Lagrangian in its minimal, as well as in its extended form, and which 
includes vector and axial vector degrees of freedom. The only exceptions to this statement 
occur for the zero moment sum rules [Eqs. (|3.34 ), ( |3.53j ) and ( [3.54 )1 that acquire a finite 



value proportional to the quark condensate density that represents the order parameter of 
the model. We suggest that this finding might perhaps be an example of a more general 
feature of any chiral Hamiltonian that energetically favors a chirally broken groundstate 
with a finite quark condensate. 

The NJL model in either form is both non-renormalizable and non-confining. Thus the 
statements regarding sum rules are subject to stating what sort of regularization precedure 
has been followed. We show in particular that the sum rules hold for spectral densities com- 
puted from polarization amplitudes which have been regulated according to the Pauli-Villars 
prescription. As might be anticipated, the PV-scheme introduces an unphysical behavior 
into the spectral densities at high energies, where these densities become negative instead 
of remaining positive definite. Nevertheless this happens in such a way that the sum rules 
continue to be saturated. We interpret this rather surprizing conclusion to be a consequence 
of the underlying chiral symmetry properties of the Lagrangian and its faithful representa- 
tion under the PV-scheme. As a corrollary to this statement, we find that the violation, due 
to PV-regularization, of the positive definite nature of the model densities becomes less and 
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less important for the sum rules involving higher and higher inverse moments, such as the 
sum rule determining the pion polarizability, for example. 

While the sum rules alluded to above also hold for the spectral densities belonging to 
the minimal NJL model that does not contain any explicit vector meson modes, this is 
perhaps less interesting from the point of view of physical applications. The sum rules 
involve inter alia explicit reference to both vector and axial vector densities, and so it is of 
direct physical interest to consider such degrees of freedom explicitly. With this in mind, 
we have revisited the problem of incoporating these modes into the NJL model in order 
to determine the behavior of the associated spectral densities and their masses. The basic 
effects of these modes are (a) to renormalize the pion decay constant and nqq coupling 
constant in a complementary fashion such that the Goldberger-Treiman relation is retained, 
and (b) to concentrate the vector and axial vector strengths of the one-loop spectral densities 
of the minimal NJL model into two well-defined peaks at low energies that lie just above the 
qq continuum threshold. This redistribution of strength takes place in such a way, however, 
that all the sum rules continue to be obeyed. If the positions of these peaks are taken to 
define the p and a\ meson masses of the ENJL model, their numerical positions are found 
to be in qualitative accord with the experimental masses of the physical p and a\ mesons, 
suggesting that this identification has physical meaning. 

We have combined the sum rules techniques developed in this paper, together with other 
calculated information to identify ENJL expressions for the empirical coupling constants 
introduced by Gasser and Leutwyler into chiral perturbation theory, to evaluate the current 
algebra expression of Das et al. for the charged to neutral electromagnetic mass splitting 
of the pion, and to compute the charged pion's polarizability from Holstein's sum rule. In 
extracting the coupling constants of chiral perturbation theory, it has been borne in mind 
that the physically relevant scale-invariant constants, k, are not well-defined in the chiral 
limit because of the chiral logarithms they contain. However, the combination ^ + In {rn\j p?) 
is well-defined, being proportional to the renormalized l\ of Gasser and Leutwyler, and it is 
this combination that one can extract from the ENJL model. In fact, it is possible to go 
further and extract the Zj's themselves due to the following circumstance. One can show that 
the leading order non-chiral corrections to the ENJL amplitudes are precisely in the form of 
chiral logarithms that occur with the correct coefficients to reproduce the combination given 
above (at the scale of the sigma meson mass of the model). The appearance of the correct 
chiral logarithm has been demonstrated explicitly for the case of 1$, for example, which must 
be identified from a calculation of the pion radius. One is thus fully justified in solving for 
the U themselves from the ENJL predictions of the l\. This we have done and have obtained 
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reasonable agreement with the empirical values of the l^s as set out in Table |HI . 

The sum rule calculation of the pion polarizability made direct use of the ENJL 
model generated spectral densities, as well as the ENJL expression for the pion radius. 
The resulting expression for oe, which coincides exactly in form with the CHPT result for 
this quantity, only depends on the difference Iq NJL — l^ NJL . This allows common terms in 
the two ENJL values for the coupling constants to cancel, and leads to an upper limit of 
a/(8ir 2 m n f 2 ) ~ 6 x 10 _4 fm 3 on in terms of the physical mass and decay constant of 
the pion. This is an important result, since it represents a prediction of the ENJL model 
approach that may be compared directly with the chiral perturbation theory prediction 
(a C E Pt = 2.7 x 10" 4 fm 3 ), as well as the result of Ref. |5j| (a Q E CD = 5.6 ± 0.5 x l(T 4 fm 3 ) that 
is based on an evaluation of Holstein's sum rule using approximate QCD spectral densities 
and methods. We also showed that the limiting value of cue quoted above is reduced by g\, 
the quark axial form factor squared, in the presence of axial vector mesons accompanying 
the pion. Such meson clouds lead to a less polarizable system: for g& = 0.75, one finds 
ctE = 3.4 x 10~ 4 fm 3 that is nicely bracketed by the chiral perturbation theory and QCD 
sum rule predictions. By contrast, the p mesons do not change the pion polarization at 
all. Their contributions cancel out completely, in complete agreement with the findings 



of Holstein p2| , who employed delta function spectral densities together with Weinberg's 
vector meson mass relation and the KSRF relation. 

Finally we must caution that the the limiting value on «j quoted above is subject 
to a multiplicative correction of the form [1 + 0(m 2 ) + 0(1/N C )] where the additional 
modifications arise from neglected non-chiral effects of 0(m 2 ) in the sum rule itself, and 
neglected chiral corrections of 0(1/N C ) in the evaluation of the spectral densities that enter 
the sum rule. The calculation of such higher order correction terms within the framework 
of the ENJL presents an important but formidable challenge, especially since the chiral 
perturbation theory results for ole to two loops, consisting of the calculation and summation 
of over 100 diagrams(!), are now available [p3| for direct comparison. 
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APPENDIX A: REGULARIZATION OF FEYNMAN DIAGRAMS 



In this appendix, we discuss the regularization of the various Feynman loop diagrams 
that have appeared in the main text. 

1. Gap equation for the light quark masses 

The gap equation in the chiral limit rh = that generates the quark mass self-consistently 
in both the NJL and ENJL versions of the model Lagrangian in the mean field, or Hartree 
approximation has been written down in Eq. ( [2.34 ), and its extension to include a non- 



vanishing current quark mass in Eq. ( 2.89Q . Either equation contains a single quark loop 



that diverges quadratically. Regulating this loop according to the Pauli-Villars prescription 
for the reasons explained in the main text, we get 



m — m 



V a J ( 2 7r)V-M2J 

(!6GiiV c m){ j-y 2 £ C a Ml In — }, (Al) 



provided that 



£C a = and ]TC a M a 2 = 0, (A2) 

a=0 a=0 



that may be satisfied by the standard choice [fTS 



C a = (1, 1, -2) and M 2 a = m 2 + a a A 2 , a a = (0, 2, 1), (A3) 

where A is the regulating parameter of the PV scheme. We designated the solution of 
Eq. (|Al|) in the presence of rh as m = m* in the main text. One has the useful relation 

EC ^_4 = ^ (m _* ) = _£<M, (A4) 
„ m z GiJM c m JS C m 

after using the mean field relation m — rh = —2G\(ipip) between m and the condensate 
density (ipijj) of the quarks. 



2. The one-loop integral I(q 2 ) 
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a. Evaluation 



The integral I(q 2 ) given in Eq. ( |2.33| ) appears in the pion weak decay constant f 2 = 
—4:N c m 2 iI(0), in the pion form factor Fp(q 2 ) = il(q 2 )/il(0), and thus the pseudoscalar 
polarization loop as given by Eq. ( |2.30a|) . To evaluate it we introduce the standard Feynmann 



parametrization |18| to find 

d A p 



(2vr) 4 \{p + q) 2 -m 2 ][p 2 - m 2 ] 

poo 

/ dalAWn 2 - q 2 a{l - a)}, (A5) 
Jo 



(47T) 



where 



h(q 2 )= / ^P e -ipl™ 2 -zMl-<*)\. ( A6 ) 
Jo p 

The last integral is logarithmically divergent. We regulate it using the Pauli-Villars [PV] 
prescription again. Denote the corresponding PV-regulated quantity with a bar. Then 

dp 



h(q 2 ) = lim ^2 C a I" ^e^ 1 *-^ 1 -^ = - £ C a \n[M 2 - q 2 a(l - a)]. (A7) 



Hence 



z/( S ) = z7(0) + JL Y,C a £daHl ~ - «)], (A8) 

where 

1 M 2 
(4vr) 2 ^ m 2 

after calling q 2 = s. The remaining integral is elementary. For space- like momentum 
transfers squared — oo < s < 0, it is given by 

/ 1 daln[l--^ra(l-a)] = 2J(z)-2, (A10) 
Jo m 2 

where 



J( z ) = Jf-coth- 1 y / 7= yi- ^ lnfv^ + v / T 3 ^] (All) 

~m^ (1 _^__L + ... ) + ln2 _ ( i + i ln2) I 

+ (^-^ln2)^ + ---, |z|-oo, (A12) 
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after setting z = s/4m 2 and 



/=1- 



The PV-regulated version of il(s) for space-like q is thus 

= 7^[|> ^ + ^ coth- ^- 2)} 



(A13) 



(A14) 



where the PV-subscript signifies the replacement ]1S 



{(2^// coth- 1 y/f - 2)} PV - £ C a (2^// a coth- 1 ^ - 2) 

a=0 



in terms of the variable 



AM 2 

fa=l -■ 



(A15) 



(A16) 



In accordance with Feynman's causal prescription [IS], i/(s) has to be continued to time- 
like q 2 = s > by passing along the upper lip of the cut starting at the branch point at 
s = Am 2 and running to oo. This in turn means that the square roots in logarithmic function 
determining J(z) have to be chosen such that 



— = \ z \l/2 e i(#-v)/2 



|1 _ z \W e W-*)fi 



< (f)< 2tt 
< 9 < 2tt 



(A17) 
(A18) 



on the "first", or physical sheet of J(z), in terms of the complex variables z = \z\ exp iip and 
z — 1 = \z — l\ exp iO. The further branch points at z — s/4M 2 , a = 1, 2, that are introduced 
into J{s/Am 2 )pv and thus il{s) by the regulating masses of the Pauli-Villars scheme are 
handled in the same way. The schedule of forms for the multivalued complex function J(z) 
can now be written down from 



J{z) 



l/2 e <(0-0/2 ln ^|l/2 e i(<^-7r)/2 + |1 _ ^| 1 /2 e i(0-7r)/2]_ 



(A19) 



Along the real axis z = x of the first sheet J(x) thus behaves like 

' (i_ I)i/2 ln ( v /3^ +v /T^) = //coth" 1 V7, -oo < x <0 



J(x) = (i - 1)V2 ^-1 ^ = cot- 1 v^J, 



< a; < 1 



I (i - i) 1/2 (MVS if 



A/7(tanh 1 V7-^f) ; 1 < a; < oo, 



(A20) 
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where / = 1 — 1/x. The appearance of an imaginary part for x > 1 is associated with the 
(unphysical) qq decay threshold that is one of the drawbacks of the NJL model. 

In order to discuss the possible poles of the vector and axial mode propagators in ENJL, 
we also require the analytic continuation, J(z) of J(z) through the cut onto the second 
sheet. This sheet is defined by letting the angles <fr an d 9 assume values in the angular 
interval — 2n < 9, <p < 0. Except for the sector 1 < x < oo along the real axis of the second 
sheet where the two functions coincide by construction, the analytic continuation J(x) is a 
different function on the second sheet. In particular, its behavior along the real axis of that 
sheet is as follows: 

J(x) — vk\J~$, — oo < x < 

J(x) = I J(x) - TTy/^f, < X < 1 



J(x), 



1 < X < oo. 



(A21) 



We note in particular that J(x) like J(x), is also real for < x < 1 but picks up an 
imaginary part in the sector — oo < x < 0. Their behavior at the origin as x — > + is 
also different: J(x) — > x~ 1//2 sin x 1//2 — > 1, but J(x) — > — nx^ 1 ^ 2 — > —oo. Thus ReJ(x) and 
ReJ(x) coincide for — oo < x < and 1 > x > oo, but differ for < x < 1, while ImJ(x) 
and ImJ(x) differ in the first interval but coincide in the latter two. All of these properties 
are essential for appreciating the analytic structure of the various polarization functions to 
follow. The real parts of J(x) and J(x) are plotted for real arguments in Fig. [F[ 

Returning to the problem at hand, we can now write down a general form for the PV- 
regulated version of H(q 2 ) for all q 2 = s as 

1 



" 00 = 7^-2 [ E Ca In + {(2 J{s/Am 2 ) - 2)} 



(47T) 



PV 



a=0 



(A22) 



where 



{(2J(s/4m 2 ) -2)} 



PV 



Y.C a {2J(s/4M 2 a )-2}, 



(A23) 



a=0 



as before. The form of J appropriate to the sheet, and interval that its argument falls in 
on that sheet, has to be selected from Eqs. (|A20| ) or ( |A21 ). In particular, for time-like 



q 2 = s > 4m 2 , the explicit form of Im[il(s + ie)] along the upper lip of the cut plane is 
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1 / Am 2 

Im[il(s)\ = -—{\ 1 9(s- Am 2 )} 

lo7r V s 



PV 



1 2 AM 2 

E^V 1 a -9{s-AMl). (A24) 



167T ^ V S 

Thus, apart from the contribution starting at s — Am 2 , the imaginary part of il(s) also con- 
tains additional contributions of different weights coming from the artificial thresholds that 
are introduced by the regulating masses of the Pauli-Villars scheme at s — Am 2 (l + A 2 /m 2 ) 
and s = Am 2 (1 + 2A 2 /m 2 ). This in turn means that Im[il(s)} need not remain negative 
definite for s large enough, as could perhaps have been anticipated from the introduction of 



indefinite metric sectors of Hilbert space by the PV- regulation procedure [18 



b. Asymptotic behavior and dispersion relations 



The asymptotic behavior of the function J[z) is given in Eq. (|A12 ) which shows that it 



diverges logarithmically Thus if the PV instruction were removed from the second factor 
in Eq. (|A22Q , then the resulting expression for il(s) would also diverge logarithmically, 
il(s) ~ log(— s), with s — > oo in the complex plane. The PV regularization instruction 
alters this behavior drastically to a 1/s convergent one as we show below. The change-over 
is, not surprizingly, related to keeping the imaginary part PV regulated as in Eq. (|A24 ) 



(which is finite and does not need regularization anyway), which then determines what 
sort of dispersion relation the expression for il(s) can satisfy in a consistent fashion. The 
change in asymptotic behavior in the present case comes about because upon PV regulating 
J(s/Am 2 ), the conditions that J2 a C a and J2 a C a M 2 must vanish, suppress both the ln(— s) 
and the (1/s) ln(— s) terms in Eq. ( |A12| ), leaving only a 1/s and the sub-leading logarithm 



^l/s 2 )ln(— s) to survive: 

1 M 2 1 M 2 

J{s/Am 2 ) = E C a J(s/ AM 2 ) ~ - - £ CJn -± + -(£ C a M 2 a In ^) 

- iln(-s/4m 2 ) (ZCaMt) + 0(l/s 2 ). (A25) 

S a 

Inserting this form of J into Eq. ( |A22| ) and using the result that J2 a C a M^ = 2A 4 to simplify 
the coefficient of the (1/s 2 ) ln(— s) term, one obtains the behavior of il{s) at large s as 

1 /If 2 A 4 1 

= ^E^ln-f - ^Hs) + 0(1/* 2 ), (A26) 

instead of ~ ln(— s). Thus the function il(s) is rendered sufficiently convergent to satisfy 
an unsubtracted dispersion relation of the form 
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il(s) 



1 f°° Im[il(t)] 
t - s 



71 J Am 2 



(A27) 



It is instructive to verify this result explicitly, for the calculation shows how the PV- 
regularization of the imaginary part conspires to allow the dispersion relation to hold. One 
finds, using Eq.( |A24Q , that 



il(s) 



1 



(4tt) 



2 /"A* 



dt 



4M? t 



l- A -^6{t-AMl). 



(A28) 



The auxiliary cutoff \i? has been introduced as a calculational device since the individual 
integrals diverge before carrying out the PV summation. Thse can be evaluated explicitly 
using the change of variable t = 4M 2 cosh 2 (j). Assuming s < for definiteness, one obtains 



il(s) 



1 



(4< 



V 



lim ^C a \2J&coth- l Jf a (l 



a=0 



4Ae 



1 — 2 cosh 



-i 



\4M 2 



(A29) 



Taking the limit /i 2 — > oo is harmless in the argument of the inverse hyperbolic cotangent 
function, but the second term diverges logarithmically. The PV-instruction removes this 
divergence and replaces it by a finite constant in the standard way as follows, 



2 lim ^2 Ca cosh 1 



a=0 



\4M 2 



^oc a=0 Mi 

2 u 2 2 M 2 

- lim £CUn^ + £C< a ln a 



a=0 



a=0 



(A30) 



since, according to the PV regulation rules, the instructions of Eqs. ( |A2|) take precedence 
over taking the limit fi 2 — ► oo. Inserting this result back into Eq. ( |A29| ) and reinstating an 
irrelevant —2 under the PV-regulation bracket, we regain Eq. (|A14|) . 

On the other hand, if the PV instruction on the imaginary part given by Eq. (|A24|) 
is dropped, then il(s) diverges logarithmically as we have seen. Then, not il(s), but only 
il(s)/s, is sufficiently convergent to give no contribution upon integration around a closed 
contour at infinity in the complex s plane. Hence the dispersion relation for this repre- 
sentation of il(s) requires one subtraction, 



il(s) = i/(0) + 



00 Im[il(t)} 

(Jjv~ 



71 J4m 2 



t{t - s) 



(A31) 



where only the (formally) infinite constant il(0) now requires regularization. This "disper- 
sive" type of regularization prescription was followed in Ref. [^3| for example, in conjunction 
with a proper time cutoff. If we opt for PV-regulating il(0), then Eq. ( |A31|) leads back to 
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Eq (|A"8|), or equivalently to Eq. ( [A 14]) again, but with the PV- instruction on the second term 
in either of these equations removed, if the unregulated imaginary part of H{s) as given by 
Eq. (|A24j) is used. One then has 



s f°° , Im\il(t)} s r°° dt \ I - , , , , ,. 

' dt 1 \ = — T~j — \n / — -^7 h = 7T^2 da ln[l -a(l - a)} (A32) 



4m 2 
t 

71 hm? ~" t(t - S) (An) 2 J Am? t (t - S) (Alx) 2 Jo "^"^ m 2- 

under the successive transformation of variables u = (1— 4m 2 /t) 1 / 2 , followed by a = (l—u)/2 
after an integration by parts. Hence 

il(s) = z/(0) + — / X rf « ln t 1 - AaC 1 - «)] ( A33 ) 
(47r)^ Jo m A 

that reproduces Eq. (|A8|) . 

The origin of the difference between Eq. ( |A8| ) and the above expression for il(s) thus 
lies entirely in the way the imaginary part is treated. This in turn means that either an 
unsubtracted, or a once-subtracted dispersion relation is satisfied by the regulated function. 
Not surprisingly, these two cases can lead to regulated functions that behave very differently 
at large momentum transfers. A summary of the asymptotic behavior of functions that 
involve amplitudes regulated via these two methods is given in Table [VT] . 

3. Pion weak decay constant, form factor and pseudoscalar polarization loop in the 

NJL model 

It is now a simple matter to write down expressions for the regulated pion weak decay 
constant, the form factor, and thus the one-loop pseudoscalar polarization function Il pp (g 2 ) 
in terms of il(q 2 ). One has 

fl = -AN c mHm = -^Y,Caln^ (A34) 



and 



»/Y«A N m 2 

Fp{s) = W) = [1 ~ 4^ {(2J(s/4m2) " 2)}pv] (A35a) 
N c m 2 f-^„ „ „ M 2 \ 1 N c m 2 A 4 . . 



2vr 2 / 2 V a a rn 2 J s n 2 f 2 s 2 

(A35b) 

With the aid of the relation (A~4) between the quark mass and the condensate density that 



follows from the PV-regulated version of the gap equation, the asymptotic behavior of Fp(s) 
can be recast in the transparent form 
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f2i?n m(i/jip) N c m 2 A 4 

f p F P (s) + — = -ln(-s), s ->oo. (A36) 

For time-like s > 4m 2 , the imaginary part of Fp(s) is obtained from Eq. ( |A24j ) : 



ImF P {s) = ^Jiif- ^fo(s - 4m 2 )} py . (A37) 

The dispersion relation Eq. (|A27|) for il(s) immediately translates into a dispersion relation 
for the form factor: 

p, w = - r (A38) 

Note that using the PV- regulated form ImF P (t) automatically guarantees the condition 
Fp(0) = 1, since we have explicitly demonstrated that the dispersion relation satisfied by 
il(s) correctly reproduces il{0) according to Eq ( |A29| ) . 

The pion electromagnetic radius is determined by the derivative F' P (s) at s — 0. This 
derivative is most simply obtained by evaluating 

using the same change of variable t = 4:M 2 cosh 2 <fi as before. No auxiliary cutoff is necessary 
in this case since each integral converges separately due to the extra t in the denominator 
of the integrand. For N c = 3, the mean field approximation to the pion radius is thus found 
to be given by 

<#"=i^? c -ir5%- (A40) 

that essentially coincides with the result of Tarrach given on the right of this equation. 

Finally we can write down the explicit form of the pseudoscalar polarization and its 
asymptotic behavior using Eq. ( |A35b|) . We take the non-chiral case given in Eq. (|2.90 ): 

\PP(„\ (m-m) / 2 (s) 



F ^ =\L - l^ c 4 - (A39) 



2dm m 2 

ln(-s), s ->cx3. (A41) 

TX l S 

Notice that the constant term of the combination (s/ 2 (s)/m 2 ) — > sf 2 Fp(s) — > (ipip)/m = 
— (m — m)/2G\m just cancels the constant term in I\. pp to leading order according to 
Eq. (|A36|) , so that one has to include the next (logarithmic) order in Fp(s) as determined 
by Eq. ( A26 ) to obtain the stated result. The asymptotic behavior of the regulated pseu- 
doscalar polarizations in the chiral limit for either the NJL model [Eq. ( |2.30a| )] or the ENJL 
model [Eq ( p.4S| )1 are also given by Eq. ( |A41| ). 
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4. Evaluation and properties of the vector and axial vector polarization loops 



a. The vector polarization Uj> (q ) 

We set n = v in Eq. ( |2.28b| ) and then contract on \x in order to exploit the fact that the 
vector polarization is expected to be purely transverse. Then 

= 311^6 = 2iN c 5 ab J -^i*r[fS(p + q)l,S(p)}, (A42) 
after doing the color and flavor traces. Using the fact that 

tr[^(^+ 4 + m)^^ + m)} = 16m 2 - 8p ■ (p + q), (A43) 

one has 

16 f d A p 2m 2 — p ■ (p + q) 



(27r) 4 \{p + q 2 — m 2 ][p 2 — m 2 ] 

(An) 2 Jo Jo p 

32i N c f 1 , r°° , d e-vK-? 2 ^ 1 -")] 



, , da dp—- . (AAA) 

3 (An) 2 Jo Jo dp p 

The p integrals diverge like hip and 1/p respectively at the origin. Under PV-regularization, 
the first integral is just I\(q 2 ) as given by Eq. ( |A7[ ), while the second integral makes no 
contribution at all since 

2 i-oo A P -ip\M^-q 2 a(l-a)] 2 1 

lim n ECW dp — =-\imJ2C a [--z[M 2 -q 2 a(l-a)}}^0. (A45) 

v^o^ Jr, dpp ^°^o V 

Thus, with s = q 2 again, 

N r 1 

U^ v (s) = -16s— 4- / daa(l-a)h(s) (A46a) 
(An 2 ) Jo 

N r s^„ , M 2 A 



Qn 2 q m 3 

^£VC a ln — = --S 

6n 2 a a m 2 3 m? 
1 



2 [£ C a In -§■ + {- + (3 - f)(y/f coth- 1 y/f - 1)} PV ] (A46b) 



^C a ln— = --^ S , s^O (A46c) 
„ m 2 3 

ln(-s) , \s\ -> oo. (A46d) 



with the indicated limiting behavior near the origin and on a great circle in the complex 
plane. Eq. ( |A34T ) has been used to re-introduce the mean field pion decay constant. In 
obtaining Eq. ( |A46b|) , we have used the intermediate result that 
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f daa{\ - a) ln[l - -^a(l -a)] = \[\ + $- /)( V / 7coth- 1 yff - 1)] (A47) 

J Tft O O 

for space-like s, but which may immediately be extended to cover all s, real or complex, 
using the definition of the function J(z) in Eq. ( |A11| ). 

Thus like il(s), only the PV- regulated vector polarization is sufficiently convergent on 
the great circle to satisfy an unsubtracted dispersion relation, 

in- W = -L r * Im[W(t)] r * — < A48 > 

4 47T 7 4 m 2 t - S J4m 2 t - S 

where 

pf )V (s) = -^Im[t^ v (s)} (A49) 

is the vector spectral density function. Its explicit form is found with the help of Eq. (|A20|) 
to be 



(o)v, , N c s r I 4m 2 ,, 14m 2 . 2 .-. 
Pi 00 = 1 - — (1 + 3—)^ _ 4m ) W ( A5 °) 



Again we comment that this expression is no longer positive definite for all s for the same 
reason as before, due to the PV-instruction. One has a linearly divergent behavior p\{s) ~ s 
as s increases if the PV-instruction is removed. If it is kept, however, a short calculation 
shows that 

Pi ( s ) ' s ^°° ; ( A51 ) 

Z7T Z S 

which is convergent at the expense of becoming negative at large s. However, such a behavior 
seems unavoidable if at the same time one wants the vector polarization as given by the NJL 
model to obey an unsubtracted dispersion relation. Indeed, if the PV-instruction is removed 
in Eq. ( A46b| ) (and therefore also in Eq. ( A50 ), rendering pf^ v (s) positive definite) then 
Wf v (s) diverges like s ln(— s), so that two subtractions are required in the dispersion relation 
for this Uj- V (s) (or equivalently, a once-subtracted one for the combination Il^ y (s)/s - see, 
for example, Ref. fll3|), 

n-M = s n? l "' ( o) + tr ( A52) 



7T Jim 2 t 2 (t 

where the prime denotes the value of the derivative at the origin. As with the case of il(s) 
and il(s), one can again demonstrate explicitly by direct integration using either the PV- 
regulated, or the unregulated pf^ V \s), that the dispersion relations ( |A48 ) or (|A52j ) lead back 
to the expression for Il^ y (s) in Eq. ( |A46 b| ) with the PV-instruction either present or absent 
on the second factor. 
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b. The longitudinal and transverse axial vector polarization H AA (q 2 ) and Tl AA (q 2 ) 

We proved in Eq. ( j2.47| ) that the longitudinal axial vector polarization must vanish when 
all interactions are considered. At one-loop level, however, it does not vanish and is given 
by Eq. (P7J), i.e. 



n£ V) = iftFpW) (A53) 



after regularization, where Fp is taken from Eq. ( |A35b| ) and f 2 from Eq. ( |A34[ ) . We use 
Eq. ( |2.37| ) to determine the regulated transverse axial polarization as 

fl AA (s) = U^ v (s) + Af p F P {s) (A54a) 
= 4/ p 2 , s^O (A54b) 

~ — ln(— s) , s —>■ oo. (A54c) 

Again Yl^, A (s) will satisfy an unsubtracted dispersion relation like Eq. (|A48|) containing the 
axial vector spectral density that is related to (s). This can be seen by taking the 
imaginary part of Eq. ( A54a ) and using Eq. ( A.37 ): 



pr<s) - pr<s) = f ImlFAs)] = - ^0(s - 



N c m 2 A 



2 \ ! 

, s — > oo. (A55) 



7T 2 S 2 



Hence the density difference that enters into the two Weinberg sum rules converges like 
1/s 2 , i.e. faster than either density separately and renders both sum rules convergent for the 
PV-regulated spectral densities. By contrast, with the PV instruction removed, pi V (s) — 
p^ A {s) ~ N c m 2 /(A7c 2 ) asymptotically, and neither sum rule would converge. The regulated 
axial density itself is given by 



pTM = ^ ± C.Jl - %1 - ^) + 6 -(M! - ra'Ms - 4Mf )}. (A56) 

/47T a=Q V S S S 

The first Weinberg sum rule is explicitly satisfied by the spectral densities of the minimum 
NJL model. We have 

r - u o)v w - pr w) = £ r ^/m[F P( t)] 

JO S 71 Jim 2 t 

= f p FM = fl (A57) 



after using the dispersion relation, Eq. ( |A38j ), at s = to deduce the value of the integral. 
To check the second sum rule we require 
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, ds („f >"M - pf" W) = ^ »m £ C. I'* dsJl-'-Ml 6(s - AMD 

,2Af S/f2 . ,,2 



or 



— ^ im Y, c *\l* Y 2 ~ M 2 a cos^ 1 {Z- s 

'"^ECo^ln^fi, /x 2 -oo (A58) 



2tt 2 V a m 2 



ds (p { i )V (s) - pf )A (s)) = -m(^>, (A59) 
where the last step follows from Eq. 



o 



5. Electromagnetic form factor of the pion in the ENJL model 

In the absence of vector and axial vector meson fields, the electromagnetic form factor 
of the pion is given by the Fp(q 2 ) of Eq. ( 2.32c ) to O(N ) in the mean field approximation. 



In the presence of the vector fields, this is no longer the case, since the axial form factor 
GU(0) = <7a < 1 differs from unity. This renormalizes the irqq coupling constant, g p — > 
9-nqq — 9p/ \/~9a-> that in turn renormalizes the pion charge so that -Fp(O) ^ 1. However, 
as we now show, a compensatory renormalization of the pion field also occurs due to the 
pseudovector coupling and pseudoscalar-pseudovector interference terms introduced by the 
axial vector meson field that restores the total charge of the pion to its physical value. The 
inclusion of the vector meson fields does not alter this conclusion since -FV(O) = 1. 

The set of diagrams that determines the charged pion's electromagnetic form factor in 
the ENJL model is obtained by coupling the photon to both the quark and antiquark lines 
making up the 'jirn triangle diagram. Since the electromagnetic current is renormalized as 
per Eq. ( |4.1|) , the photon vertex contains a contribution from the p° exchange diagram too. 



The resulting jim vertex thus has two contributions as depicted in Fig. [T5|. We work in the 
chiral limit of zero pion mass. Then the contribution from the first diagram of Fig. [l^ to 
7rqq vertex can be written as a four vector of the form 

\T^ m \k', k) = ( ig7Tqq ) 2 K(0)-y^ m \k', k). (A60) 

Adding in the second diagram replaces T^ m \k\ k) by 

r ^)(k', fc)- T^ m \k\ k) - U v J(q 2 )D vv ^(q 2 )Ti em \k', k) 
= [L lw + F v (q 2 )T^]T^(k',k) 

= F v (q 2 )T ( ~r\k',k). (A61) 
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The last expression is only valid for on-shell pions. The electromagnetic form factor F n (q 2 ) 
of the pion can then be identified from 

e^{q 2 ){k + fc% = T^ m \k', k)F v (q 2 ) (A62) 

We thus only need to evaluate the contribution from the first diagram. Its spinor structure 
can be read off from Eq. ( |2.67| ) as P a cos 9 ± i(k • A a ) sin 9 for outgoing (incoming) pions of 
isospin a using the abbreviations P a = i^j^T a and A a ^ = 7^751"°. We first construct the fully 
off-shell version of —iV^ em \k',k) and then specialize to on-shell chiral pions to obtain the 
charge form factor. One finds 

W^ em) (k',k) = cos 9' cos 9^V pp (k', k) + cos 9 sin 9' V A » p (k',k)k w 

- cos 9' sin 9 V PAv (k', k)k v + sin 9' sin 9 k' u -V A " A "{k' , k)k p (A63) 

for absorbing a photon of four momentum q = k' — k. The individual amplitudes in this 
expression arise from the various combinations of pseudoscalar and axial vector vertices that 
are introduced through the revised nqq interaction. They are given by 
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V pp (k', k) = (-2ie„N c ) J P t r[ l5 S(p + k')j^S(p + k) l5 S(p)) (A64a) 



(2*0 

V PA ^(k' } k) = {-2e„N c ) J -0-tr[ l5 S(p + k') lfl S(p + k) lul5 S(p)} (A64b) 



V AuAp (k', k) = (2ie^N c ) / — -t-tr[ lu75 S(p + k')^S(p + k) lpl5 S(p)]. (A64c) 



d A p 

" v '"'" v J {2m)A- 

The symbol = eT 3 is the charge carried by the pion. 
A direct calculation gives jT5[] 

p 



V pp (k\ k) = -e n ^-F P (q 2 ) (k + fc% , onshell, (A65) 



for the first amplitude, a result that also extends to the 1/2 off-shell amplitude in the chiral 
limit. On the other hand, repeated use of the Ward identity 

^7 5 = 2m 75 + S' 1 (p + % 5 + 75<S'- 1 (p) = - [2m 75 + l5 S-\p + k) + S^phs] ( A66) 

leads to the results 

V PA »{k', k)k v = -2imV pp (k', k) + e n U PA {k' 2 ) 

= -2imV pp (k', k) - ^f 2 p (k') k^, (A67) 



SO 



using Eq. ( g305[) for U PA , and 



k'»V^ A >(k', k) = 2imV p \k', k) + e.[n^(fc 2 ) - H^g 2 )]. (A68) 
It follows that 

W^tf, k)k" = 2imV PA »{k\ k)k» + e n [U AA (k 2 ) - (A69) 

Only the longitudinal part of the axial polarization contributes to this expression, since 
H AA (k 2 )k v — Ii AA {k 2 )k^ = Af 2 {k 2 )k^. However, there is also a contribution from the (purely 
transverse) vector polarization, since its argument involves the momentum transfer g M = 
{k'—k)n that is different from k^. Then T^ u (q 2 )k v = T^ u (q 2 )k' v ', since T^(q 2 ) = g^-q^/q 2 , 
and 

nJZV)*" = nr( g 2 )T p ( 9 2 )r = n^VX? 2 )** = ^n£V)W)(* + k'y. (A?o) 

Hence 

WfA'tf, k)kf = 4m 2 V pp (k', k) + Ae^[f p {k' 2 )k^ + / 2 (A; 2 )^] 

+ ^X l/ (9 2 )^(9 2 )(H^T. (A71) 

Collecting terms, one finds that the final form of the off-shell electromagnetic vertex is 

V( em) (k',k) = (cos 6' - ^Lsine')V pp (k',k)(cose - -^Lsintf) 

V k' 2 v A; 2 



_ 2 ^-f 2 (k' 2 )k f ll ( C os6'- ^=sm6" Sin6 



Ctt .9/, o, , , „ 2m „. sin^' 
- 2— f 2 (k)k ^(cos 6 -j= sin 9)—r= 
m v Jk 2 vk' 2 



- \e^V{q 2 )T^ q 2 ){k + k'Y ^^=. (A72) 

after using the symmetry property V AuP (k', k) = —V PA "(k,k'). The on-shell value for 
chiral pions, k' 2 = k 2 = m 2 = 0, is obtained by inserting the limiting values cos# — > 1 and 
(2m/\/~k?) sin# — > 8G2/ 2 — (1 — 9a) which come from Eqs. ( |2.71| ) and hold for both the 
primed and unprimed quantities. Then 

Vjr\K, k) = -e^{g A F P (q 2 ) + (1 - g A ) + (1 ^^ ^ v (q 2 )}(k + k% (A73) 



is manifestly gauge invariant. In conjunction with Eqs. (A60) and (A62), it leads to the 
identification 
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FM 2 ) = {9AF P (q 2 ) + (1 - g A ) + (1 - g A f^U^ v (q 2 )}F v (q 2 ), (ATA) 

after invoking the Goldberger-Treiman relation f^g^qq = m once more. Thus the electro- 
magnetic form factor of the pion has the required property, F v (0) = 1, that avoids charge 
renor malizat ion . 



6. Evaluation of the ENJL Feynman diagrams that enter into AraJ 



In this section, we obtain closed forms for the Feynman diagrams shown in Fig. [20 
that enter into Am^. We start with the scattering diagram in Fig. pp|(a). Translating the 
diagram, one finds in the chiral limit for the external pion that 

i{ n -(°)L it = ^ I j^M T h T ^ 



(2tt)V 

+tr(T / t^r l |)(nr( ? 2 ) J Dr(g 2 ))V / 



2\ 



(A75) 



with 



I" u (q 2 ) = J -0-,tr[ l5 S(p + lh,S(ph 5 S(p) lu S(p + q)]. (A76) 

As before, e q = |(| + T3) is the quark electric charge in units of e, and Tj and Tj are 
isospin creation and destruction operators for the initial and final pion channels. The photon 
propagator is in the Feynman gauge of Eq. ( |4.16| ). In obtaining these results, we have used 
the fact that the effective ^qq vertex, including p° exchange, is given by Eq. ( |4.1[ ), i.e. 

- ^ em V) = -ie[e q l, + ^ v (q^D™ (q 2 )T^r 3 Y)}- (A77) 

The values of the isospin traces in II^ 3 M of course depend on the charge channel T3 in 
question. They are listed in Table |V11|. We are, however, only interested in the difference 



II em - U° EM . With the help of the isospin table, one finds the surprisingly simple result for 
this difference, 



i{n 



d 4 q 



™ ~ ^"L* = (^hi L ^ + T^ v tf)}l^{q% (A78) 



(2vr) 4 g 

after using g^ = L^ u + T M „ and inserting the explicit form of the transverse p° meson mode 
propagator from Eq. fl2.52|) , plus the definition of the vector form factor from Eq. ( (2.88|) . 
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Finally we compute I^ u (q 2 ). From its spinor structure, I^ v (q 2 ) has to be a linear combination 
of the longitudinal and transverse projections with coefficients that are only functions of q 2 . 
Upon evaluating the integral, one in fact finds that the two coefficients are equal, and 

I^(q 2 ) = -AI{q 2 )g» v (A79) 
where I(q 2 ) is defined by Eq. ( |A5|) . Consequently, 

i{nt„ - ntu}_ = p^/tf) + v*?».f jSk Iiq2)FHq2) - 

We eliminate I(q 2 ) in favor of the pion form factor Fp(q 2 ) and pion decay constant f p , 
Eqs. ( p.32c| ) and ( |2.32a| ), that refer to the mean field calculation, to obtain the result quoted 
in the main text as Eq. fl4.17|) , 

(2tt)V 

-3e 2 4 I 7^F P (q 2 )F 2 ( q 2 ) (A80) 



m 2 J (2tt)V 

The next task is to evaluate the dumbbell diagrams of Fig. pO|(b). Due to the presence 
of the electromagnetic vertices they contain, there is only a contribution in the charged 
channels from these diagrams. Let us start with the diagram where a pion is exchanged in 
the intermediate state. Label this diagram n^2r- Moving immediately to the chiral limit 
for the external pions, one obtains 

\^eIM = J j0^V^ m \O,q)XAq 2 )9" U Dpp(q 2 )V^ m \q,O). (A81) 

The 77T7T electromagnetic vertex is given by Eq. (|A72 ) , and the pion eigenmode propagator 



can be identified from Eq. ( |2.67 ). This expression does not yet contain the contributions from 



the three remaining diagrams of Fig. |2rj(b) that involve the p° meson to photon conversion 
in the intermediate state. However, none of these diagrams contribute in the chiral limit as 
we now show. From Eq. ( A72j ), the half off-shell behavior of the electromagnetic vertex is 



f 2 ?m 
Vf m \O jq ) = Vf m \q,0) = -e n ^F P (q 2 )(cos6 - s fn%, (A82) 

in the chiral limit, which is purely longitudinal. Consequently the additional p° meson 
diagrams, which are automatically included together with the direct photon exchange by 
writing 

l*«")( , q) - [Lfjj, + F v {q 2 )T^]V^ v {^ q), (A83) 



S3 



as in Eq. ( |A61| ), fail to contribute in this case due to the longitudinal nature of V^ em )(0, q); 
the transverse piece has been projected out completely. The corresponding diagram for the 
a>L exchange has exactly the same form and properties. We simply must substitute the 
eigenvalue \ aL (q 2 ) for A^-(g 2 ) in Eq. ( |A81| ). Then the n and exchanges taken together 
contribute an amount 



ing? £) (0) = e 2 ( Jr f J A.(g 2 ) [(cos 2 9 + ^ sin 2 9)-^ sin 9 cos 9] 



+\ aL (q 2 ) [(sin 2 9 + ^ cos 2 9) + ^= sin 9 cos 9] } D pp (q 2 )F 2 (q 
q v'/ 



(A84) 



the factor 1/g 2 from the photon propagation having cancelled against a similar factor coming 
from the contraction of the two electromagnetic vertices. The combinations appearing in 
the integrand have the values 



A,,- cos 2 9 + \ aL sin 2 9 



+ A,-A a , cos2 9 = A 



(\ w — X aL ) sin 9 cos 9 



K XaL sin29 = B 



\ • 2 n , \ 2 n "I" \ aL \ aL 

X n sin 9 + X ar cos 6 1 = — - cos 29 = C 



(A85) 



from the properties of the eigenvalues and eigenvectors given by Eqs. ( |2.64|) and (|2.65| ). The 
value of the contents of the curly brackets in the integrand of Eq. (|A84|) is thus 

{ ' ' ■ } = W ^ - ^B) = 1 - SG^tf) = TT ^ ?? (A86) 

after inserting the values of A, B and C from Eqs. (|2.59| ) to (|2.60| ) and using the definition 
of f 2 {q 2 ) from Eq. ( 2.41|) . Since the renormalized pion propagator of Eq. ( p. 61 ) contains the 
inverse of the same factor, 

1 



D PP (q 2 )=g;[l + 8G 2 f^q 2 )} 



(A87) 



, q 2 F P (q 2 ) 

this factor cancels out in the final expression for the integrand, together with one of the 
Fp(g 2 )'s. Hence we obtain the simple expression 



■yr(Tr + a L ) 

' LL EM 



(0) 



f 2 

2 Jp 



d 4 q 



m 2 J (2tt)V 



:F P {q 2 ) 



(A88) 



as the final result after using the Goldberger-Treiman relation f 2 g 2 = m 2 for the mean field 



quantities once again. The final answer in fact coincides with the NJL result where no vector 



mesons are present |T5 
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The calculations of the dumbbell diagrams involving the a% and p mesons follow an iden- 
tical pattern. The basic electromagnetic vertex for these two cases is obtained by replacing 
one of the incoming (or outgoing) 7r's in Fig. [TBI by an ai or a p respectively. In the first 
case this leads to 

V^\k\ k) = V? A »{k>, k){cos6' - sine') + ie^(k 2 ) - n^(g 2 )]^| (A89) 



for an incoming ai(k) and outgoing ir(k'); here g M = (k f — k)^ is the momentum transferred 
by the photon. The required half off-shell value in the chiral limit is 



V^\0, q) = V^(0, q)[l - 8G 2 f 2 ] + i^8G 2 f 2 [U AA (q 2 ) - U v J(q% (A90) 



due to the limiting values cos 6' — ► 1 and 2m sin 9' / yk/ 2 — > 8G2/ 2 for the emitted chiral 
pion. But the axial vector minus vector polarization equals 4f 2 Fp(q 2 ) according to the 
Ward identity ( |2.37|) , and 

V, PA »(0,q) = 2ime^F P (q 2 )g^, (A91) 

that follows immediately from the relation (|A67 ) and the expression ( |A65| ) for V pp (k', k). 
The terms involving G 2 cancel, and one is left with V^ ai \Q,q) = V PA -(0,q). The addi- 
tional three diagrams that include the p° degrees of freedom can be incorporated into this 
expression by replacing g^ u with [L^ + Fi/(g 2 )T M ^] as before. The final expression for the 
777a! vertex including p° contributions is thus 

V^(0, q) = 2ime^F P {q 2 )[L, v + F v {q 2 )T, u ] = -V^(q, 0). (A92) 



The structure of the a\ dumbbell diagram is thus 

Tng&(0) = / (^^ ai) (°' q)9^'D£ A (q 2 )T™'V^\q, 0). (A93) 

In contrast with the case of it + exchange in Eq. (|A84 ) where the transverse part of the 
vertex involving Fy{q 2 ) fell away, here only the transverse part of V^2 nai \0, q) survives due 
to the transverse nature of the a\ mode propagator; now the longitudinal part falls away 
completely. One has 

[V + F v (q 2 )T^}T uv '{L^ u , + F v (q 2 )T^ u ,} = F 2 {q 2 )T^. 

Using g^'Tppi = 3, the electromagnetic polarization contribution due to the a\ exchange 
becomes 
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^(O) = 12eW(-|) 2 / -^i% 2 )W)^V), (A94) 



(27T)4g 



after inserting Eq. ( |A92| ) 



The basic •ynp vertex is obtained by replacing the spinor A a v = ^v^t 0, by V u = •y u r a 
in Eq. (|A64b| ) . Then this diagram assumes the same structure as the anomalous it — > 27 



vertex, see for example Ref. ||18|| , p 552. After attaching the photon to both the quark and 



antiquark lines, and adding, one finds 

V™>(V,k) = (-eN c )tr({T},T t }e q ) J ^-tr[ lh S{p + k') lfM S(p + k)j„S(p)} 

= (~eN c )(Aims^ pa k' p k°) x J(k\ k), (A95) 

since the trace that arises in the numerator of the integrand after inserting the quark prop- 
agators has the value Aime llV p ( jk' p k a ' . The factor J(k', k) is the integral 

J ^ k ' ^ = / (2vr) 4 [(p + k') 2 - m 2 } [(p + k) 2 - m 2 } [p 2 - m 2 ] ' 

The three additional diagrams containing the p° can be incorporated by including the factor 
[L^i + Fv(q 2 )T^r] as before, and we obtain the final result 

V^o\k', k) = ~ieN c mJ{k\ k)[L^ + F v {q 2 )T^' vpa W. (A97) 

The half off-shell value of this vertex is zero in the chiral limit, V$ np) (Q,q) = 0, so the 
dumbbell diagram with a p replacing the a\ vanishes. 



7. Dashen's theorem 

The complete set of diagrams that determine the corrections to the pion polarization 



to 0(aN c ) are shown in Figs. [19] and |20j. We have already calculated the scattering and 



meson pole, or dumbbell contributions of Fig. EUJ in the previous section. The remaining 



diagrams in Fig. [19] dress either the quark or antiquark line, i.e. they change the quark mass 
circulating in the polarization loop. Since the quark and the antiquark lines are dressed 
identically by the same interaction to 0(a), both diagrams have the same structure apart 



from their isospin weights. One finds for Fig. [TpJ (a) + (b) that 

r4T(0) = -e 2 N c J J^{2tr(e 2 q )9» u + (W) - 1)T^}m^ (A98) 



where 
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M, 



d A p 
12k) 4 



tr 



But 



gTM^ = 4[7(g 2 ) + 7(0) - {q 2 + 2m 2 )K(q 2 )] and L^M„ V = 2I(q 2 



where we have defined 

K(q 2 ) 

Hence 



d 4 p 



(27r) 4 [(p + q) 2 — m 2 } [p 2 — m 



212 ' 



(A99) 



(A100) 



(A101) 



n 



(a+b) 
EM 



(0) = e 2 N c 



d A q 



(27r)V {8tr(eJ)[-7(g 2 ) - 7(0) + (q 2 + 2m 2 )K(q 2 )]. 
-(F 2 (q 2 ) - l)[27(g 2 ) + 47(0) - 4(g 2 + 2m 2 )K{q 2 )]}. 



(A102) 



Diagrams 18(c) and 18(d) involve the dressing of the quark and antiquark lines via the 
exchange of a a mode with the condensate. The sum of these diagrams can be expressed in 
terms of the quark electromagnetic self-energy T>em, since the a propagator terminates on 
an electromagnetic self-energy insertion of the full quark self-energy: 



n (c+d) 
11 EM 



(0) 



J EM 

2m 



(A103) 



One has 



Zem = (^GiJVc) / (^^{[Me 2 ) + (F 2 (q 2 ) - l)]7/(g 2 )} 

= {%e 2 mG 1 N c ) J ^f^{[2tr(ej) + (F 2 (g 2 ) - 1)] 
x[7(g 2 )-7(0) + (g 2 + 2m 2 )^(g 2 )]} 



(A104) 



using 



d A p 



trh,S(p + q)rS(p)S(p)} 



(2vr) 4 ~' l,m 
8m[7(g 2 ) - 7(0) + (q 2 + 2m 2 )7<V)] 



(A105) 



The total contribution to the pion polarization from Fig. [19| is then given by the sum 

ngj?(0) + ngj?(0)] = -e 2 N c J j0^{Wtr(e 2 q ) + 6(F 2 (g 2 ) - l)}7(g 2 ) (A106) 

On the other hand the polarization contribution due to the interaction of like-flavored qq 
pairs in the T 3 = channel can be found from Eq. (|A75 ). It is 



^7 



{n EM(o)} sca u = **cf j0^{ mr «) + 6 (W) - i)}i(q 2 )- (Mo?) 

This is precisely the negative of the sum in Eq. ( |A106[ ). If we convert these two expressions 
into pion self-energies by supplying the missing coupling constant (ig n qq) 2 and using F P (q 2 ) = 
I(q 2 )/I(0) plus the GT relation f 2 qq gl qq = m 2 of the ENJL model, one reaches the final 
conclusion that 

E self = -Y? scatt = ie 2 g A l J ^S_{ur{e 2 q ) + \(F 2 (q 2 ) - l)}F P (q 2 ). (A108) 

This shows the exact cancellation between the self-energy of the neutral pion that arises from 
the electromagnetic self-energies of the quarks it contains, and the self-energy due to their 
electromagnetic interaction, and illustrates how Dashen's theorem is realized in the ENJL 



model. A direct comparison of the various self-energies is given in Fig. |2Tj. One sees there 
that £ se j/ > 0, and that the scattering contribution is the negative of this, in agreement 
with the attractive electromagnetic interaction between like-flavored qq pairs. 
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FIGURES 

FIG. 1. Modification of the quark axial isospin current vertex due to the exchange of the pion 
plus longitudinal axial eigenmode, and the transverse a\ mode. Quark propagation is shown as 
solid lines, mesons as broken lines. The wavy line represents an isospin current component of 
momentum q entering the vertex. 

FIG. 2. Self-energy diagrams entering the non-chiral gap equation to 0{rh) for m* . (a) 
The Hartree contribution — iTiH in Eq. Q2.34 ). In this figure, (b)+(c) represent the non-chiral 



corrections arising from a finite current quark mass, — irh, denoted by the cross. The a mode 
exchange is denoted by the double solid line. 

FIG. 3. Behavior of the quark mass and condensate density as a function of the ratio 
x = A 2 /m 2 . The upper pair of curves labelled m and — (qq) 1 ' 3 refer to the ENJL calculation 
with f n = 93MeV and gA = 0.75. The lower pair of curves refer to the results for the NJL model. 
They were obtained by scaling the upper pair by yfgX = 0.866. 

FIG. 4. The vector and axial vector spectral densities of the chiral ENJL model versus the 
momentum transfer squared s = q 2 . The interaction strengths and regulating cutoff are those given 
in the second line of Table | The qq threshold then lies at Am 2 = (2 x 0.264 GeV) 2 = 0.279 GeV 2 . 
The peak values of the vector and axial strengths lie at m 2 = (0.713 GeV) 2 and m 2 x = (1.027 GeV) 2 
respectively that are taken to identify the p and a\ masses of the ENJL model. Their ratio 
m a 1 / m p = V2.07 lies close to the original Weinberg estimate || of \pl. 

FIG. 5. Comparison of the PV-regulated density functions of the ENJL and NJL model respec- 
tively. The parameters are the same as for Fig. ||. The position and spacings of the qq threshold 
at 4m 2 = 0.279 GeV 2 , and the two thresholds introduced by the Pauli-Villars regularization at 
4m 2 (l + A 2 /m 2 ) = 4.8 GeV 2 and 4m 2 (l + 2A 2 /m 2 ) = 9.2 GeV 2 for this parameter choice are also 
shown. Note (a) the strong redistribution of strength to lower energies of the vector mode, and (b) 
the strong suppression of the unphysical negative density at large s for the ENJL case. 
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FIG. 6. Comparison of the PV-regulated ENJL and NJL strength function differences p\ — p± 
and p^ V — p^ A as obtained from Eqs. ( |3.15|) and ( 3.16 ) and Eqs. ( A50f) and ( |A55 ). The pa- 



rameters are the same as for Fig. |j. The PV induced thresholds are as described in the previous 
figure. The strong concentration of strength at much lower energies in the ENJL case is due to the 
presence of the vector and axial mode form factors that essentially damp out the PV-induced un- 
physical behavior of the NJL spectral density at energies considerably below the first PV-generated 
threshold. 

FIG. 7. The behavior of the real part of the function J{z) along the real axis on the first 
Riemann sheet, and its analytic continuation J(z) through the cut z = 1 to oo, onto the second 
sheet. 

FIG. 8. The behavior of the real part of the vector polarization Y[)f v (s) as a function of the 
variable z = s/4m 2 on the first and second sheets of the cut s-plane with the cut running from 
z = 1 to oo along the real axis. The values of the polarization along the upper lip of the cut on 
the first sheet join smoothly with their analytic continuation onto the lower lip of the cut on the 
second sheet. At the branch point these values bifurcate into two branches of the function that 
assume different values along the sector < z < 1 of the two sheets. As explained in more detail 



in connection with Eq. (A21), the values on the first and second sheets are reached by allowing the 
variable z to pass either infinitesimally above or through the cut infinitesimally below the branch 
point to reach this sector of the real axis. The solid circles indicate the values of s/4m 2 where 
the polarization function equals — l/(8G2m 2 ) for the indicated interaction strengths G*2. These 
intercepts determine the position of the real poles of vector mode propagator lying in the interval 
< s/Am 2 < 1. 



FIG. 9. The switch-over of the real pole of the p propagator Dj. v of Eq. (2.52) from the second 



to the first Riemann sheet as the interaction strength G2 increases through the values 2 (1) 8, 10 
and 16 GeV -2 (solid circles). The critical value of G2 at which the switch-over from the second to 
the first sheet occurs is G 2 c) = 7.40 GeV" 2 in this illustration. 

FIG. 10. The positions of the complex poles on the second Riemann sheet of the p and a\ 
propagators in Eqs. ( 2.52; ) and ( 2.54 ) in units of 4m 2 . The value of the interaction strength G2 is 



recorded in GeV 2 next to each pole. 
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FIG. 11. Top figure (a): Comparison of the behavior on the first Riemann sheet of ReILj, J (s) for 
J = V, A as obtained from an unsubtracted dispersion relation, plotted as a function of s = q 2 /Am 2 . 
The bottom figure, (b), gives the same information for ReHj> J (s), now obtained from a once 
subtracted dispersion relation. Note the occurrence of a space-like zero in both functions in this 
case. 

FIG. 12. Landau ghost poles in the p, a±, and tt channels. These poles correspond to the 
space-like roots of the combination l + 2GII(s) where the polarization function has been calculated 
via a once-subtracted dispersion relation. For the vector and axial vector modes, G = G2 and II(s) 
comes from Eqs. ( 3.2l| ) or ( |3.23| ) with the bars removed. The roots corresponding to the resulting 



p and a\ ghost modes are almost degenerate at s/4m 2 = —18.30 and —18.04 respectively for the 
standard parameter set of Table |[ For the pion mode, II = ti pp of Eq. ( 2.48j ) and G = —G\. 



Apart from the expected Goldstone mode for the pion at zero, a second, ghost pole appears at 
s ng /Am 2 = —13.48. The extensions of the curves beyond the cusp at s/4m 2 = 1 refer to the real 
part of the corresponding polarization function. 

FIG. 13. Comparison of the modulus squared of the exact vector form factor with the single 
pole approximation. The p pole is located at s p = 4m 2 (0.522 — 1.326 i) = (0.146 — 0.370 i)GeV 2 on 
the second Riemann sheet, as illustrated in detail in Fig. [l5|. Only that part of the approximate 
curve that lies to the right of the threshold at 4m 2 = 0.279GeV 2 , marked by the cusp in the exact 
curve, is accessible on the real axis of the physical sheet. 

FIG. 14. The same comparison as in Fig. [13| between the exact and single pole ap- 
proximation, but for the axial form factor |Gvi(s)| 2 . In this figure, the a\ pole lies at 
Sai = 4m 2 (2.123 - 2.141*) = (0.592 - 0.597i)GeV 2 . 

FIG. 15. Closed integration contour on the first, or physical, Riemann sheet for evaluating the 
sum rules. The physical cut starts at s = Am 2 . Additional poles can appear on the real axis of this 
sheet at bound states of the qq system in < s < Am 2 if G2 is large enough, or as ghost poles in 
— 00 < s < 0. As discussed in the main text, such ghosts only appear if once-subtracted dispersion 
relations are used in the construction of the propagators. The crosses indicate the pole positions of 
the virtual p mode at 0.1658, and complex p and a± modes at 0.5221 — 1.326i and 2.1227 — 2.1405i 
respectively in units of 4m 2 for the standard parameter set of Table |[ As these singularities all lie 
on the second sheet, they do not fall within the integration contour. 
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FIG. 16. The ^tttt electromagnetic vertex for absorbing a photon of four momentum q = k' — k, 
where k and k! are the initial and final pion four momenta. Quark propagators are shown as 
solid lines, pions as broken lines and photons as wavy lines. In this figure, the small open circle 
represents the qq vertex of Fig. [I] with the isospin current replaced by the electromagnetic current. 
Thus the pseudoscalar and pseudovector exchange contributions in that diagram are to be replaced 
by p° exchange, indicated by the heavy solid line. 

FIG. 17. The space- and time-like behavior of the electromagnetic form factor \F 7T (q 2 ))\ 2 of the 
pion. The experimental data (open circles) have been compiled from the various sources listed 
under Ref. [^TJ . The calculated value of this quantity from the EN JL model is given by Eq. ( f4.2| ) . 
The modulus squared of this expression has been plotted as a solid line. Although the ENJL model 
is restricted to low energies as an effective theory, where the agreement is quantitatively significant, 
the plot has been extended to cover the full range of available data. The possibility of p — uj mixing 
as evidenced by the interference structure around q 2 ~ 3GeV 2 has not been considered in the ENJL 
model employed in this paper. 

FIG. 18. The calculated electromagnetic mass squared splitting of the pion as a function of the 
photon cutoff squared, Ap feoton for the standard parameter set of Table |. The observed value of 
the splitting, Am 2 = 1261MeV 2 is indicated by the horizontal line. This value of the splitting is 
obtained for X 2 photon = (1.27MeV) 2 . 

FIG. 19. The complete set of pion polarization diagrams to 0(N C ), where the mass of either the 



quark or the antiquark line is renormalized to 0(a). The conventions are the same as in Fig. 16- 
The double solid line indicates the exchange of the a mode of the model. 

FIG. 20. The two classes of irreducible pion polarization diagrams involving single photon 
exchange that determine the electromagnetic mass splitting of the charged to neutral pions to 
0(a): (a) diagrams that describe the scattering of qq pairs in the intermediate state, and (b) 
meson pole, or "dumbbell" diagrams where the exchange of the tt + and a\ meson eigenmodes 
together with the photon between identical vertices occurs. The conventions are the same as in 
Fig. [16|. There is no contribution from the dumbbell diagram with a p substituted for an a\ in the 
chiral limit. The small circle at a photon vertex indicates that p° exchange has to be included as 



in Fig. 16. As a consequence, both diagrams are actually a sum of four diagrams. In (c) we have 



written out this sum explicitly for the dumbbell diagram. 
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FIG. 21. Breakdown of the U v (l) x S77l(2) x S77ij(2) chiral symmetry of the original ENJL 
Lagrangian into Ua(1) x £/y(l) by the photon gauge field. This the origin of the electromagnetic 
mass of the Goldstone pions of the ENJL model. The physical mechanism that keeps the neutral 
chiral pion massless in agreement with Dashen's theorem is illustrated explicitly, with the mass 
shift 

^seif being exactly cancelled out by the internal electromagnetic interaction E sca # of the 
neutral pion. 
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TABLES 

TABLE I. Comparison of the NJL and ENJL model coupling constants, regulating cutoff, as 
well as the quark properties they give at x = 16. The input parameters are f p = 93MeV for the 
NJL case (first row), and f w = 93MeV together with a form factor g A = 0.75 in the ENJL case 
(second row). We refer to the latter set of parameters as the standard set in the text. 



Gi(GeV" 2 ) G 2 (GeV- 2 ) A(GeV) 


m(MeV) 


-(qq)*(MeV) 


3.29 0.916 
2.47 3.61 1.058 


229 
264 


299 


TABLE II. The combination k + ln(m%/ ' n 2 ) and lj of the chiral perturbation 
constants as given by the ENJL or NJL model. 


theory coupling 


i ENJL 




NJL 


3 {87T 2 ^(2g A - 1) - 4g A (3g A - l)}g A 




87r 2j % - 8 


4 {4vr 2 ^ - 3g A }g\ 




4tt 2 4t - 3 


5 1 + 6^ + 8^5^) 




7 


6 1 + 12^ + 8^^(^1) 




13 


7 2rn 2 ft 

9 'A rn% 




2m 2 '^j 
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TABLE III. Comparison of the Gasser-Leutwyler chiral perturbation theory parameters with 
the predictions of both the NJL and ENJL model calculations. The relevant input parameters are: 
f n = 93MeV, m n = 135MeV, x = 16, and quark masses m = 229MeV and 264MeV for the NJL or 
ENJL cases as listed in Table |. The contribution from — \n.(m\/ p 2 ) with p = 2m is displayed as 
an intermediate step in both cases. 



Parameter 



Empirical value, Ref. 



NJL 


5.0 + 2.4 = 


= 7.4 


3.5 + 2.4 = 


= 5.9 


7.0 + 2.4 = 


= 9.4 


13 + 2.4 = 


15.4 


1.4 x 10 


-3 



ENJL 



h 
k 
h 
h 
h 



2.9 + 2.4 
4.3 + 0.9 
13.9 ± 1.3 
16.5 ± 1.1 

~5x 10~ 3 



4.5 + 2.7 = 7.2 
1.5 + 2.7 = 4.2 
10.1 + 2.7= 12.8 
13.5 + 2.7= 16.2 
0.5 x 10~ 3 



TABLE IV. Comparison of second sheet real and complex poles of the p and a\ propagators in 
units of 4m 2 . The type of dispersion relation used to generate the polarization amplitude employed 
in Eqs. (|3.22j ) and fl3.25| ) is indicated on the left. The ghost poles only occur if subtracted dispersion 
relations are used, and lie on the first, or physical sheet. The input parameters are the standard 
set as given in the second line of Table |[ 

p pole p pole cti pole 

Unsubtr. 0.1658 0.5221 - 1.3260t 2.1227 - 2.1405i 

Subtr. 0.1659 0.5557 - 1.3473i 2.2340 - 2.1728i 

Ghost poles - -18.3038 -18.0418 



TABLE V. Comparison of the second sheet complex poles of the p and a\ propagators and 
the peak positions in their spectral densities. The complex poles are recorded in the usual form 
s/M = M — iT in order to display a mass and a width contribution. We have used a capital letter 
for the mass that appears as the real part of the complex pole in order to distinguish it from the 
lower case m p and m ai that we associate with the position of the peak value of the corresponding 



spectral function. Parameters are as for Table IV. The experimental p(770) and ai(1260) meson 



masses are given in brackets. The energy unit is MeV. 





p-mode 


ai-mode 


Pole position 


382 - 970i 


770 - 777i 


Peak position 


713 (768+0.6) 


1027 (1230+40) 
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TABLE VI. Comparison of the asymptotic behavior as \s\ — > oo in the complex s-plane of 
Pauli-Villars regulated Feynman loop diagrams, employing either unsubtracted or once-subtracted 
dispersion relations. The amplitudes that are generated by the first method carry a bar in the 



main text. 


AmnHt.nHp 


T In fin Irtvrf with Har^ 


Suhtrf without hart 




1 /c 

i/s 


ln{-s) 


Tl pp (<j) 
11 [hj 


l^i/ft; 111^ a; 


— e In 1 — 


u)f v ( s ) ~ n^ A (s) 


-(l/s)ln(-s) 


s ln(— s) 


F P (s) 




-ln(-s) 


F v (s) ~ G A (s) 


1 


l/(s ln(-a)) 


F P (s)Fy(*)<3U(*) 




-l/0 2 ln(-s)) 


(0)F/ \ {0)A, \ 




s 




-1/s 


l/s(ln|s|) 2 


pr( S )-pS o)A ( 5 ) 


-1/s 2 


positive constant 




-1/s 2 


-l/(s In |s|) 2 


TABLE VII. Values of isospin traces 


Trace of 


7T + 7T 






4 

9 ! 


1 5 
! 9 


1 f 2 ±l 2 


1 1 1 

2 2 2 


J- fZ-q-l-i 2 


1 ; 

3 i 


! 1 
1 2 


1 / 2 J * e 9 


2 1 1 

3 3 2 




8 I 

9 < 


! 5 
) 9 


T ] Tp 2 
1 f ±ie q 


2 J 
9 < 


i 5 
) 9 


T^p T 3T. 
J- ft q 2 -tj 


2 1 1 

3 3 2 


J- f-Li^q 2 


1 I 

3 5 


! 1 

! 2 




1 1 1 

2 2 2 
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